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The body of research on (III,Mn)V diluted magnetic semiconductors initiated during the 1990's 
has concentrated on three major fronts: i) the microscopic origins and fundamental physics of 
the ferromagnetism that occurs in these systems, ii) the materials science of growth and defects 
and iii) the development of spintronic devices with new functionalities. This article reviews the 
current status of the field, concentrating on the first two, more mature research directions. From 
the fundamental point of view, (Ga,Mn)As and several other (III,Mn)V DMSs are now regarded 
as textbook examples of a rare class of robust ferromagnets with dilute magnetic moments cou- 
pled by delocalizod charge carriers. Both local moments and itinerant holes are provided by Mn, 
which makes the systems particularly favorable for realizing this unusual ordered state. Advances 
in growth and post-growth treatment techniques have played a central role in the field, often 
pushing the limits of dilute Mn moment densities and the uniformity and purity of materials far 
beyond those allowed by equilibrium thermodynamics. In (III,Mn)V compounds, material quality 
and magnetic properties are intimately connected. In the review we focus on the theoretical un- 
derstanding of the origins of ferromagnetism and basic structural, magnetic, magneto-transport, 
and magneto-optical characteristics of simple (III,Mn)V epilayers, with the main emphasis on 
(Ga,Mn)As. The conclusions we arrive at are based on an extensive literature covering results of 
complementary ab initio and effective Hamiltonian computational techniques, and on comparisons 
between theory and experiment. The applicability of ferromagnetic semiconductors in microelec- 
tronic technologies requires increasing Curie temperatures from the current record of 173 K in 
(Ga,Mn)As epilayers to beyond room temperature. The issue of whether or not this is a realistic 
expectation for (III,Mn)V DMSs is a central question in the field and motivates many of the 
analyses presented in this review. 
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conduction-band electrons or valence-band holes. This 
definition implies that in ferromagnetic semiconductors, 
magnetic properties can be influenced by the same assort- 
ment of engineering variables that are available for other 
more conventional semiconductor electronic properties. 
In the best understood arsenide DMSs, semiconductor 
valence-band carriers participate in the magnetic order, 
-rr: : — I ,^ The materials require participation of valence band holes 

DiscussionI 46 ^ . . ^ ^. „ „ 

Ia. Magnetic interactions in systems with coupled local and itiifeihl^f'iAjMfl lg" of a ferromagnetic State. Efforts to m- 

crease their critical temperatures further run into incom- 

!B. Ferromagnctism in the diluted magnetic semiconductor fanijt^fe flely understood fundamental limits On the ratio of the 

magnetic transition temperature to the Fermi tempera- 
ture of the free-carrier systems and the role of disorder 
in these heavily doped materials. The tension between 
achieving high Curie temperatures and the desire for low, 
and therefore gateable, carrier densities, is among the 
major issues in the study of these materials. 

In this article we review the considerable theoretical 
progress that has been made in understanding the very 
broad range of properties that occur in (III,Mn)V fer- 
romagnetic semiconductor epilayers in diS'erent regimes 
of Mn content and defect density. The main focus of 
this article is on the extensively studied (Ga,Mn)As fer- 
romagnetic semiconductor, but we also make frequent 
comments on other (III,Mn)V DMSs. Comparisons to 
experimental data are made throughout the article. In 
Section I we review progress that has been achieved in 
the effort to realize useful DMS materials for spintron- 
ics (or magnetoelectronics) . In Section II we discuss the 
properties of dilute Mn atoms in a (III,V) crystal, and 
the various mechanisms that can couple the orientations 
of distinct moments and lead to ferromagnctism. In Sec- 
tion III we discuss several different strategies that can be 
used to elevate material modeling from a qualitative to a 
more quantitative level. Sections IV- VII address a vari- 
ety of different characteristics of (III,Mn)V layers, includ- 
ing their structural, magnetic, magneto-transport, and 
magneto-optical properties. Finally in Section VIII we 
discuss the ferromagnetic ordering physics in (III,Mn)V 
DMSs in the broad context of magnetic interactions in 
systems with coupled local and itinerant moments, and 
then extrapolate from (III,Mn)V materials to comment 
on the effort to find high temperature ferromagnctism in 
other DMS materials. We conclude in Section IX with a 
brief summary. 

To partially remedy omissions in the bibliography that 
originate from our incomplete coverage of this topic, we 
refer to an extended database of published work and 
preprints maintained at http://unixl2.fzu.cz/ms. The 
structure of the database is similar to the structure of 
this review and we encourage the reader in need of a 
more detailed bibliography to use this resource. 

A number of review articles on various aspects of the 
physics of DMSs have been published previously and 
may help the reader who seeks a broader scope than we 
are able to supply in this review. The extensive body 
of research on DMSs in the 198 0's, fo c used mostly on 



Semiconductor physics and magnetism are established 
subfields of condensed matter physics that continue to 
reveal a rich variety of new phenomena, often in new 
types of solid state materials. The properties of semicon- 
ductors are extraordinarily sensitive to impurity atoms, 
defects, and charges on external gates. Magnetism is a 
collective electronic phenomenon with an ordered state 
that is often stable to exceptionally high temperatures. 
Magnetic order, when it is present, has a large impact 
on other material properties including transport and op- 
tical properties. In both semiconductor and magnetic 
cases, sophisticated and economically important tech- 
nologies have been developed to exploit the unique elec- 
tronic properties, mainly for information processing in 
the case of semiconductors and for information storage 
and retrieval in the case of magnetism. 

The realization of materials that combine semicon- 
ducting behavior with robust magnetism has long been 
a dream of material physics. One strategy for creat- 
ing systems that arc simultaneously se miconducting and 
magnetic , initia t ed in the late 1970's ijGai et all Il978t 
Uaczvnski et all Il978|) . is to introduce local moments 
into well-understood semiconductors. The result is a new 
class of materials now known as diluted magnetic semi- 
conductors (DMSs). Over the past fifteen years, build- 
ing on a series of pioneering publications in the 1990's 
llHavashi et a?J.ll997tlMunekata et a?J.ll989llT993lOhnol 
119981: lOhno et adll992lll996bHVan Esch et oi.l.ll997^ ■ it 
has been established that several (III,V) compound semi- 
conductors become ferromagnetic when heavily doped 
with Mn, and that the ferromagnetic transition tem- 
peratures can be well above 100 K. In semiconductors 
like GaAs and InAs, Mn has been shown to act both as 
an acceptor and as a source of local moments. These 
(III,Mn)V materials are examples of ferromagnetic semi- 
conductors, a phrase we reserve for magnetic systems in 
which ferromagnctism is due primarily to coupling be- 
tween magnetic clement moments that is mediated by 



(II,Mn)VI alloys, is reviewed in 



iFurdvnal 
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Il988t iFurdv na and Kossutl Il988|) . Several extended pa- 
pers cover the experimental properties of (III,Mn)V 
DMSs, particularly (Ga,Mn)As and (In,Mn)As, in- 
terpret ed within the car rier-medi ated ferromagnetism 
model (iMacDonald et all 12005; .Matsukura et all l2002t 
lOhnoi I1999I) . Theoretical predictions based on this 
model for a number of properties of b ulk DMSs and 
heterostructures are reviewed in dPiet l. 2002, 200i 
iKonig et all l2003t iLee et all l2002|) . A detailed de- 



found in 


Fukumura et aZ.l.l2005 




Graf et 


a/J.l2003hl 


iLiu et al. 


. l2005l IPearton et al. 


. I2OO.I 


). We 


also men- 



tion here several specialized theoretical reviews focus- 
ing on the predictions of density functi onal first prin- 
ciples calculations for (III,Mn)V DMSs ijSanvito et all 
l2002HSato and Katavama-Yoshidal |2002|) . on Mn-doped 
II- VI and III-V DMSs in the low carrier density 
regime (iBhatt et ffld | 200 2D. and on effects of disorder in 
fGa.Mn) As ^Timml . 12003^) . 



A. Functional (lll,Mn)V material requirements 

III-V materials are among the most widely used semi- 
conductors. There is little doubt that ferromagnetism in 
these materials would enable a host of new microelectron- 
ics device applications if the following criteria were met: 
i) the ferromagnetic transition temperature should safely 
exceed room temperature, ii) the mobile charge carriers 
should respond strongly to changes in the ordered mag- 
netic state, and iii) the material should retain fundamen- 
tal semiconductor characteristics, including sensitivity to 
doping, electric fields produced by gate charges, and sen- 
sitivity to light. For more than a decade these three 
key issues have been the focus of intense experimental, 
and theoretical research into the material properties of 
Mn-doped III-V compounds. At first sight, fundamen- 
tal obstacles appear to make the simultaneous achieve- 
ment of these objectives unlikely. Nevertheless, inter- 
est in this quest remains high because of the surprising 
progress that has been achieved. Highlights of this sci- 
entific endeavor are briefly reviewed in this introductory 
chapter. 



B. Search for high transition temperatures 

Under equilibrium growth conditions the incorporation 
of magnetic Mn ions into III- As semiconductor crystals is 
limited to approximately 0.1%. Beyond this doping level, 
surface segregation and phase separation occur. To cir- 
cumvent the solubility problem a non-equilibrium, low- 
temperature molecular-beam-epitaxy (LT-MBE) tech- 
nique was applied and led to first successful growth 
of (In,Mn)As and (Ga,Mn)As DMS ternary alloys with 
more than 1% Mn. Since the first report in 1992 of a 
ferromagnetic transition in p-type (In,Mn)As a t a criti- 
cal temperature Tc = 7.5 K (jOhno et a^J . ll99'^ . the un- 



folding story of critical temperature limits in (III,Mn)V 
DMSs has gone through different stages. Initial exper- 
iments in (In,Mn)As suggested an intimate relation be- 
tween the ferromagnetic transition and carrier localiza- 
tion, reminiscent of the behavior of manganites (per- 
ovskite (La,A)Mn03 with A=Ca, Sr, or Ba) in which 
ferromagnetism arises from a Zcner double exchange 
proc ess associated with d-electron hopping between Mn 
ions (jCoev et a/J . ll999(l . (We comment at greater depth 
on qualitative pictures of the ferromagnetic coupling 
in Sections III.AI and IVIILAI ) This scenario was cor- 
roborated by a pioneering t heoret i cal ab -initio study of 
the (In,Mn)As ferromagnet ijAkail . Il998j) and the mech- 
anism was also held responsible for mediating ferro- 
magnetic Mn-Mn coupling in some of the first fcrro- 
magnetic (Ga,Mn)As samples with T^s close to 50 K 
(|Van Esch adll997|) . 

In 1998 the Tohoku University group annou nced a 
jump of Tc in p-type (Ga,Mn)As to 110 K ijOhnol 
I1998D and pointed out that the critical temperature 
value was consistent with the kinetic-exchange mecha- 
nism for ferromagnetic coupling, also first proposed by 
Zener (see Section Hi. A() . In its simplest form, ferr omag- 
netism in this picture follows l)Dietl et all Il997t) from 
Ruderman-Kittel-Kasuya-Yosida (RKKY) indirect cou- 
pling between Mn d-shell moments mediated by induced 
spin-polarization in a free-hole itinerant carrier system 
(sec Sections III.AI IV.A.2I and IVIII.A|) . Zener proposed 
this mechanism originally for transition metal ferromag- 
nets for which the applicability of this picture is now 
known to be doubtful because of the itinerant character of 
transition metal d-electrons. The model of Mn((i^) local 
moments that are exchange-coupled to itinerant sp-band 
carriers does however provide a good desc ription of Mn- 
doped IV- VI and n-VI DMSs l|Dietlll994|) . The key dif- 
ference between (III,Mn)V materials like (Ga,Mn)As and 
IV- VI and II- VI DMSs is that Mn substituting for the 
trivalent cation (Ga) is simultaneously an acceptor and a 
source of magnetic moments. Theoretical critical temper- 
ature calculations based on the kinetic-exchange model 
predict room temperature ferromagnetism in (Ga,Mn)As 
with 10% Mn content. In spite of these optimistic predic- 
tions, the goal of breaking the 110 K record in (Ga,Mn)As 
remained elusive for nearly four years. Only recently 
has progress in MBE growth and i n the development of 
post-growth annealing techniques JCliib^^^Zjj_|2003£ ; 
[Edmonds et al gOOga , Eid et aA l2005t iHavashi et al. , 
boOlTlKu et~d\ . l2003t IYu et allmO^ made it possible 
to suppres s extrinsic effects, pushing Tr in (Ga,Mn)As u p 
to 173 K ijJungwirth et aZ.I. i2005bt IWang adl2005a;i . 
Tc trends in current high quality (Ga,Mn)As epilayers 
are consistent w i th the Zener kinetic-exchange model 
ijjungwirth et alj l2005b() . The current Tc record should 
be broken if DMS material with a higher concentration 
of substitutional Mn ions can be grown. 

Based on the few experimental and theoretical stud- 
ies reported to date, (III,Mn)Sb DMSs are expected 
to fall into the same category as (Ga,Mn)As and 
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(In,Mn)As DMS. The kinetic- exchange model calcula- 
tions predict T^s th at are smal l compared to their ar- 
senide counterparts l|Dietl et all 1200(1 I.Tungwirth et all 
l2002bt). This difference, confirmed by experiment 
llAbe all l2000t ICsontos et all l2005t iPanguluri et all 
l2004HWoitowicz et a/.l . l200,'^ . is caused by the weaker p— 
d exchange and smaller magnetic susceptibility (smaller 
effective mass) of itinerant holes in the larger unit-cell 
antimonides. Also consistent with the kinetic- exchange 
model is the remarkable obseryation of an increase of 
in (In,M n)Sb by 25% induced by the applied hydrostatic 
pressure (|Csontos et a?J . l2005|) . 

Moying in the opposite direction in the periodic ta- 
ble towards (III,Mn)P and (III,Mn)N appears to be the 
natural route to high ferromagnetic semiconductors. 
The kinetic-exchange model predicts T^s far above room 
temperature in these small er lattice constant materials, 
in particular in (Ga,Mn)N (|Dietl et allhOOCt) . Also the 
solubility limit of Mn is much larger than in arsenides, 
making it possible in principle to grow highly Mn-doped 
DMSs under or close to equilibrium conditions. How- 
ever, the nature of the magnetic interactions in Mn- 
doped phosphides and nitrides is not complet ely under- 
stood either theoretically or experimentally ()Liu et all 
|2005|) . As the valence band edge moves closer to the 
Mn d-level and the p — d hybridization increases with 
increasing semiconductor gap width and decreasing lat- 
tice constan t, charge fluctuations of the d-states may be- 
come large (ISa ndrats kii et all l2004t ISanval et ail l2003t 
IWierzbowska et al., . .2004t) . With increasing ionicity of 
the host crystal, the Mn impurity may also undergo a 
transition from a divalent acceptor to a trivalent 
neutral impurity (,KrcissLe,it al._, 1996t iLuo and Martini 
l2005HSchulthess et all\200^ . In either case, the picture 
of ferromagnetism based on the Zener kinetic- exchange 
model needs to be revisited in these materials. 

Experimental critical temperatures close to 1000 K 
have been r eported in some (Ga.Mn)N samples 
l|Sasaki et all l2002a|) . It is still unclear, how- 
ever, whether the high-temperature ferromagnetic phase 
should be attributed to a (Ga,Mn)N ternary alloy or to 
the presence of ferromagnetic metal precipitates embed- 
ded in the host GaN lattice. Reports of (Ga,Mn)N epi- 
layers synthesized in cubic and hexagonal crystal struc- 
tures, of p-typc and n-type ferromagnetic (Ga,Mn)N, 
and of multiple ferromagnetic phases in one material 
all add to t he complex p henomenology of these wide- 
gap DMSs iArkun et aJ] l2004 Isdmonds et all l2004c[ 
2005bt IC^f et all l2002l l2003at [ Hwang et all 1200,^ 
Korotkov et a^J . l200lt ISawicki et g/J . l2005a|) . 



Uncertainties apply also to the interpretation of fer- 
romagnetism seen in the (Ga,Mn)P samples studied 
to date, which have been prepared by post MBE 
ion-implantati o n of Mn followed by rapi d ther mal 
l|Poddar et all l2005t iTheodoropoulou et all l2002l) o r 
pulse laser melting annealing ijScarpulla et al I 1200,^ . 
Experiments in these materials have not yet established 
unambiguously the nature of magnetic interactions in the 



(III,Mn)P compounds. However, a comparative study of 
(Ga,Mn)P and (Ga,Mn)As prepared by the post MBE 
ion-implantation and pulse laser melting annealing sug- 
gests carrier mediat ed origin of ferromag netism in the 
(Ga,Mn)P material ijScarpuUa et all\2003l} . 

Our current understanding of the material physics of 
(III,Mn)V DMS epilayers suggests that synthesis of room 
temperature ferromagnetic semiconductor will require a 
level of doping and defect control comparable to what has 
now been achieved in high quality (Ga,Mn)As samples, 
Mn densities of order 10%, and may require the use of 
wider gap III-V alloys. 

Finally, we note that efforts to enhance Curie tem- 
perature in Mn-doped (III,V) semiconductors have also 
led to material research in more complex semicon- 
ductor heterostructures with highly Mn-doped mono- 

showing promising r esults 



layers (i5-doped layers), 
(IChen et all '2002*. 'Fcrna 

mi et QL.200Q: A-h 

20031 ISanvitol l2003t IVurgaftman and Mevei 



Fernandez- Rossier and Sham, 200J ; 
Kawakami et Q^j.|200ll: .MYcrs a?.,.^2004^:,Nazmul et al. 



200E, 



2001) 



C. Conventional spintronics 

Spintronic devices exploit the electron spin to manip- 
ulate the flow of electrons and therefore require ma- 
terials in which the charge and spin degrees of free- 
dom of carriers a re strongly coupled llDe Boeck et all 
l2002t IWolf et all l200lt IZutic et all l2004^ . The most 
robust, and currently the most useful, spintronic de- 
vices rely on the collective behavior of many spins in 
ferromagnetic materials to amplify the coupling of ex- 
ternal magnetic fields to electronic spins, a coupling 
that is very weak for individual electrons. The in- 
trinsically large spin-orbit interaction in Ill-Sb and 
III-As valence band states makes these hosts ideal 
candidates for exploring various spintronic functionali- 
ties. In (Ga,Mn)As DMS epilayers, for example, the 
measured anisotropic magneto-resistance (AMR) effect 
(the relative difference between longitudinal resistivi- 



~ 10% (iBaxter et all 


20021 iGoennenwein et al. 


Jungwirth et all l2002a 


. I2003bl iMatsukura et al 


Tang et aZ.I. l2003HWang et a?J.l2005d 


1200211. 





2005 




2004 



A particularly strong manifestation of valence-band 
spin-orbit coupling occurs in the antisymmetric off- 
diagonal element of the resistivity tensor. The anoma- 
lous Hall effect (AHE) shown in Fig.^ which completely 
dominates the low-field Hall response in (Ga,Mn)As and 
some other III-V DMSs, has become one of the key tools 
used to detect the paramagn e tic /fe rromagnetic transi- 
tion l)Ohnd. Il998t lOhno et all Il992(l . Its large value is 
due to the spin-polarization of holes and provides strong 
evidence for the participation of mobile charge carriers 
in the ordered magnetic state of these DMSs. 

In metals, the current response to changes in the 
magnetic state is strongly enhanced in layered struc- 
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0.0 
B(T) 

FIG. 1 Hall resistance versus external magnetic field for an- 
nealed Gao.94Mno.06 As at 110 K (black), 1 30 K (red), and 
140K (green). From ijEdmonds et a^.Ll2002a^ . 



tures consisting of alternating ferromagnetic and non- 
magnetic materials. The giant magneto- resistance effect 
()Baibich et alL\l98^ which is widely exploited in current 
technology, for example in field-sensors and magnetic 
random-access memories, reflects the large difference be- 
tween resistivities in configurations with parallel and an- 
tiparallel polarizations of ferromagnetic layers in mag- 
netic superlattices or trilayers like spin- v alves and mag- 
netic tunnel junctions l|G!regg et a/J . 12003) . The effect re- 
lies on transporting spin information between layers and 
therefore is sensitive to spin-coherence times in the sys- 
tem. Despite strong spin-orbit coupling which reduces 
spin coherence in DMSs, functional spintronic trilayer de- 
vices can be built, as demonstrated by the measured large 
MR effects in (Ga,Mn)As based tunneling structures 
llChiba et all l2004at iMattana et all l2005l: ISaito et all 
l2005t iTanaka and Higd . l20nit) . The coercivities of indi- 
vidual DMS layers can be tuned via exchange-biasing to 
an antiferromagnet (|Eid et fflZ.l . l200l which is a standard 
technique used in m etal giant-magnetoresistance devices 
jCreiiieranilQQl). 



D. Magneto-semiconducting properties and related new 
spintronics effects 

DMS ferromagnets possess all the properties that are 
exploited in conventional spintronics. They qualify as fer- 
romagnetic semiconductors to the extent that their mag- 
netic and other properties can be altered by the usual 
semiconductor electronics engineering variables. The 
achievement of ferromagnetism in an ordinary III-V semi- 
conductor that includes several per cent of Mn demon- 
strates on its own the sensitivity of magnetic properties 
to doping. Remarkably, doping profiles and, correspond- 
ingly magnetic properties can be grossly changed, even 
after growth, by annealing. Early studies of (Ga,Mn)As 
indicated that annealing at temperatures above the 
growth temperature leads to a reduction of magnetically 
and electrically active Mn ions and, at high enough an- 
nealing temperatures, to the formation of MnAs clusters 
ijVan Esch et a/J . ll997(l . On the other hand, annealing at 




-150 -100 -50 50 100 150 

Magnetic Field (mT) 

FIG. 2 Top panel: False-color SEM picture (side-view) of a 
double constriction showing part of the outer wires with the 
voltage leads. The insets show the relative magnetization of 
the different parts (left) and the resulting schematic magne- 
toresistance trace for sweep-up (solid line) and sweep-down 
(dashed line). Bottom panel: Measured magnetoresistance 
in a sample with tun nel barriers at the constrictions. From 
llRuester et a;.ll2003D . 



temperatures below the growth temperature can substan- 
tially improve magnetic and transport properties of the 
thin DMS layers due to the out-diffusion of charge and 
moment co mpensating defects, now identified as inter- 
stitial Mn JChiba et all l2003at lEdmonds all l2002at 
lEid et oil 1200,4 iKu et alU2m?k\Yn'et a/J . l2002D . 

(In,Mn)As based field effect transistors were built to 
study electric field control of ferromagnetism in DMSs. 
It has been demonstrated that changes in the carrier 
density and distribution in thin film DMS systems due 
to an applied bias voltage can revers ibly induce the fer - 
romagnetic/paramagnetic transition l)Ohno et aZ.U200 0|). 
Another remarkable effect observed in this magnetic 
transistor is electric field assisted magnetization reversal 
l|Chiba et all l2003bj) . This novel functionality is based 
on the dependence of the width of the hysteresis loop on 
bias voltage, again through the modified charge density 
profile in the ferromagnetic semiconductor thin film. 

Experiments in which ferromagnetism in a (III,Mn)V 
DMS system is turned on and off optically add to the 
list of functionalities that result from the realization of 
carrier-in duced ferromagnetism i n a semiconductor hos t 
material ijKoshihara et a;J . ll997tlMunekata et a/J . ll997(l . 
The observed emission of circularly polarized light from a 
semiconductor heterostructure, in which electrons (holes) 
injected from one side of the structure recombine with 
spin-p olarized holes (electrons) emitted f rom a DMS 
layer l)FiederUng et all Il999t lOhno et all \199^ . is an 
example of phenomena that may lead to novel magneto- 
optics applications. 
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Tunneling anisotropic magnetoresistance (TAMR) is 
another nov e l spin t ronic effect o b served in (Ga,Mn)As 
(iRrev Pi all IM lOonId et all HqHI m.nester et dl 
l2005HSaito et alVhm^ . TAMR, like AMR, arises from 
spin-orbit coupling and reflects the dependence of the 
tunneling density of states of the ferromagnetic layer on 
the orientation of the magnetization with respect to the 
current direction or the crystallographic axes. 

The larger characteristic electronic length scales in 
DMSs compared to ferromagnetic metals make it pos- 
sible to lithographically define lateral structures with 
independent magnetic areas coupled through depleted 
regions that act as tunnel barriers and magnetic weak 
links. The electrical response to magnetization rever- 
sals in these spintronic nanodevices can lead to MR ef- 
fects with magnitudes of order 1000% llRiiester et al\ . 
l200ril). as shown i n Fig . El and with a rich phenomenology 
ijGiddings et alV l2005|) . Wider lateral constrictions have 
been used to demonstrate controlled do main-wall nucle- 
ation and propagation in DMS stripes ijHonolka et alV 
l2005t iRuester et al\ . l2003j) . a prerequisite for develop- 
ing semicondu ctor logic gates based on magnetic domain 
mani pulation ijAUwood et alV l2005t iGate and Reeisteil 
|2002() . (Ga,Mn)As nanoconstrictions with lateral side 
gates have revealed a new effect. Coulomb blockade 
anisotropic magnetoresistance, which reflects the mag- 
netization orient ation dependence of the single-electron 
charging energy ijWunderlich et al\. |2006() . These spin- 
tronic single-electron transistors offer a route to non- 
volatile, low-field, and highly electro- and magneto- 
sensitive operation. 



II. THE ORIGIN OF FERROMAGNETISM 

A. General remarks 

The magnetic dipole-dipole interaction strength be- 
tween two discrete moments separated by a lattice con- 
stant in a typical solid is only ^ 1 Kelvin, relegating di- 
rect magnetic interactions to a minor role in the physics 
of condensed matter magnetic order. Relativistic effects 
that lead to spin-orbit coupling terms in the Hamiltonian 
provide more plausible source of phenomena that are po- 
tentially useful for spintronics. Although these terms are 
critical for specific properties like magnetic anisotropy, 
they are rarely if ever crucial for the onset of the magnetic 
order itself. Instead the universal ultimate origin of fer- 
romagnetism is almost always the interplay between the 
electronic spin degree of freedom, repulsive Coulomb in- 
teractions between electrons, and the fermionic quantum 
statistics of electrons. The Pauli exclusion principle cor- 
relates spin and orbital parts of the many-electron wave- 
function by requiring the total wavefunction to be anti- 
symmetric under particle exchange. Whenever groups of 
electrons share the same spin state, the orbital part of the 
many-body wavefunction is locally antisymmetric, low- 
ering the probability of finding electrons close together 



and hence the interaction energy of the system. Because 
magnetic order is associated with the strong repulsive 
Coulomb interactions between electrons, it can persist to 
very high temperatures, often to temperatures compara- 
ble to those at which crystalline order occurs. Ferromag- 
netism can be as strong as chemical bonds. Very often 
the quantum ground state of a many electron system has 
non-zero local spin-density, either aligned in the same di- 
rection in space at every point in the system as in simple 
ferromagnets, or in non-coUinear, ferrimagnetic, or an- 
tiferromagnetic materials in configurations in which the 
spin direction varies spatially. 

Although this statement on the origin of magnetic 
order has very general validity, its consequences for a 
system of nuclei with a particular spatial arrangement, 
are extraordinarily difficult to judge. Because ferromag- 
netism is a strong-coupling phenomenon, rigorous the- 
oretical analyses are usually not possible. There is no 
useful universal theory of magnetism. Understanding 
magnetic order in a particular system or class of sys- 
tems can be amo ng the most ch allenging of solid state 
physi cs problems ijAshcroft and M ermin. 1976; MardeJ. 
119991) . For most systems, it is necessary to proceed in 
a partially phenomenological way, by identifying the lo- 
cal spins that order, and determining the magnitude and 
sign of the exchange interactions that couple them by 
comparing the properties of simplified (often 'spin-only') 
model Hamiltonians with experimental observations. 

One approach that is totally free of phenomenologi- 
cal parameters is density functional theory (DFT), in- 
cluding its spin-density functional (SDF) generalizations 
in which energy functionals depend on charge and spin 
densities. Although DFT theory is exact in principle, 
its application requires that the formalism's exchange- 
correlation energy functional be approximated. Approx- 
imate forms for this functional can be partially phe- 
nomenological (making a pragmatic retreat from the 
ah initio aspiration of this approach) and arc normally 
based in part on microscopic calculations of correlation 
effects in the electron gas model system. This is the 
case for the oft en-used local (spin) density approxima- 
tion (L(S)DA) l|von Barth and HedinL Il972t) . For many 
magnetic metals, in which correlations are somewhat sim- 
ilar to those in the electron gas model system, ah ini- 
tio LSDA theory provides a practical and sufficiently 
accurate solution of the mag netic many-body problem 
ijjones and Gunnarssonl Il989j) . This is particularly true 
for elemental transi tion metal ferromagnets Fe, Co, and 
Ni an d their alloys ijMarderl Il999t iMoruzzi and Marcusl 
Il993() . In practice LSDA theory functions as a mean- 
field theory in which the exchange-energy at each point in 
space increases with the self-consistently determined lo- 
cal spin density. With increasing computer power, LSDA 
theory has been applied to more complex materials, in- 
cluding DMSs. 

As we discuss below, both phenomenological and DFT 
approaches provide valuable insight into (III,Mn)V fcr- 
romagnetism. Model Hamiltonian theories are likely 
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to remain indispensable because, when applicable, they 
provide more transparent physical pictures of ferromag- 
netism and often enable predictions of thermodynamic, 
transport, and other properties that are sometimes (de- 
pending on material complexity), beyond the reach of 
ab initio theory techniques. Of particular importance 
for DMSs is the capability of model Hamiltonians to 
describe localized electronic levels coincident with an 
itinerant-electron band which, strictly speaking, is be- 
yond the reach of the effectively one-particle band theo- 
ries of solids that emerge f r om L SDA theory l)Andersonl 
Il96lt ISchrieffer and Wolfl Il96fil^ . Interpreting experi- 
ments with model Hamiltonian approaches can, on the 
other hand, be misleading if the model is too simplified 
and important aspects of the physics are absent from 
the model. What's more, even simplified models usu- 
ally leave complex many-body problems that cannot be 
completely solved. 

Magnetism in (Ga,Mn)As and some other (III,Mn)V 
ferromagnets originates from the Mn local moments. (As 
already pointed out, Mn-doped phosphides and nitrides 
DMSs are less well understood, however, local Mn mo- 
ments are likely to play an important role in these materi- 
als as well.) The dependence of the energy of the system 
on the relative orientation of Mn moments is generally 
referred to as an exchange interaction. This terminology 
is part of the jargon of magnetism and recognizes that 
Fermi statistics is the ultimate origin. Several types of 
qualitative effects that lead to exchange interactions can 
be separately identified when addressing magnetic order 
in (III,Mn)V DMSs; the applicability of each and the 
relative importance of different effects may depend on 
the doping regime and on the host semiconductor mate- 
rial. In this section we first introduce the jargon that is 
commonly used in the magnetism literature, by brieffy 
reviewing some of the effects that can lead to magnetic 
coupling, most of which have been recognized since near 
the dawn of quantum mechanics. 

For spins carried by itinerant electrons, exchange inter- 
actions are often most simply viewed from a momentum 
space rather than a real sp ace point of vi ew. S toner's 
itinerant- exchange ijAshcroft and Mermlnl . Il976j) favors 
spontaneous spin-polarization of the entire electron gas 
because electrons are less likely to be close together and 
have strongly repulsive interactions when they are more 
likely to have the same spin. Because the band energy is 
minimized by double-occupation of each Bloch state, the 
Stoner ferromagnetic instability occurs in systems with 
a large density of states at the Fermi energy. This helps 
to explain, for example, why fcrromagnetism occurs in 
the late 3d transition elements. A large density of states 
makes it possible to gain exchange energy by moving elec- 
trons from one spin-band to the other while keeping the 
kinetic energy cost sufficiently low. Since the key spins in 
many (III,Mn)V DMS materials are localized the Stoner 
mechanism does not drive fcrromagnetism, although we 
will see in Section IV. Al that it still plays a minor sup- 
porting role. 



In many systems, including (III,Mn)V DMSs, both the 
local nature of the moments and strong local Coulomb 
interactions that suppress charge (valence) fluctuations 
play a key role and have to be included even in a minimal 
model. Many mechanisms have been identified that cou- 
ple localized spin s in a solid. The origin of He isenberg's 
direct-exchange ijAshcroft and Merminl . Il976(l between 
two local spins is the difference between the Coulomb en- 
ergy of a symmetric orbital wavefunction (antisymmetric 
singlet spin wavefunction) state and an antisymmetric or- 
bital wavefunction (symmetric triplet spin wa vefunction) 
state . Kramer's super- exchange interaction jAndersonl 
Il950() . applies to local moments that are separated by a 
non-magnetic atom. In a crystal environment, an elec- 
tron can be transferred from the non-magnetic atom to 
an empty shell of the magnetic atom and interact, via di- 
rect exchange, with electrons forming its local moment. 
The non-magnetic atom is polarized and is coupled via 
direct-exchange with all its magnetic neighbors. Whether 
the resulting super-exchange interaction between the lo- 
cal moments is ferromagnetic or antiferromagnetic de- 
pends on th e relative sign of the two dire c t-exch ange 
interactions ^Goodenoughl Il958l: iKanamoril Il959(l . In 
(III,Mn)V materials, superexchange gives an antiferro- 
magnetic contribution to the interaction between Mn mo- 
ments located on neighboring cation sites. 

Zener's double-exchange mechanism ijZeneil Il951bj) 
also assumes an intermediate non-magnetic atom. In 
its usual form, this interaction occurs when the two iso- 
lated magnetic atoms have a different number of electrons 
in the magnetic shell and hopping through the inter- 
mediate non-magnetic atom involves magnetic-shell elec- 
trons. Combined with the on-shell Hund's rule, double- 
exchange couples the magnetic moments ferromagneti- 
cally. Parallel spin alignment is favored because it in- 
creases the hopping probability and therefore decreases 
the kinetic energy of the spin-polarized electrons. A ver- 
sion of double-exchange, in which Mn acceptor states 
form an impurity band with a mixed spd-character, has 
often been referred to in the (III,Mn)V literature. In this 
picture electrical conduction and Mn-Mn exchange cou- 
pling are both realized through hopping within an impu- 
rity band. The potential importance of double-exchange 
is greater at lower Mn doping and in wider-gap (III,Mn)V 
materials. 

Fina lly, we identify Zener's kinetic- exchange ijZeneil 
Il951a^ or indirect exchange interaction. It arises in mod- 
els with local, usually d-shcU or /-shell, moments whose 
coupling is mediated by s- or p-band itinerant carriers. 
The local moments can have either a ferromagnetic direct 
exchange interaction with band electrons on the same 
site and/or an antiferromagnetic interaction due to hy- 
bridization between th e local moment and band el e ctrons 
on n eighboring sites jBhattachariee et Ml, Il98,'^ iDieti 
Il994|) . Polarization of band electrons due to the inter- 
action at one site is propagated to neighboring sites. 
When the coupling is weak (the band carrier polariza- 
tion is weak, e.g., at temperatures near Curie tempera- 
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ture) , the effect is described by RKKY theory which was 
originally applied to carrier mediated indirect coupling 
betw'een nuclear moments |BlQemla£r££iL.aiid-Ji<Iwian(: , 
19551 'Frohlich anc l Nabarrolll940HRuderman and Kitte , 
1954: .Yosidat. .1957|) and between local d-shell moments 
in meta,ls (|Ka,si" ll lYosidal . Il957t IZenerl Il951al] . 
The range of this interaction can be long and the in- 
teractions between separate local moments can be either 
ferromagnetic or antiferromagnetic and tend to vary in 
space on the length scale of the itinerant band's Fermi 
wavelength. Unlike the double-exchange case, magnetic 
order in this case does not lead to a significant change 
in the width of the conducting band. This type of mech- 
anism certainly does play a role in (III,Mn)V ferromag- 
netism, likely dominating in the case of strongly metal- 
lic (Ga,Mn)As, (In,Mn)As, and Mn-dopcd antimonides. 
There is no sharp distinction between impurity band 
double-exchange and kinetic-exchange interactions; the 
former is simply a strong coupling, narrow band limit of 
the latter. 

The starting point for developing a useful predictive 
model of (III,Mn)V ferromagnetism is achieving a full 
understanding of the electronic state of a single Mn im- 
purity in the host lattice. We need to fully understand 
the character of the isolated local moments before we can 
critically discuss how they are coupled. The character of 
the local moment need not be the same in all (III,Mn)V 
materials. The remaining subsections will focus on prop- 
erties of a Mn impurity in GaAs and on the nature of fer- 
romagnetic coupling in (Ga,Mn)As and related arsenide 
and antimonidc DMSs. At the end of this section, we 
comment on how things might change in wider-gap hosts 
like GaP and GaN. 



B. Substitutional IVIn impurity in GaAs 

Among all (III,V) hosts, Mn impurity has been 
studied most extensively in GaAs. The elements in 
the (Ga,Mn)As compound have nominal atomic struc- 
tures [Ar]3fii°4s2pi for Ga, [Ar]3d^4s^ for Mn, and 
[Ar]3(i^"4s^p'^ for As. This circumstance correctly sug- 
gests that the most stable and, therefore, most common 
position of Mn in the GaAs host lattice is on the Ga 
site where its two 4s-electrons can participate in crys- 
tal bonding in much the same way as the two Ga 4s- 
electrons. The substitutional Mnca, and the less com- 
mon interstitial Mnj, positions are illustrated in Fig. O 
Because of the missing valence 4p-electron, the Mnca im- 
purity acts as an acceptor. In the electrically neutral 
state, labeled as A^{d^ + hole), MuQa has the charac- 
ter of a local moment with zero angular momentum and 
spin S = 5/2 (Lande g- factor g = 2) and a moderately 
bound hole. The local moment is formed by three occu- 
pied sp — d bonding states with dominant t2g (Sd^y, -id^z, 
3dyz) character and by two occupied eg (3^3.2 _j,2, 3^22) 
orbitals that arc split from the t2g states by the tetrahe- 
dral crystal field and do not strongly hybridize with the 




FIG. 3 Top panel: Substitutional Mnoa and interstitial Mni 
in GaAs. Bottom panel: two eg 3ci-orbitals and three t2g 
3d-orbitals of Mn. 



sp-orbitals. All occupied d-orbitals have the same spin 
orientation and together comprise the S = 5/2 local mo- 
ment. The weakly bound hole occupies one of the three 
antibonding sp — d levels with dominant As Ap character. 
The charge — e ionized MuQa acceptor center, labeled as 
has just the 5" = 5/2 local spin character. 

Electron paramagnetic resonance (EPR) and ferro- 
magnetic resonance (FMR) experiments confirm the 
presence of the A~{d^) center through the entire 
range of M n concentrations in both bulk and epilayer 
(Ga,Mn)As llAlmeleh and Goldsteiiil Il962t ISasaki all 
l2002bl: ISzczvtko et o,l\ . Il999b^ . The S = 5/2 lo- 
cal moment on Mn was detected through a resonance 
line centered at g = 2 and, in low Mn-density sam- 
ples, through a sextet splitting of the line due to 
the hyperfine interaction with the / = 5/2 ^^Mn nu- 
clear spin. The neutral MuQa centers arc more elu- 
sive because of nearly full compensation by uninten- 
tional donor impurities at low Mn concentrations and 
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FIG. 4 STM imaging of a Mnca impurity. Top panel: energy 
band diagram for the negative (left) and positive (right) bias. 
Bottom panel: Mnca impurity in the ionized A~(c?^) (left) 
and the neutra l yl°(d^ + hole) (right) state. Adapted from 
(|Yakunin et aZ.L l2004bh . 



because of the metal-insulator transition at high Mn 
concentrations. Nevertheless, a multitude of exp e rimen - 



^^^.^^ ^^^^^V-^^ Q 

infrared fIR) spectroscopy 


(IChaDman and Hutchinsonl 


119671: iLinnarsson et al\ '1997 


1), and magnetization mea- 


surements (,Frev et at 


1988|) 


, have detected the A°{d^ + 



hole) center in (Ga,Mn)As. Strik ingly direct ev i dence 
was give n by an STM e xperiment l)Kitchen et all l2005t 
lYakunin et all. |200433), shown in Fig. 01 in which the 
state of a single impurity atom was switched between the 
ionized A~{d^) and the neutral A'^{(f + hole) by apply- 
ing a bias voltage that corresponded to a binding energy 
Eb^QA eV. 



The binding energy Eb 



112.4 meV inferred 



from IR spectros c opy JChapman and Hutchinsonl . Il967t 
ILinnarsson et al] . Il997t) is consistent with the above 
STM measurement and with inferences based on pho- 
toluminescence exper i ments (|Lee and AndersonL Il964t 
ISchairer and Schmidtl . ^^). These observations iden- 
tify Mn as a moderately shallow acceptor in GaAs whose 
band gap is Eg = 1.52 eV. The binding energy, which 
governs the electrical behavior of the Mn impurity, has 
contributions from both Coulomb attraction between the 
hole and the A~{(f) core (and a central cell correction) 
and spin-dependent p—d hybridization. The latter effect 
is responsible for the exchange interaction on which this 
review centers. We now discuss this parameter in more 
detail. 



1. p-d exchange coupling 

The top of the GaAs valence band is dominated by 
4p-levels which arc more heavily weighted on As than 
on Ga sites. Direct exchange between the holes near 
the top of the band and the localized Mn d-electrons 
is weak since Mnca and As belong to different sublat- 
tices. This fact allows p—d hybridization to domi- 
nate, exp laining the antiferromagne tic sign of this in- 
teraction ijEhattacharjee et all Il983|) seen in experiment 
(|Okabavashi et all \199^ . 

There is a simple physical picture of the p — d exchange 
interaction which applies when interactions are treated in 
a mean-field way, and therefore also applies as an inter- 
pretation of LSDA calculations. Given that the filled, 
say spin-down, Mn d-shell level is deep in the valence 
band and that the empty spin-up d-level is above the 
Fermi level and high in the conduction band, hybridiza- 
tion (level repulsion of like-spin states) pushes the energy 
of spin-down valence band states up relative to the energy 
of spin-up valence band states. The resulting antiferro- 
magnetic coupling between valence band states and local 
Mn spins is illustrated schematically in Fig. jS] The same 
basic picture applies for itinerant valence band states in 
a heavily-doped metallic DMS and for the acceptor state 
of an isolated MuQa impurity. Note that the cartoon 
band structure in Fig.jSjis plotted in the electron picture 
while the DMS literature usually refers to the antiferro- 
magnetic p — d coupling between holes and local Mn mo- 
ments. We comment in detail on the equivalent notions 
of p — d exchange in the physically direct electron-picture 
and the computationally more convenient hole-picture for 
these p-type DMSs in Section IV.BI 



Single Mn(d + hole) 

d\ 




hybridization 

Many-Mn system 
— d\ 



■d\ 




FIG. 5 Electron-picture cartoon: splitting of the isolated Mn 
acceptor level (top panel) and of the top of the valence-band 
in the many-Mn system (bottom panel) due to p-d hybridiza- 
tion. 

We have already mentioned in Section III.AI the con- 
ceptional inadequacy of effective single-particle theories, 
including LSDA, in dealing with local moment levels co- 
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incident with itinerant electron bands. Anderson sug- 
gested a many-body model Hamiltonian that circumvents 
this problem by including explicitly the Coulomb cor- 
relation inte gral of the localiz ed electron states i n the 
HamiltonianlA'n derson. .l961: .FTeurov and KikoinLll976t 
iHaldane and AndersonL Il97(1^ . The problem is that the 
change in the effective potential when the number of 
occupied localized orbitals changes by one, the Hub- 
bard constant U, can be comparable to or larger than 
other band parameters, invalidating any mean- field- like 
approach. The consequences of this fact can be captured 
at a qu alitative level in mod els that include the Hub- 
bard U l)Krstaiic et all l2004() . In these phenomenolog- 
ical models, the localized orbital part of the Hamilto- 
nian generally has an additional parameter, the Hund's 
rule constant Jh- This parameter captures the local di- 
rect exchange physics which favors spin-polarized open 
shell atomic states. For the case of the Mn((i'^) config- 
uration, Jh forces all five singly-occupied d-orbitals to 
align their spins in the ground state. Recently consid- 
erable effort has been devoted to developing approaches 
that combine the local correlation effects induced by the 
U and Jh terms in phen o menol o gical models, wi t h SDF 
theorv llAiiisimoy et al\. Il99lt iFihppetti et ail I200E : 
Park et all. l200nHPerdew a,nd Zungeri ll 98lHPetit"^ . 
2006"; 'Sandratskii et a/i l2004t ISchiilthess et all 12005 : 
Wierzbowska et al., . .20041) . We comment on these ab ini- 
tio techniques in Section IlII.AI 

When hybr idizat ion jHarrisonl Il980t 

ISlater and Kosteil Il954 between the local moment 
and band electron states is weak it can be treated 
perturbatively. The Schrieffer- Wolff transformation of 
the Anderson Hamiltonian, 

Ha = CfeJifcs + ^ edUds + U Ud^Udi 

k.s s 

+ Y.^Vkdcl,Cds+C.C.) , (1) 

ks 

removes the hybridization (last term in Eq. |^ and leads 
to a model in which the local moment spin interacts 
with the valence band via a spin-spin interaction only, 
ifc'fcSrf • Sfc'fci with the number of electrons in each 
band fixed iSchrieffer a,nd WolfUFigfifil) . Here we assume 
for simplicity a single localized orbital and a single itin- 
erant band, we use k to represent band states and d 
to represent the localized impurity state, s labels spin, 
e„ a single particle energy, and n„ = cJ^Cq, and are 
the standard second quantization operators. This proce- 
dure is normally useful only if the hybridization is rela- 
tively weak, in which case it is not usually a problem to 
forget that the canonical transformation should also be 
applied to operators representing observables. Strictly 
speaking, the Schrieffer- Wolf tran sformation also leads 
to a spin-independent interaction ijSchrieffer and WolfJ . 
Il966j) which is normally neglected in comparison with 
the stronger spin-independent long-range part of the 
Coulomb potential. 



Since the valence band states of interest in GaAs, even 
in heavily doped samples, are near the Brillouin-zone 
center F, the single phenomenological constant extracted 
from experiment for this interaction should be thought of 
as its value when both initial and final states are at the 
F-point. The quantity '^f.i^jk'k^d^k'k is then approxi- 
mated by JoSrf • s/j=o, where 

Vpd represents the As p-orbital - Mn d-orbital hybridiza- 
tion potential (neglecting again the multiple orbital na- 
ture of the p and d- levels for simplicity), < is the 
single-particle atomic-level energy of the occupied Mn d- 
state measured from the top of the valence band, and 
ed + U > is the e nergy cost of adding a seco nd elec- 
tron to this orbital ijSchrieffer and Wolfl. 1196^ . When 
the hole density is large or the holes are more strongly 
localized near Mn acceptors, the crystal- momentum de- 
pendence of this interaction parameter ca nnot be entirely 
neglected l|Timm and MacDonaldL 120051) . 

Since Jq originates from hopping between Mn and the 
neighboring As atoms, the p ~ d exchange potential, 
J(Ri — r), produced by Mn impurity at site Ri has a 
range of order one lattice constant, and 

J ^ Lc drJ{R,-r)u*„{r)un{r) 
dr<(r)^„(r) ' 

In Eq. Q we assumed that the perfect crystal Bloch func- 
tion, ?/'ri,k(r) = exp(ikr)M„,k(r), is composed of a slowly 
varying envelope function and a periodic function u,i_k(r) 
with the normalizations 1/1/ J rfri/;* j,(r)7/'„'.k'(r) = 
Sn,n'Sk,k' and 1/ilu.c. /„ ^ dr<(r)u„'(r) = Sn,n'- Here V 
is the crystal volume and ^u.c. is the unit cell volume. In 
GaAs, r^M.c. = af^/4=0.045 nm'^ and the lattice constant 
aic = 0.565 nm. These wavefunctions can be obtained 
from k • p theory which treats the band Hamiltonian of 
the system perturbatively expanding arou nd the F-point 
ijAbolfath et a?J . l2001atlDietl adl2001hj) . Jo inEq. @ 
corresponds to the average value of J(Ri — r) seen by the 
k = Bloch state over the i-th unit cell. 

The slowly varying envelope function sees an effective 
zero-range p — d exchange potential, since 

(-0„,k| J(R.j - 
= y"drexp [i(k- k')r] J(Ri - r)u*_k(r)w„',k'(i") 

« -r^ ^ exp [lik - k')R,] / drJ(R, - r)<(r)u„(r) 

^ JV„... 

= — ^ exp [i(k - k')Rj] Jon„.c.(5jj 

~ ^ j dv exp [i(k - k')r] Jpd^(R. - r) , (4) 

where iV„ c. is the number of unit cells in the crystal 
volume V and Jpd = Jq^u.c- Note that the assumption 
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of zero-range in the k • p treatment corresponds to the 
neglect of momentum dependence of this interaction for 
wavevectors close to the F-point in a microscopic picture. 
Positive sign of the Jpd (or Jq) constant corresponds to 
an antiferromagnetic p — d coupling. We also note that in 
the II- VI DMS literature a symbol (3 is often used instead 
of Jpd and NqP instead of Jo, where A^o — l/^n.r, is the 
number of unit cells (or cation s ites) per volume ijPietj 
ll994HFurdvna and Kossul Il988() . 

The value of the Jpd constant is often c onsidered to 
be independ e nt of the host se i iiicondu ctor ijPietl et al\ . 
I2001bl l2000t lJungwirth eloll. l2002bl) Indeed the in- 
crease of ~ ar" l|Harrisoni ll98(Tl) in Eq. © with 
decreasing lattice constant is partly compensated by 
smaller r2„.c. (~ afc)j and the increase of |l/ed| in a wider 
gap host is partly compensated by the decrease of the 
term \ / {ed + U) . Although it may have similar values in 
many materials, Jpd will tend to be larger in larger gap, 
smaller lattice constant hosts. 

In a virtual crystal mean-field approximation, the 
p — d exchange potential due the Mn impurities in a 
Gai_:,Mn^As DMS, xNoVLu.c. T.r^,,, >/(R«.c. - r)(S) • s, 
has the periodicity of the host crystal. (Here (5) is 
the mean-field Mn spin.) The valence band states 
in this approximation experience an effective single- 
particle kinetic-exchange field, h.MF = NMnJpdi^) , 
where Nmu = ^Nq is the Mnca density. 

Finally, we discuss the relationship between the ex- 
change constant Jpd and the exchange constant, £, used 
to provide a k • p interpretation of spectroscopic stud- 
ies of the neutral A'^{d^ + hole) center. We emphasize 
that the use of a k • p approach assumes that the bound 
hole is spread over at least several lattice constants in 
each direction. The fact that it is possible to achieve a 
reasonably consistent interpretation of the very detailed 
spectroscopic data in this way is in itself strong sup- 
port for the validity of this assumption. The coupling 
between the weakly bound hole moment J and the lo- 
cal spin S of the Mnca core is expressed in the form 
eS • J, where J = j + L, j is the (atomic scale) total 
angular momentum operator of the band-hole at the F- 
point (j = 3/2 or 1/2 for the As 4p-orbitals forming the 
band states near k = 0), and L is the additional (hole 
binding radius scale) angular momentum acquired by the 
hole upon bind ing to the Mnpia i mpuri ty. The IR spec- 
trosco pv data jLinnarsson et all Il997t) have b een ana- 
lyzed l)Bhattachariee and a la Guillaum^ |2000() within a 
spherical approximation, i.e., considering only the L = 0, 
s-like bound state. (Note that a sizable anisotropic d-like 
component in the bound-hole ground st ate has been iden- 
tified in the analysis of the S TM data ijTang and Flattd . 
I2004t lYakunin et all hoOAlt} .) Further simplification is 
achieved by the neglect of the admixture of the two 
j = 1/2 {jz = il/2) F-point states which is justified 
by the large spin-orbit splitting, Aso = 341 meV, of 
these states from the two heavy- hole states (j = 3/2, 
jz = ±3/2) and the two light-hole states (j — 3/2, 
J. = ±1/2). 



Writing the ground-state wavefunction in the form 
(r) = Fj^ (r)ui- (r) with a spherically symmetric enve- 
lope function ^^.(r) and for = ±3/2, ±1/2 {j = 3/2), 
the expectation value of the exchange potential reads 

(V'.J J(R7 - r)S • s\i',J = \f{Ri)\^-Jpd{jz\S ■ s\jz) 

= L7(R7F^(j.|s-j|j.)(^) 



where |/(R/)p =< F*^{r)Fj^{r) >u.c. is the mean 
value of the square of the slowly varying envelope func- 
tion within the unit cell containing the MuQa impurity. 
Eq. jSJ implies that 



£=^|/(R,)P, (6) 

i.e., the ratio between e and Jpd is determined by the 
strength of the binding of the hole to the A^{d^) Mnca 
core and is larger for more localized holes. 

A combination of IR data and theoretical cal- 
culations has been used to analyze this in more 
detail llBhattachariee and a la Giiillauin3 . l200flt 
iLinnarsson et alV Il997t) First of all, the value of 
the g-factor g = 2.77 of the neutral A°{d^ -\- hole) MuQa 
complex obtained from IR spectroscopy measurements 
is in agreement with the theoretical value expected for 
the total angular momentum state F = S — J = \ 
of the complex, confirming the antiferromagnetic 
character of the p — d coupling between the hole and 
the local Mn spin. The contribution of the p ~ d 
potential to the binding energy is then given by 
eS • J = e[F{F + 1) - S{S + 1) - J(J + l)]/2 which 
for the F — 1 ground state gives — 21e/4. The IR 
spectroscopy measurement of the splitting 2e between 
the F = 1 and F ~ 2 states gives e w 5 mcV, i.e., 
the contribution to the binding energy from the p—d 
interaction is approximately 26.25 meV. The remaining 
binding energy, 112.4-26.25=86.15 meV, is due to the 
central field potential of the impurity. Bhatta c harjee 
and Benoit l|Bhattachariee and a la GuiUaiimel. 120001) 
used the hydrogenic-impurity model with a screened 
Coulomb potential and a central cell correction, whose 
strength was tuned to reproduce the value 86.15 meV, to 
obtain a theoretical estimate for |/(R/)P w 0.35 nm""^. 
From Eq. © they then obtained Jpd ~ 40 meV nm"^. 
Given the level of approximation used in the theoretical 
description of the A°{d^ + hole) state, this value is in a 
reasonably good agreement with the exchange constant 
value Jpd = 54 ± 9 meV nm^ (Na = 1.2 ± 0.2 eV) 
inferred from photoemission data ijOkabavashi et ali . 
Il998l). We no te here that photoemissio n spect r oscopy 
jHwang et aO, '20 051 lOkabavashi et ali Il999l Il998l 
I2OO I. 2002; .Radcr et all |2004|) has represented one 
of the key experimental tools to study the properties 
of Mn impurity in DMSs, in particular, the position 
of the Mn d-level and the strength and sign of the 
p — d coupling. An indirect measurement of the Jpd 
constant is performed by fitting the photoemission data 
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to a theoretical spectrum of an isolated MnAs4 cluster 
(|Okabavashi et all Il998j) . This procedure is justified 
by the short-range character of the p — d exchange 
interaction. 



C. Other common impurities in (Ga,Mn)As 

Most of the single Mn-impurity spectroscopic studies 
mentioned above were performed in samples with doping 
levels X < 0.1% for which the Gai_2.Mn3,As random al- 
loy can be grown under equilibrium conditions. In these 
materials Mn can be expected to occupy almost exclu- 
sively the low energy Ga-substitutional position. Ferro- 
magnetism, however, is observed only for x > 1% which 
is well above the equilibrium Mn solubility limit in GaAs 
and, therefore, requires a non-equilibrium growth tech- 
nique (in practice low-temperature MBE) to avoid Mn 
precipitation. The price paid for this is the occurrence 
of a large number of metastable impurity states. The 
most important additional defects are interstitial Mn ions 
and As atoms on cation sites (antisite defects). Both act 
as donors and can have a severe impact on the electric 
and magnetic properties of the DMS epilayers. More un- 
intended defects form at higher Mn doping because of 
the tendency of the material, even under non-equilibrium 
growth conditions, toward self-compensation. 



1. Interstitial Mn 

Direct experimental evidence for Mn impurities occu- 
pying interstitial (Mni) rather than substitutional po- 
sitions was uncovered by combined channeling Ruther- 
ford backscattering and part i cle in duced x-ray emis- 
sion measurements dYu" et al 1 l2002|) . This technique 
can distinguish between Mnj and Mnoa by counting 
the relative number of exposed Mn atoms and the 
ones shadowed by lattice site host atoms at differ- 
ent channeling angles. In highly doped as-grown sam- 
ples, the experiment identified nearly 20% of Mn as 
residing on interstitial positions. The metastable na- 
ture of these impurities is manifested by the substan- 
tial decrease in their density upon post-growth anneal- 
ing at te mperatures very clo s e to the growth tem- 
peratures fChiba et all l2003at lEdmonds et al\. l2002at 
iKu et a^J . i2003i: iStone et oi.L l200,'4 lYu et fl^.Ll2002^ ■ De- 
tailed resistance-monitored annealing studies combined 
with Auger surface analysis established the out-diffusion 
of Mni impurities towards the free DMS e pilayer surface 
during annealing l)Edmonds et g^J . l2004a|) . The charac- 
teristic energy barrier of this diffusion process is esti- 
mated to be 1.4 eV. (Note that a factor of 2 was omitted 
in the original estimate of this energy in (jEdmonds et all 
l2004aD .1 

Isolated Mnj spectroscopy data are not available in 
(Ga,Mn)As, underlying the importance of theoretical 
work on the electric and magnetic nature of this impurity. 



Density functional calculations jErnst et all l2005(l sug- 
gest, e.g., that minority-spin Mni d-states form a weakly 
dispersive band at ^ 0.5 eV below Fermi energy, a fea- 
ture which is absent in the theor etical Mnn^ spectra. Ab 
initio total energ y calculations jMaca and Masefl l2002t 
iMasek and Macat l2003(l showed that Mn can occupy two 
metastable interstitial positions, both with a compara- 
ble energy, one surrounded by four Ga atoms (see Fig.|31) 
and the other surrounded by four As atoms. The two Mnj 
states have similar local magnetic moments and electro- 
negativity. 

The calculations have confirmed that Mnj acts as a 
double-donor, as expected for a divalent metal atom oc- 
cupying an interstitial position. Each interstitial Mn 
therefore compensates two substitutional Mn acceptors. 
It seems likely that because of the strong Coulombic at- 
traction between positively charged Mni and negatively 
charged Mnca defects, the mobile interstitials pair-up 
with substituti onal Mn during the grow th, as illustrated 
in Fig. 151 ( Blin owski and Kacmanll2003(l . The total spin 
of a MnQa-Mni pair, inferred from ab initio calculations 
is much smaller than the local spin S = 5/2 of the 
isolated MuQa acceptor. This property, interpreted as 
a consequence of short-range antiferromagnetic interac- 
tions between two local moment defects that have com- 
parab le local moments, has been confirmed experimen- 
tally ijEdmonds et a/J . r2005a|l . According to theory, the 
strength of this magnetic coupling c ontribution to the 
Mnoa-Mni binding energy is 26 meV ijMasek and Macal 

Ab-initio calculations of the spin-splitting of the va- 
lence band indicate that the Jpd coupling constants 
of interstitial a nd substitutional Mn are comparable 
(jMasek and Mac a. 2003). This would suggest a neg- 
ligible net p — d coupling between the antiferromag- 
nctically coupled Mnca-Mni pair and the valence band 
holes. We do note, however, that these LSDA calcu- 
lations give Jpd ~ 140 meV nm"^ which is more than 
twice as large as the experimental value, so this conclu- 
sion must be regarded somewhat cautiously. It is gen- 
erally accepted that this discrepancy reflects the general 
tendency of the DFT-LSDA and similar theories to sys- 
tematically underestimate the splitting between the oc- 
cupied Mn d states in the valence band continuum and 
the empty Mn c?-states. Eq. illustrates how this defi- 
ciency of SDF band theories translates into an overesti- 
mated strength of the p — d exchange interaction. Apart 
from this quantitative inaccuracy, the relative strength of 
the p — d interaction of Mnj compared to MuQa is still a 
somewhat controversial issue in the theoretical lit erature 
l)Bhnowski and Kacma^ l2003t iMasek and Macal l2003(l . 
which has not yet been settled experimentally. This prop- 
erty, and many others related to magnetism are sensitive 
to details of the SDF implementation, to the way in which 
disorder is accounted for, and even to technical details 
associated with the way in which the Bloch-Schrodinger 
equation is solved numerically. For example, the positive 
charge of the Mnca-Mnj (single-acceptor double donor) 
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pair, which will tend to reduce its exchange coupling with 
valence band holes due to Coulomb repulsion, is not in- 
cluded in all approaches. 

In addition to the direct hole and local moment com- 
pensation effects of Mni defects on ferromagnetism in 
(Ga,Mn)As, the structural changes they induce in the 
crystal arc indirectly related to important magnetic prop- 
erties, particularly to the various magnetic and transport 
anisotropics. Ab initio theory predicts that the separa- 
tion of the four nearest As neighbors surrounding the 
Mni in a relaxed lattice is increased by 1.5% compared 
to the clean GaAs lattice ijMasek et a^J . [2003(l . Because 
of the Coulomb repulsion between Ga cations and Mni 
defects, an even larger lattice expansion (~ 2.5%) is 
found for the Mnj in the four Ga tetrahedral position. 
When grown on a GaAs substrate, this effect of inter- 
stitial Mn leads to a lattice-matching compressive strain 
in the (Ga,Mn)As thin layers that induces a large uni- 
axial magnetic anisotropy, as we discuss in detail in Sec- 
tion ESI 



2. As antisites 

Low-temperature growth of GaAs is known to lead to 
the incorporation of high levels of As antisite defects. 
This property is a combined consequence of the non- 
equilibrium growth conditions and the As overpressure 
often used in the MBE process to assure the 2D growth 
mode. These double-donor defects are likely to be also 
present in the (Ga,Mn)As epilayers and may contribute 
to hole compensation. 

Unlike Mnj impur i ties. A s antisites are stable up to 
~450 °C l)Bliss et all\l99'^ . This is well above the tem- 
perature at which Mn precipitation starts to dominate 
the properties of (Ga,Mn)As and, therefore, the As an- 
tisites cannot be removed from the epilayer by a post- 
growth annealing treatment. Experimental studies sug- 
gest that the degradation of (Ga,Mn)As magnetic prop- 
erties due to hole compensation by As antisites can be re- 
duced by using As2 dimers instead of As4 tetramers and 
by maintainting a stric tly stoichiometric growth mode 
(|Gampion et ad l2003bj) . 



D. Qualitative picture of ferromagnetism in (Ga,l\/ln)As 
and other (lll,i\/ln)V materials 

The following elements of the qualitative picture of fer- 
romagnetism in (Ga,Mn)As emerge from the experimen- 
tal data and theoretical interpretations discussed in this 
section (see also (|Dietll . l2002|) '). The low-energy degrees 
of freedom in (Ga,Mn)As materials are the orientations 
of Mn local moments and the occupation numbers of ac- 
ceptor levels near the top of the valence band. The num- 
ber of local moments participating in the ordered state 
and the number of holes may differ from the number of 
MuGa impurities in the III-V host due to the presence of 



charge and moment compensating defects. Hybridization 
between Mn d-orbitals and valence band orbitals, mainly 
on neighboring As sites, leads to an antiferromagnetic 
interaction between the spins that they carry. 

At low concentrations of substitutional Mn, the 
average distance between Mn impurities (or between 
holes bound to Mn ions), rc = {S/AnNMnY^^ , is much 
larger than the size of the bound hole characterized 
approximately by the impurity effective Bohr radius, 
a* ~ efi/m*e'^. Here Nmu = 4x/af^ is the number 
of Mn impurities per unit volume, e and aic are the 
semiconductor dielectric function and lattice constant, 
respectively, and m* is the effective mass near the top 
of the valence band. For this very dilute insulating 
limit, a theoretical concept was introduced in the late 
1970's in which a ferromagnetic exchange interaction 
between Mn local moment s is mediated by thermally 
activ ated band carriers (|Pashitskii and Rvabchenka 
Il979() . Experimentally, ferromagnetism in (Ga,Mn)As 
is observed whe n Mn doping reaches a pproximately 1% 
llGarnpion et alV l2003at lOhnd Il999t IPotashnik et all 
I2002D and the system is near the Mot t insul ator-to-metal 
transition, i.e.. Tr ~ a * llMardeil Il999ft . At these 
larger Mn concentrations, the localization length of the 
impurity band states is extended to a degree that allows 
them to mediate ferromagnetic exchange interaction 
between Mn moments, even though the moments are 
dilute. Several approaches have been used to address 
ferromagnetism in DMSs near the metal-in sulator 
transition (for a review see ijPietll. l2002t iTimnJ. l2003i|) . 
including finite-size exact diagonalization studies of 
hole-hole an d hole-impurity Coulomb interaction effects 
(jTimm et al . 2002: Yang and MacDonald ,, 2003) , and 



holes hopping within the impurity band (lAlvarez et al. 


20021 iBerciu and Bhatt, ,2001; Bhatt et all 


2002 


Chudnovskiv and Pfannkuchel l2002l iDurst et al\. 


2002 


Inoue et all l2000l ll 


Caminski and Das Sarmal 


2002 


Litvinov and Dueaevl 


200 ll). Because of the hoppine 



nature of conduction and the mixed spd character of 
impurity band states, the regime is sometimes regarded 
as an example of double-exchange ferromagnetism. 

At even higher Mn concentrations, the impurity band 
gradu ally merges with the valence band ijKrstaiic et all 
I2004D and the impurity states delocalize. In these 
metallic (Ga,Mn)As ferromagnets, which we focus on 
in the following sections, the coupling between Mn lo- 
cal moments is mediated by the p — d kinet ic- ex c hange 
mechanis m llDietl et all Il997l l2000t iJain et all 1200 ll 
Uunawirt h et aZ.I. ll999HMatsukura et all\l99^ A qual- 
itatively similar picture applies for (In,Mn)As and Mn- 
doped antimonides. In the metallic limit the influence 
of Coulomb and exchange disorder on perfect crystal 
valence-band states can be treated as perturbatively. 

The crossover from impurity-band mediated to Bloch 
valence-band mediated interactions between Mn mo- 
ments is a gradual one. In the middle of the crossover 
regime, it is not obvious which picture to use for a 
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qualitative analysis and quantitative calculations are 
not possible within either picture. Strongly localized 
impurity-band states away from Fermi energy may play a 
role in spectroscopic properties l)Okabavashi et all . l200l|) . 
even when they play a weaker role in magnetic and 
transport properties. The crossover is controlled not 
only by the Mn density but (because of the impor- 
tance of Coulomb interaction screening) also by the car- 
rier density. There is a stark distinction between the 
compensation dependence predicted by impurity-band 
and Bloch valence-band pictures. When the impurity- 
band picture applies, ferroma gnetism does not occur 
in the absence of compen s ationllDas Sarma et a H. l2003t 



iKaminski and Das Sarmai l2fl0.'^ ScarpuUa et oil 1200,^ . 

because the impurity band is filled. Given this, we 
can conclude from experiment that the impurity band 
picture does not apply to optimally annealed (weakly- 
compensated) samples which exhibit robust ferromag- 
netism. 

The phenomenology of Mn-doped phosphide or nitride 
DMSs is more complex, with many aspects that prob- 
ably ca nnot be captured by the Zener kinetic -exchange 
model l|Dietl et ad . l2002HKrstaiic et aZ.ll200^ . Experi- 
mentally, the nature of the Mn impurity is very sensitive 
to the presence of other impurities or defects i n the lattice 
(lArkun et all l2004t lEdmonds et all l2004d iGraf et all 



120021 12003al: iHwang et ad l20n,'THKorotkov et all\20dlt) . 

This substantially complicates the development of a con- 
sistent picture of ferromagnetism in these materials. 
Larger p—d coupling in the wide ga p DMSs (reported e. g. 
in the photoemission experiment l|Hwang et all l2005j) '). 
and stronger bonding of the hole to the Mn ion, might 
shift the metal-insulator transition to higher Mn densi- 
ties. At typical dopings of several per cent of Mn, the 
impvirity band is st ill detached from the valence band 
iKronik et a/.j2002(l and ferromagnetic Mn-Mn coupling 
is mediated by holes hopping within the impurity band. 
Recent experiments indicate that this s cenario may ap- 
ply t o (Ga,Mn)P with 6% Mn dopin g (IScarpulla et al 
2005 ) . LSDA calculations suggest (ISandratskii et al. 



2004[lSanval all l2003t IWierzbowska e7a?J . l2004|) that 



the p—d hybridization can be so strong that the ad- 
mixture of Mn 3d spectral weight at the Fermi energy 
reaches a level at which the system effectively turns into 
a d-band metal. To illustrate this trend wc show in 
Fig. El the LSDA and LDA-I-U calculations of the spin- 
split total density of states (DOS), and in Fig. [7| the re- 
sults for the M n d-states projected DO S in (Ga,Mn)As 
and (Ga,Mn)N (|Sandratskii et a?J . l200^ . Indeed, in the 
wider gap (Ga,Mn)N the spectral weight of Mn d-orbitals 
at the Fermi energy is large and is not significantly sup- 
pressed even if strong on-site correlations are accounted 
for by introducing the phenomenological Hubbard pa- 
rameter in the LDA-I-U method (see next section for de- 
tails on combined Hubbard model and SDF techniques). 

Another possible scenario for these more ionic III-P 
and III-N semiconductors, supported by EPR and opti- 
cal absorption measurements and ab initio calculations 
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FIG. 6 The spin-resolved total DOS of (GaMn)As and 
(GaMn)N for a Mn concentration of 6.25%. Calculations 
are performed within LSDA and LDA+U approaches. For 
(GaMn)N, both zinc-blende and wurzite structures are pre- 
sented. The spin-up/spin-down DOS is shown above/below 
the abscissa axis. The total DO S is given per chemical unit 
cell of the semiconductor. From iSandratskii et all 1200411 . 



(iGraf et all l2002l iKreissl et all Il996t iLuo and Martini 
l2005HSchulthess et a/J . l200,5l) . involves a transition of the 
substitutional Mn from a divalent (d^) impurity to a 
trivalent (d^) impurity. This strongly correlated d^ cen- 
ter, with four occupied d-orbitals and a non-degenerate 
empty d-level shifted deep into the host band gap, may 
form as a result of a spontaneous (Jahn- Teller) lowering 
of the cubic symmetry near the Mn site. If the energy 
difference between divalent and trivalent Mn impurity 
states is small, the DMS will have a mixed Mn-valence 
which evokes the conventional double-exchange mecha- 
nism. Systems with dominant d^ character of Mn im- 
purities, reminiscent of a charge transfer insulator, will 
inevitably require additional charge co-doping to provide 
for ferromagnetic coupl ing between dilute Mn moments 
ijSchulthess adl2005|) . 



III. THEORETICAL APPROACHES 

Our focus is on the theory of (III,Mn)V ferromagnets 
and we therefore present in this section an overview of 
the different approaches that can be used to interpret the 
existing experimental literature and to predict the prop- 
erties of materials that might be realized in the future. 
Since the electronic and magnetic properties of ferromag- 
netic semiconductors are extremely sensitive to defects 
that are difficult to control in real materials and may 
not be completely characterized, the ability to make rea- 
sonably reliable theoretical predictions that are informed 
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FIG. 7 Same as in Fig. M but for partial Mn 3d DOS. In 
LDA+U the p — d hybridization in (Ga,Mn)As is weak for 
both majority band and minority band states. 



by as many relevant considerations as possible can be 
extremely valuable to the effort to discover useful new 
materials. For example we would like to be able to make 
confident predictions of the ferromagnetic transition tem- 
perature of a (III,Mn)V material as a function of the 
density of substitutional Mn, interstitial Mn, co-dopants, 
antisite defects, and any other defects whose importance 
might be appreciated in the future. This ability is devel- 
oping, we believe, although there is no simple silver bullet 
that solves all difficulties for all host materials, and there 
may still be some considerations that are important for 
less studied host materials and are not yet part of the 
discussion. In this section we address in a general way 
the strengths and weaknesses of some of the different the- 
oretical approaches. The following sections of the review 
compare predictions made with different types of theo- 
retical approaches with experimental data on a variety of 
important electronic and magnetic properties. 



A. First-principles calculations 



In SDF th e ory ijHohenberg and KohnL 119641: 
iKohn and Shanj . Il965^ all many-body effects are 
buried in a complex exchange-correlation energy 
functional. Once an approximation is made for this 
functional, predictions for electronic and magnetic 
properties depend only on the particular arrangement 
of atomic nuclei under consideration. In principle 
nuclear positions can be relaxed to make sure that 
the spatial distribution of nuclei is mctastable and 
therefore realizable. The exchange-correlation energy 
functional leads to a self-consistcntly determined spin- 



dependent exchange-correlation potential that appears 
in an effective independent-particle Hamiltonian. The 
main technical challenge in DPT theory applications is 
the development of numerically efficient methods that 
provide accurate solutions of s ingle-body Schrodinger 
equations (see revi e w art icles ijJones and Gunnarssonl 
Il989t ISanvito et all [200^ 1. DFT theory is established 
as a flexible and valuable tool for studying the micro- 
scopic origins of magnetism and for predicting electronic, 
magnetic, and ground-st ate structural properties in a 
wide- variety of mate ri als llJones and Gunnarssonl 119891: 
iMoruzzi and Marcu j. Il993|) . It has the advantage 
that it is a first principles approach without any phe- 
nomenological parameters. DFT theory falls short of 
being a complete and general solution to the many- 
electron problem only because the exact form of the 
exchange-correlation energy functional is unknown. A 
simple and widely successful approximation is the LSDA 
(|von B arth and Hcdin, 1972). 

The problem of solving LSDA equations with adequate 
accuracy remains a challenge even in perfectly ordered 
crystals. In DMSs the degrees of freedom that are im- 
portant for ferromagnetism, the orientations of the Mn 
local moments, typically reside on approximately 1/40 of 
the atomic sites which further complicates numerical im- 
plementation of the LSDA technique. Other length scales 
that are characteristic for the physics of interest in these 
materials, like the Fermi wavelengths of the valence band 
carriers, are also longer than the atomic length scale on 
which DFT theory interrogates matter. This property 
limits the number of independent magnetic degrees of 
freedom that can be included in a DFT simulation of 
DMS materials. The problem is exacerbated by the alloy 
disorder in (Ga,Mn)As. Even if all Mn atoms substitute 
for randomly chosen lattice sites, it is necessary to find 
some way to average over microrealizations of the alloy. 

Disor der-ave r aging coherent-potential ap proximation 
(CPA) (ISoveni Il967t IVelickv et all Il968ii and super- 
cell approaches, have been used successfully in combi- 
nation with DFT calculations to address those phys- 
ical parameters of (III,Mn)V DMSs that are de- 
rived from total-energy calculations, such as the lat- 
tice constants, formation and binding energies of var- 



(see e.a. llEdmonds et a 


IL l2004al lErwin and Petukhov 


20021: iLuo and Martin'. ' 


20051 iMaca and Masekl 12002 


Mahadevan et all 12004: 


Park et all. l2000l iPetit et al 


20061 ISandratskii et all 120041 ISanvito et all 


2005 


Sato et all l2003l 


van Schilfeaardc and Mrvasovl 


2001 


Wierzbowska et al 


,,2004: Xu et al, 2005)). 



Supercell calculations have usually studied interactions 
between Mn moment orientations by comparing the ener- 
gies of parallel spin and opposite spin orientation states 
in supercells that contain two Mn atoms. An effective 
spin Hamiltonian can be extracted from this approach 
if it is assumed that the interactions are pairwise and 
of Heisenberg form. Even when these assumptions are 
valid, the interaction extracted from these calculations is 
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the sum of the interactions at a set of separations con- 
nected by the superceU lattice vector. If the Mn-Mn spin 
interaction has a range larger than a couple of lattice con- 
stants, this poses a problem for the superceU approach. 
Longer range interactions can, however, be estimated us- 
ing a spin wave approach which allows spin-orientation 
variations that arc incommensurate with the superceU. 
Effective interactions extracted in this way lead to a clas- 
sical Heisenberg model from which the critical tempera- 
ture can usually be calculated withou t substantial fur- 
ther approximation ijSandratskii et "ail, 120041 \Xu et all 
|2005() . The CPA approach can estimate the energy cost 
of flipping a single spin in the ferromagnetic ground state, 
which is proportional to the mean-field-approximation 
for the critical temperature of the effective He isenberg 
model (|Sandratskii et ad 120041: ISato e^ aZ.l.l2003|) . and in 
this sense is limited in its predictive powers when mean- 
field theory is not reliable. Alternatively, a more de- 
tailed picture of magnetic interactions is obtained by di- 
rect mapping of the CPA tota l energy to the Heisenberg 
Ham i ltonian llBer ggvist e t aZ.I. l2004t iKudrnovskv et all . 
l2004HLiechtenstein et aZ.l.ll987D . 

LSDA predictions for spectral properties, like the local 
DOS, are less reliable than predictions for total energy 
related properties. This is especially true for states above 
the Fermi energy, and is manifested by a notorious inac- 
curacy in predicted semiconductor band gaps. From a 
DFT theory point of view, this inconsistency arises from 
attempting to address the physics of quasiparticle excita- 
tions using ground-state DFT. In Mn-doped DMSs, the 
LSDA also fails to account for strong correlations that 
suppress fluctuations in the number of electrons in the 
d-shell. One generally accepted consequence is that the 
energy splitting between the occupied and empty c?-states 
is underestimated in SDF theory, leading to an unrcalis- 
tically large d-statc local DOS near the top of the valence 
band and to an overestimate of the strength of the p — d 
exchange. 

LDA-hU H nisimov et all Il99ll) and self-interaction 
corrected ('SIC'l LSDA llPerdew and Zungeii ll98lD 
schemes have been used to obtain more realis- 
tic energy spectra and help to establish theoreti- 
cally the microscopic origins of ferromagneti sm in 

(III,Mn)V s emiconductor allo y s (iFilippetti et all \200l ; 



Park et alL l2000t iPetit a^J. 120061 ISandratskii al. 



20041: ISchulthess all 1200,^ IShick et all 12004 : 
WierzhowskeTe^ all . 1200'^ . LDA-I-U schemes used 



in studies of (III,Mn)V DMSs combine SDF theory with 
the Hubbard description of strongly correlated localized 
orbitals. Additional parameters from the Hubbard 
model are added to the energy functional which are 
obtained by fitting to experiment or, in principle, can be 
calculated self-consistently (.Anis imov et a/,. . 1991]) . The 
SIC-LSDA method is based on realizing that spurious 
self-interactions present in the SDF lead to unphysically 
large energy penalties for occupying localized states. 
Subtracting these interactions of a particle with itself 
from the density functional suppresses the tendency 



of the LSDA to delocalizc strongly correlated atomic 
orbitals. 



B. Microscopic tight-binding models 

A practical approach that circumvents some of the 
complexities of this strongly-correlated many-body prob- 
lem is base d on the Ande r son many-body Hamilto- 
nian theory ( Anderson', 'l96l':'Flcurov and Kikoin. Il976t 
iHaldane and Anderson,. 1976i:.Krstaiic et al. . 2004) and a 
tigh t-binding-approximation (TEA) band-str ucture the- 
ory (|Harrisonl[l980tlSlater and Kosteilll954|) . The TEA 
Hamiltonian includes the 8x8 sp^ term with second- 
neighbor-i nteraction integrals desc ribing the host semi- 
conductor l)Talwar and Tinelll982(l and terms describing 
hybridization with non-magnetic impurities and Mn. Ef- 
fective single-particle TEA theory is obtained from the 
Anderson Hamiltonian by replacing the density opera- 
tors in the Hubbard term in Eq. ^ with their mean 
values ^MasekL[T99l|) . In the TEA model, local changes 
of the crystal potential at Mn and other impurity sites 
are represented by shifted atomic levels. The parameters 
chosen for the atomic level shifts and the hopping am- 
plitudes between atoms can be inferred from experiment 
in a manner which corrects for some of the limitations 
of LSDA theory. The parameterizati on, summarized in 
l)Masekl [T99lt iTalwar and Tind. Il982|) . provides the cor- 
rect band gap for the host crystal and the appropriate 
exchange splitting of the Mn c?-states. In the calcula- 
tions, the hole density can be varied independently of 
Mn doping by adding non-magnetic donors {e.g. Si or Se 
in GaAs) or acceptors {e.g. C or Ee in GaAs). 

Although the TEA model is a semi-phenomcnological 
theory, it shares with first principles theories the 
advantage of treating disorder microscopically. A 
disadvantage of the tight-binding model approach, 
which is often combined with the CPA, is that it 
normally neglects Coulomb interaction effects which 
influence the charge and spin densities over several 
lattice constants surrounding the Mn ion positions. 
Curie temperatures, magnetizations, the lifetimes of 
Bloch quasiparticle states, the effects of doping and 
disorder on the strength of the p ~ d exchange coupling, 
and the effective Mn-Mn magnetic interaction are 
among the problems th at have been analy zed usi ng this 
tool 



tnc problems tn at nave been analy zed usi ng tni 
Elinowsk i and K acmanl l2003t lJungw irth et al. 
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C. k • p effective Hamiltonian theories 

The highest critical temperatures in (Ga,Mn)As DMSs 
arc achieved in optimally annealed samples and at Mn 
doping levels above 1.5% for which the band holes 
arc itinerant, as evidenced by the metallic conductiv- 
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itics ()Canipion et all l2003a|) . In this regime, semi- 
phenomenological models that are built on crystal Bloch 
states rather than localized basis states for the band 
quasiparticles might be expected to provide more useful 
insights into magnetic and magneto-transport properties. 
A practical approach to this type of modeling starts from 
recognizing that the length scales associated with holes 
in the DMS compounds arc still long enough that a k • p, 
envelope function description of the semiconductor va- 
lence bands is appropriate. Since for many properties it is 
necessary to incorporate spin-orbit coupling in a realistic 
way, six- or eight-band Kohn-Luttinger (KL) k • p Hamil- 
tonians t hat include the spin-orbit split-off band are 
desir able l|Luttinger and KohnLll95,4lVurgaftman et all 
\2nml . 

The kinetic-exchange eff e ctive Hamil t onian approach 
llBhattachariee et Il98,'4 IPietil . 11994 IfTirdynal Il988t 
IZeneil llOSlalT asserts the localized character of the five 
MuQa d-orbitals forming a moment S = 5/2 and de- 
scribes hole states in the valence band using the KL 
Hamiltonian and assuming the p — d exchange interac- 
tion between Mnca and hole spins. As discussed in Sec- 
tion III. B. II the exchange interaction follows from hy- 
bridization between Mn c?-orbitals and valence band p- 
orbitals. The approach implicitly assumes that a canoni- 
cal transformation has been performed \y h ich el ir ninated 
the hybridization l)Schrieffer and Wo3 Il966t iTimmL 
|2003|) . The k • p approximation applies when all relevant 
wavevectors are near the Brillouin-zone center and the 
model also assumes from the outset that the states near 
the Fermi energy mainly have the character of the host 
semiconductor valence band, even in the neighborhood of 
a substitutional Mn. When these assumptions are valid 
it follows from symmetry considerations that the spin- 
dependent part of the effective coupling between Mn and 
band spins is an isotropic Heisenberg interaction char- 
acterized by a single parameter. If the KL Hamiltonian 
parameters are ta ken from the known valu es for the host 
III-V compound ijVurgaftman et a?] . l200l[l . the strength 
of this exchange interaction Jpd can be extracted from 
one set of data, for example from spectroscopic studies 
of isolated Mn acceptors as explained in Se ction [ II.B.l l 
and used t o predict all oth e r properties JPietl et all 
19971 l2000t iJungwirth et all Il999t iKonig et all l2000t 
Matsukura et a/j . ll998|) . Since the value of Jp^ can be ob- 



tained from experiments in a paramagnetic state the ap- 
proach uses no free parameters to model ferromagnetism 
in these systems. In the absence of an external magnetic 
field the KL kinetic-exchange Hamiltonian has the gen- 
eral form: 

H^nhoies + Jpd^Si -StSir.-Rj), (7) 

i,I 

where Tihoies includes the k • p KL Hamiltonian and the 
interactions of holes with the random disorder potential 
and with other holes. The second term in Eq. rep- 
resents the p — d exchange interaction between local Mn 
spins S/ and hole spins s^. 



The k • p approach has the advantage that it fo- 
cuses strongly on the magnetic degrees of freedom in- 
troduced by the dilute moments, which can simplify 
analysis of the model's properties. Disorder can be 
treated in the model by introducing Born approxima- 
tion lifetimes for the Bloch states or by more sophis- 
ticated, cxact-diagonalization or Monte-Carlo methods 
|juiigwirth_ef al, 2002a; K onig et aL .200^ . This ap- 
proach makes it possible to use standard electron-gas 
theory tools to accoun t for h ole-hole Coulomb interac- 
tions ijJungwirth et all Il999() . The envelope function 
approximation is simply extended to model magnetic 
semiconductor heterostructures, like superlattices or 
quan t um wells (|Brev and Guineall2000HFrustaglia et'oll 



2004 iKech rakos et all 120051 iLee et all l2002l l20nC 

Souma et aL ■2005n'l . This strategy will fail, however, 
if the p — d exchange is too strong and the Mn accep- 
tor level is correspondingly too spatially localized or too 
deep in the gap. For example, Mn-doped GaP and GaN 
compounds are likely less favorable for this approach 
than (III,Mn)As and (III,Mn)Sb compounds. Generally 
speaking the advantages of a fully microscopic approach 
have increasing importance for more localized acceptors, 
and hence shorter range Mn-Mn interactions, while the 
advantages of the k • p approach are more clear when the 
acceptors are more shallow and the Mn-Mn interactions 
have longer range. 

D. Impurity band and polaronic models 



There has also been theoretical work on (III,Mn)V 
DMS materials based on still simpler models in which 
holes are assumed to hop between Mn acceptor sites, 
where they interact with the Mn mo ments via phe- 
nome n ological exchange intera cti ons ( [ Alvarez et al 



2002; 



'2002'; 'Bcrciu and Bha tt). looH iBhatt et all 
Chudnovskiv and Pfamikuchel l2002t iFiete et all 
iMavr et all l2002jr " Hamiltonians used in these studies 
have a form of (or similar to), 



2003; 



n 



<ij>,a 



(8) 



where c^^. creates a hole at site i with spin a, the hole 
spin operator aj ~ cl^crapci/j, and eXafj are the Pauli 
matrices. The models apply at least qualitatively in the 
low Mn density limit and are able to directly attack the 
complex and intriguing physics of these unusual insulat- 
ing ferromagnets. Insulating ferromagnetism persists in 
this limit eve n when the carrier density is not strongly 
compensated l|Fiete et a^J . 12001^ . The dependence of fer- 
romagnetism in this regime on the degree of compensa- 
tion has not yet been systematically studied experimen- 
tally and seems certain to pose challenging theoretical 
questions. The models that have been used to study 
ferromagnetism in this regime can easily be adapted to 
include holes that are localized on ionized defects that 
may occur in addition to Mn acceptors. 
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Other related models assume that the Mn acceptors 
are strongly compensated so that the density of lo- 
calized holes is much smaller than the density of Mn 
ions, leading to a polaronic picture in which a single 
hole polarizes a cloud of Mn sp ins ijPurst et all l2002t 
iKaminski and Das Sarmal l2002|^ . The free-parameter 
nature of these phcnomenological models means that 
they have only qualitative predictive power. They are not 
appropriate for the high Tc (Ga.Mn)As materials which 
are heavily doped by weakly compensated Mn acceptors 
and are metallic. On the other hand, the impurity band 
models may represent a useful approach to address ex- 
perimental magnetic and transport properties of ferro- 
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IV. STRUCTURAL PROPERTIES 

A. Impurity formation energies and partial concentrations 

Experimental efforts to increase Mn doping in 
(Ga,Mn)As DMSs beyond the solubiUty limit of 0.1% 
have been assisted by modern ab initio theoretical 
studies of impurity formatio n energies and effects 
related to the growth kinetics (^Edmo nds et al I l2004a[ 
pirwin and Pctukhov, 2002; ..Mahadcvan and Zungei , 
boO.'^lMasek a,Tid MacaL l20oi iMasek et «/J. I2002L 120041). 
Substitution of Ga by Mn is expected, based on 
these studies, to be enhanced when the Ga chemical 
potential is kept low relative to the Mn chemical 
potentia l, i.e., under Ga-poor, Mn-ri ch growth con- 
ditions ijMahadevan and Zungeil 12003(1 . Calculations 
also suggest that one of the major drawbacks of bulk 
growth techniques is that they allow phase separated 
precipitates, such as MnAs, to attain their most stable 
free-standing lattice geometry, leading to relatively 
low formation energies for these unwanted phases 
l|Mahadevan and Zungeil l2003(l . 

In thin film epitaxy, competing phases are forced 
to adopt the crystal structure of the substrate, 
which can significantly increase their formation energies 
(|Mahadevan and Zungei. l2003j) . The non-equilibrium 
LT-MBE has been a particularly successful growth tech- 
nique which allows a synthesis of single-phase (Ga,Mn)As 
DMSs with Mn concentrations up to ~ 10%. If the 
growth temperature is precisely controlled, a 2D growth 
mode of uniform DMSs can be maintained and, at the 
same time, a large number of Mn atom s is incorporated 
in Mnca positions l|Foxon et g^J. 12004)) . 

In (Ga,Mn)As DMSs a significant fraction of the 
Mn atoms is als o incorporated in interstitial positions 
et all l2002|) . Adsorption pathways that can fun- 
nel Mn to interstitial sites have been identified theo- 
retically (see Fig. [HI using ab initio calculations of the 
potential energy surface of M n adsorbed on GaAs(OOl) 
IIErwin and Petukhovl l2002(l . First principles calcula- 




FIG. 8 The Potential-energy surface for Mn adsorption on 
GaAs(OOl), plotted in a plane normal to the surface and con- 
taining the As surface dimer. The minimum energy adsorp- 
tion site is the subsurface interstitial site labeled i; the corre- 
sponding surface geometry is shown (light gray for As, dark 
gray for Ga, yellow for Mn). Typical adsorption pathways 
funnel Mn adatoms to this interstitial site (heavy curves) or 
to a cave site, c (light curves). Inset: Binding energy of a 
Mn adatom centered on the As-dimer; for comparison, re- 
sults are also shown for a Ga adatom. When additional As is 
deposited, the metastable Mn site, s, becomes more favorable 
and leads to partial incorpor ation of substitutional Mn. From 
llErwin and Petukhovl l2002f) . 



tions have also confirmed that interstitial Mnj impurities 
are metastable in GaAs, showing that the three distinct 
positions they can occupy arc two tctrahedral T(As4 or 
Ga4) positions surrounded by four near- neighbor As or 
Ga atoms, and one hexagonal position with three Ga and 
three As nearest neighbors. Among the three intersti- 
tial sites the hexagonal position is clearly less favorable, 
especially so in aii over all p-type (Ga,Mn)As material 
l|Masek and MacaLl2003t) . The typical energy barrier for 
Mn diffu sion between i ntersti t ial sites is approx i matel y 
1 eV (lEdmonds et all l2004al iMasek and Macal |2003() . 
On the other hand, diffusion of Mn between Ga substi- 
tutional positions involves a kick-out mechanism of Mni 
+ Gaoa M uGa + Gai for which the ty pical barrier 
is about 3 eV (|Erwin and Petukhovl |200^ . Interstitial 
Mni is therefore much more mobile than substitutional 

MuGa. 

As mentioned in Section FlLCll Mnj donors are likely 
to form near-neighbor pairs with MuGa acceptors in 
as-grown materials due to the strong Coulomb attrac- 
tion. T he net magnetic moment of su c h a pair is close 
to zero llBlinowski and Kacma 1 l2003t lEdmonds et~d\ . 
I2004al l2005al iMasek and Macal l2003(l . Although Mni 
can be removed by low-temperature annealing, the num- 
ber of substitutional MuGa impurities will remain smaller 
than the total nominal Mn doping. The MuGa dop- 
ing efficiency is, therefore, one of the key parameters 
that may limit the maximum Tc that can be achieved 
in (Ga,Mn)As epilayers, as we discuss in detail in Sec- 
tion IV.AI Ab initio calculations of the formation en- 
ergies can be used to estimate the dependence of Mnca 
and Mni partial c oncentrations on t otal Mn doping in as- 
grown materials l|Jungwirth et all. l2005bt iMasek all 
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l2004l) . Similarly, correlated doping effects can be stud- 
ied for other defects that occur frequently in LT-MBE 

i Ga,Mn)As rnateria ls, such as the arsenic antisites AsQa 
Masek et dYhOO'j} . 

Correlations between acceptors (Mnca) and donors 
(Mni or AsQa) in III-V semiconductors like GaAs are 
strong due to the nearly covalent nature of bonding in 
these crystals. The cohesion energy of the covalent net- 
works has a maximum if the Fermi energy Ep lies within 
a band gap. Whenever Ep is shifted to the valence 
band or the conduction band the strength of the bonds 
is reduced because of the occurrence of unfilled bonding 
states or occupied antibonding states, respectively. 

In case of weak doping, small changes in the impu- 
rity concentration can easily move Ep across the band 
gap with a negligible influence on the energy spectrum. 
The dependence of the formation energy, i.e. the en- 
ergy cost for incorporating a particular impurity in a 
crystal, on the number of electrically active impurities 
can then be represented by the corresponding change 
in Ep multiplied by the charge state of the impurity 
ijMahadevan and Zungeil l2003() . In case of strongly 
doped and mixed crystals, the redistribution of electronic 
states in the valence band due to the impurities may play 
a more important role and should therefore be included 
in the microscopic calculations. 

In general, the formation energy as a function of 
impurity concentrations can be obtained from the 
composition-dependent cohesion energy Wcoh of the crys- 
tal. Assuming a sample consisting of unit cells of the 
impure (mixed) crystal, the formation energy Ea of an 
acceptor A replacing atom X is defined as the reaction 
energy of the substitution process, 

sample -I- A — > sample with one extra acceptor -I- X. 

The corresponding reaction energy is 

Ea{xa,Xd) = N ■ {WtotixA + i/N,XD) ~WtotixA,XD)) 
+ Eatom{,X) — Eatom{,A) , (9) 

where Wtot{xA, xd) is the total energy of the doped crys- 
tal normalized to a unit cell, xa and x d are the acceptor 
and donor concentrations, and the last two terms repre- 
sent total energies of free-standing atoms X and A. With 
increasing size of the sample, TV — > oo, the first term in 
Eq. Q approaches the derivative of Wtot{xA,XD) with 
respect to xa and the atomic energies can be absorbed by 
using the relation between the total and cohesion energy. 
As a result, the cohesion energy represents a generating 
functional for the formation energies, i.e.. 



and xp). 



OEa dEn 



Ea{xa,xd) 



dWcoh{xA,XD) 

dxA 



(10) 



The formation energy of a donor D substituting for an 
atom Y has the same form with D and X ^ Y . 

Having defined Ea and Ep), we note that in the low 
concentration regime where they depend linearly on xa 



dxr 



dxA 



= K{xa,xd), 



(11) 



i.e., that the mutual influence of the two kinds of im- 
purities is symmetric. K(xa, xp) in Eq. (|ll|l plays the 
role of a correlation energy characterizing the co-doping 
process. For positive K(xa, xp), the formation energy of 
one impurity increases in the presence of the other. In 
this case the material tends to be either n-type or p-type 
rather than a compensated semiconductor. On the other 
hand, negative correlation energy indicates that the pres- 
ence of impurities of one kind makes the incorporation of 
the other dopants easier. In this case compensation is 
favored. 

Doping correlations over a wide and continuous 
range of impurity concentrations have been studied us- 
ing the CPA, combined with either the parameter- 
ized TBA model or the ab initio linearized-muffin-tin- 
orbital (LMTO) D FT method jjungwirt h ef al, ,2005^ 
iMasek et all l2004t) . Results of these calculations are 
summarized in Figs. |^ and IIUI The zero of energy is 
set to correspond to the formation energy of a refer- 
ence Gao.96Mno.04As system with all Mn atoms occu- 
pying substitutional Ga positions. Four representative 
examples arc considered here. Scas and Sica are typical 
single donors in GaAs which occupy the anion and cation 
sublatticc, respectively. The As antisite defects AsQa and 
the Mni interstitials, are the most important native de- 
fects in (Ga,Mn)As, both acting as double donors. Fig. 
shows that the formation energy of MuQa decreases with 
an increasing number of donors. The curves are grouped 
into pairs according to the charge state of the donors, 
with only a minor influence of the particular chemical 
origin of the defect. The dependence is almost linear 
for low concentrations and the slope of the function is 
roughly proportional to the charge state of the donor. 
All this indicates that the variation of the formation en- 
ergy of MuQa is mostly determined by the above men- 
tioned Fermi- level effect, and that the redistribution of 
the density of states induced by donor defects plays a 
minor role. Formation energies of the interstitial Mnj in 
the tetrahedral T(As4) position are shown in Fig.^] For 
the donor Mni impurity, the formation energy increases 
with the density of other donors. This means that the 
creation of Mni is efficiently inhibited in the presence of 
AsGa- Analogous results are obtained for Asca antisite 
defects. 

In Fig. Illl we show the change of formation energies of 
AsGa and Mni as a function of the number of MuQa- In 
both cases, the formation energy is a decreasing function 
of the density of MuQa- This self-compensation tendency 
is an important mechanism controlling the properties of 
as-grown (Ga,Mn)As mixed crystals. It explains the ob- 
served charge compensation in as-grown materials and 
is responsible to a large degree for the lattice expansion 
of highly Mn-doped (Ga,Mn)As DMSs, as we discuss in 
Section IrTbI 
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FIG. 9 Theoretical LMTO formation energy £;(MnGa) of sub- 
stitutional MuGa in Gao.96Mno.04As as a function of th e con- 
centration of various donors. From iMasek et aZ.ll2003) . 
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FIG. 10 Theoretical LMTO formation energy £(Mni) of the 
interstitial Mn impurity in Gao.96Mno.04As a function of the 
concentration of various donors. From dMasck et o/.. .2004) . 



The formation energies can be used to theoretically 
estimate partial concentrations of substitutional Mnca, 
Xs, and interstitial Mnj, Xi, in as-grown (Ga,Mn)As ma- 
terials. An assumption is made in these calculations, 
whose validity is tested by a comparison with experimen- 
tal data, that the probabilities of Mn atoms to occupy 
substitutional or interstitial positions are determined by 
the respective formation energies, Es and even in the 
non-equilibrium LT-MBE grown materials. 

The balanced distribution of Mnpa and M i it is re ached 
when l)Jungwirth et all l2005hl: iMasek et all l2002() , 



(12) 
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FIG. 11 Theoretical LMTO variations /^E of the formation 
energies of the interstitial Mn impurity and Asca antisite de- 
fect in Gai_3;Mna;As a function of the concentration x of the 
Mn atoms substituted in the Ga sublattice. The variations 
are measured from t he values fo r the reference material with 
4 % of MuGa From jMasek et a/.. .2004) . 



as also expected from the growth point of view. Partial 
concentrations Xs,i of Mn are obtained by solving Eq. H12() 
together with the condition < < 
X. In Fig. El we show results of TBA/CPA calculations | 



ijJungwirth et all \2005\Jt) : for x > 1.5% Xs ~ 0.8x and 
Xi « 0.2x. LMTO / CPA theory calcu lations give very 
similar predictions jMasek et al\ . \200^ . 

The linear relations between Xs, Xi, and x reflect the 
fact that the difference of the formation energies of MuQa 
and Mni impurities (see inset of Fig. I12|l can be, up to 
X = 10%, approximated by a linear function of Xs and 

^(^Xs^Xi^ Ejgi^X Xi^ EiiXs^Xi^ 

« -0.1 + 5.9a;, - 15.1xi(eV). (13) 



This relation indicates that for x < 1.5%, MuQa has a 
lower formation energy than Mnj and Mn atoms tend to 
occupy substitutional positions. At and above x w 1.5%, 
A(xs,Xi) approaches zero and both Mnca and Mni are 
formed. 

The theoretical results are in a very good agreement 
with experimental data, as shown in Fig. 1121 The balance 
considerations, confirmed experimentally in samples with 
MuGa concentrations up to 6.8%, suggest that there is no 
fundamental physics barrier to increasing Mnca concen- 
tration up to 10% and beyond. Very precise control over 
the growth temperature and stoichiometry is, however, 
required for maintaining the 2D growth mode of the uni- 
form (Ga,Mn)As materials at these high doping levels. 

Finally we note that during growth the formation 
energies control incorporation of Mn atoms, assuming 
that the total amount of Mn in the material is re- 
lated to a sufficiently high chemical potential in the Mn 
source. The annealing processes, on the other hand, do 
not depend on formation energies but rather on energy 
barriers surrounding individual metastable positions of 
Mn in the lattice. The barriers are larger f or M np a 
" Erwin and Petukhovl l2002t iMasek and Macal l2003j) so 
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FIG. 12 Main panel: Theoretical TBA equilibrium partial 
concentrations of substitutional Mnca (red line) and inter- 
stitial Mni (blue line) impurities. Triangles and circles rep- 
resent experimental Mnca and Mni partial concentrations, 
respectively. Inset: Formation energies of Mnpa and Mm as 
a func tion of total Mn concentration. From llJungrwirth et all 
l2005bD . 



that post-growth low temperature annealing can be used 
to remove Mnj without changing the number of Mnca 
significantly. 



B. Lattice constant variation 

Changes in the lattice constant of (Ga,Mn)As DMSs, 
relative to the lattice constant of undoped GaAs, are too 
small to significantly suppress or enhance p — d kinetic 
exchange or other magnetic coupling mechanisms. Di- 
rect effects of doping-induced lattice distortion on the 
onset of ferromagnetism are therefore negligibly small. 
Nevertheless, variations in the lattice parameter provide 
a measure of impurity concentrations in the DMS ma- 
terial. The impurities do, of course, control ferromag- 
netism through their doping properties. Also, because 
thin film (Ga,Mn)As epilayers are not relaxed, lattice 
constant mismatch between the DMS layer and the sub- 
strate induces strains that in many cases determine mag- 
netocrystalline and magnetotransport anisotropics, as we 
discuss in detail in Sections IV. C.ll and IVI.lH Theoreti- 
cal calculations (|Masek et aL, . .200S) of the dependence of 
the lattice constant on the density of the most common 
impurities in DMSs represent, therefore, another piece 
to the mosaic of our understanding of ferromagnetism in 
these complex systems. 

Considering the values of the atomic radii of Mn {Run 
= 1.17 A) and Ga (i?Ga = 1-25 A), the substitutional 
MuGa impurity may be expected to lead to only very 
small changes (reductions) in the lattice constant. This 



expec tation is consistent with the calculated ijZhao et all 
I2002^ lattice constant of a hypothetical zinc-blende MnAs 
crystal whose value is comparable to that of GaAs. On 
the other han d, As antisites produce an expansion of the 
GaAs lattice (|Liu all 1 19951 IStaab all l200lt) and a 
similar trend can be expected for interstitial Mn. 

Modern density-functional techniques allow one to 
move beyond intuitive theoretical considerations and 
discuss the dependence of the lattice constant on im- 
purity conce ntrations on a more quantitative level 
(|Masek et ail . |2003 ). The CPA is again a useful tool 
here for studying (Ga,Mn)As properties over a wide 
range of impurity concentrations. Some quantitative in- 
accuracies in theoretical results due to the limitations 
of the LMTO/CPA approach {e.g. neglect of local 
lattice relaxations) have been corrected by using the 
full-pote ntial linearized- augment ed-plane-wave supercell 
method ^Masek and MacaL l2005(l . Starting with an ideal 
(Ga,Mn)As mixed crystal with all Mn atoms occupying 
substitutional Ga positions, these calculations give the 
following Vegard's law expression for the doping depen- 
dence of the lattice constant: 



aicixs) = aa + asXs{A), 



(14) 



with the expansion coefficient Os ranging from -0.05 to 
0.02 depending on t he method used in the calculation 
l)Masek ad l2003t iMasek and Macal l2005fl . As ex- 
pected, a changes only weakly with the MuQa density 

A similar linear dependence is obtained for hypo- 
thetical crystals where AsQa (or Mnj) is the only im- 
purity present in the material, as shown in Fig. 1131 
According t o the more relia ble full-potential supercell 
calculations ijMasek and Maca. 2005.) , the composition- 
dependent lattice constant is found to obey: 

y)^ao- O.QdXs + 0.48a:, + 0.46j/(A), (15) 

where Xi and y are the densities of Mnj and AsQa, re- 
spectively, and ao is the lattice constant of pure GaAs. 

Recently^ several expe ri mental work s 

Ikurvliszvn-Kudelska et HI. I2004t J Potashnik et al. 
l2002t ISadowski and Domaealal l200i IZhao et all 1200 
have studied the dependence of lattice constants in 
(Ga,Mn)As materials on disorder, based on the com- 
parison between as-grown and annealed samples. The 
measurements confirmed that both AsQa and Mni 
defects lead to a significant expansion of the lattice. In 
samples grown under As-rich condition, which is ex- 
pected to inhibit formation of Mnj impurities, annealing 
has virtually no effect on the measured lattice constant. 
This is consistent with the stability of AsQa defects 
up to temperatures that are far above the annealing 
temperatures. Mni impurities, on the other hand, can 
be efficiently removed by low-temperature annealing. 
Consistently, annealing leads to a significant reduction 
of the lattice constant in materials which contain a large 
number of these defects in the as-grown form, as shown 
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FIG. 13 Theoretical LMTO/CPA lattice constant of 
(Ga,Mn)As as a function of the concentration of the most im- 
portant impurities Mnoa (circ les), Mni (full triangles) , and 
AsGa (empty triangles) . From iMasek and Macal . l2005fl . 



in Fig. 1141 On a quantitative level, experimental data 
suggest a stronger lattice expansion due to Mnca and 
a weaker expansion due to Asca and Mnj, compared to 
the theoretical predictions of Eq. IjlSfl . The quantitative 
disagreement can be attributed, in part, to the simplified 
description of the system within the theoretical model. 
Also, the presence of other lattice imperfections or in- 
accuracies in the determination of experimental doping 
values may have partly obscured the direct quantitative 
comparison between experiment and theory. 



V. MAGNETIC PROPERTIES 

A. Ferromagnetic critical temperature 

Curie temperatures in metallic (Ga,Mn)As 
have been studied theoretically starting from 
the k-p kinetic-excha ng e eff ecti ve Hamiltonian 
llBrey and Gomez-Santosl. l2003t iDas Sarina et al 



2004 iDietl et all I2000t1jungwirth et all Il999l 120021 , 



2005b() and from micr oscopic TBA o r SD F band 
structure calculations ( Berggvist et all \200 



B^zera^^^^ 
^mgwir^^^^ 



2005a 



2003a 



iHilbert and Noltin„, 

2005bt ISandratskii and Brunc , 



200 
2005t 



2002t ISa ndratskii a d 120041: \Sato et all I200S : 
Timm and Ma cDonaldl l2005t IXu et all \2003i} . For 

a more detail description of these theoretical approaches 
see Sections IIII.AIIII.CI The advantage of the k-p 
kinetic-exchange model is that it uses the experimen- 
tal value for the p — d coupling constant Jpd, i.e., 
it correctly captures the strength of the magnetic 
interaction that has been established to play the key 
role in ferromagnetism in (Ga,Mn)As. The model also 
accounts for strong spin-orbit interaction present in 
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FIG. 14 Relaxed lattice constants for the 50-nm-thick 
(Ga,Mn)As films before annealing (squares) and after anneal- 
ing (triangles) as a function of substitutional Mnca content. 
The inset shows the change of the relaxed lattice constant as 
a function of the change in the fraction of interstitial Mn due 
to the out-diffusion of interstitial Mn during annealing (i.e., 
Mni,i-Mni,f). From l|Zhao pi all\200^ . 



the host valence band which splits the three p-bands 
into a heavy-hole, light-hole, and a split-off band with 
different dispersions. The spin-orbit coupling was sho wn 
l)Brev and Gomez-Santol l2003t lK5nig et all l2001a|) to 
play an important role in suppressing magnetization 
fluctuation effects and, therefore, stabilizing the ferro- 
magnetic state up to high temperatures. On the other 
hand, describing the potentially complex behavior of 
Mnca in GaAs by a single parameter may oversimplify 
the problem. The calculations omit, for example, the 
suppression of Tc in low hole-density (Ga,Mn)As mate- 
rials due to the direct antiferromagnetic superexchange 
contribution to the coupling of near-neighbor Mn pairs. 
The whole model inevitably breaks down in DMS 
systems with holes strongly bound to Mn acceptors or 
with large charge fluctuations on Mnca d-shclls. 

The advantage of microscopic approaches to Curie 
temperature calculations is that they make no assump- 
tion about the character of MuQa impurities in GaAs and 
their magnetic coupling. They are therefore useful for 
studying material trends in Tc as a function of Mn dop- 
ing or the density of other intentional or unintentional 
impurities and defects present in real systems. Because 
spin-orbit interactions add to the numerical complexity 
of calculations that are already challenging, they have 
normally been neglected in this approach. Another short- 
coming, discussed already in Section lill.AI of the LSDA 
approaches is an overestimated strength of the p — d ex- 
change as compared to experiment. Within the mean- 
field approximation, which considers thermodynamics of 
an isolated Mn moment in an effective field and neglects 
correlated Mn-Mn fluctuations, microscopic calculations 
typically yield larger T^s than the effective Hamiltonian 
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model that uses the experimental value for Jpd- Stronger 
p—d exchange and the omission of spin-orbit coupling ef- 
fects in the DFT theories, however, also leads to a larger 
suppression of the Curie temperature due to fluctuation 
effects. (A closer agreement on the character of the Tc 
versus Mn-doping curves, calculated within the two for- 
malisms, is obtained when the deficiencies of LSDA the- 
ories are partly eliminated by using, e.g., the LDA+U.) 
Despite the above weaknesses of scmi-phenomenological 
and microscopic calculations, a qualitatively consistent 
picture is clearly emerging from these complementary 
theoretical approaches that, as we discuss below, pro- 
vides a useful framework for analyzing measured Tc's. 



1. Mean-field theory in (Ga,l\/ln)As 

Our review of theoretical Tc trends in (Ga,Mn)As 
starts with the results of the KL kinetic- exchange Hamil- 
tonian and mean-field approximation to set up a scale 
of expected Curie temperatures. These estimates which 
are not accurate in all regimes, are simplified by as- 
suming a homogeneous distribution of MuQa ions and 
neglecting the role of other defects, apart from their 
poten tial contrib u tion to hole or moment compensa- 
tion llDietl aZ.l Il997l 120001: lJungwirth et all Il999t 
iKonig et all I20n3^ . In the case of microscopic mod- 
els this assumption is equivalent to the virtual crys- 
tal approximation. Microscopic TBA calculations have 
shown very little effect of positional disorder on the 
strength of magnetic couplings in (Ga;Mn)As cpilayers 
with metallic con ductivities, partly justifyi ng the virtual- 
crystal approach jJungwirth et a/J . l2003a(l . In addition, 
detailed theoretical studies confirm the absence of any 
significant magnetic frustration associated with the ran- 
dom positions of Mnca moments in the l attice in the 
more metallic ferromagnetic semiconductors (jFiete et alV 
|2005; Timm and MacDonald. 2005}. In the very dilute 
limit, however, Tc becomes sensi tive to the dis tribution 
of the Mn moments in the lattice ravist_^^ 



_B, around (S) 



In the mean-field approximation ijPietl et al. 



I20n4^ 



lJungwirth et a/J . |l999j) . each local Mnca moment is de- 
scribed by a Hamiltonian S/-Hm_f where S/ is the Mnca 
local spin operator, Hmf = Jpd{s}, and (s) is the mean 
spin density of the valence band holes (for the definition 
of the Jpd field see Scction ril.B.lll . Hmf is an effective 
field seen by the local moments due to spin-polarization 
of the band holes, analogous to the nuclear Knight shift. 
Similarly Iimf ~ JpdNMn{S) is an effective magnetic 
field experienced by the valence band holes which is pro- 
portional to the density and the mean spin polarization 
of the MuGa local moments. The de pendence of (S) on 
temperature and field H m p is given ij Koni^ et all l2003l) 
by the Brillouin function (jAshcrofriin^MemfiirTl^T^i : 



H 



MF 



SBs{S\HMF\/kBT) 



(16) 



\Hmf\ 

The Curie temperature is found by linearizing Hmp and 



= 0: 



S+IS\Hmf\ 
3 knTr 



(17) 



Here X/ is the itinerant hole spin susceptibility given by 
d{s) d'^CT 



Xf 



dh 



MF 



"'"'MF 



(18) 



and ct is the total energy per volume of the holes. 
Eqs. and ((HI) give 



ksTc 



NunS{S + 1] 



JpdXf 



(19) 



The qualitative implications of this Tc-equation can be 
understood within a model itinerant hole system with a 
single spin-split band and an effective mass rn*. The 
kinetic energy contribution, e^, to the total energy of the 
band holes gives a susceptibility: 



Xf,k 



dh^ 

"'"'MF 



m*kp 
4^ 



(20) 



where kp is the Fermi wavevector. Within this approx- 
imation Tc is proportional to the Mnga local moment 
density, to the hole Fermi wavevector, i.e. to p^^^ where 
p is the hole density, and to the hole effective mass m*. 

A more quantitative prediction for the Curie tempera- 
ture is obtained by evaluating the itinerant hole suscep- 
tibility using a realistic band Hamiltonian, 



n = HKL + s-h 



AIF, 



(21) 



where Hkl is the six-band KL Hamiltonian of the GaAs 
host band llVurgaftman et all. l200ll) and s is the hole 
spin operator l|Abolfath et ai!.Ll2001at iDietl et ail l2001b[ 

I2000D . The results, represented by the solid black line in 
Fig. El E^re consistent with the qualitative analysis based 
on the parabolic band model, i.e., Tc roughly follows the 
~ 2;pi/3 dependence. Based on these calculations, room 
temperature ferromagnetism in (Ga,Mn)As is expected 
for 10% Mnpa, doping i n weakly compensated samples 
(lJungwirth erg/J . 12005^1) . 

Hole-hole Coulomb interaction effects can be included, 
in the lowest order of perturbation theory by ad ding the 
hole exchange contribution to the total energy l)Mahari 
Il98l|) . The red hne in Fig. UHl shows this Stoner Tc 
enhancement calculated numerically using the kinetic- 
exchange effective model with the six-band KL Hamilto- 
nian. Tc stays roughly proportional to xp^^^ even if hole- 
hole exchange interactions are included, and the enhance- 
ment of Tc due to interactions is of the order ~ 10 — 20% 
(|Dietl et a/J . l2001btljimgwirth et g/l Il999ll2005b|) . 

The mean-field effective Hamiltonian analysis above 
neglects discreteness in the random Mnca positions in the 
lattice and magnetic coupling mechanisms additional to 
the kinetic exchange contribution, particularly the near 
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FIG. 15 Ferromagnetic transition temperatures of 
(Ga,Mn)As calculated within the effective Hamiltonian 
and virtual crystal approximation: mean- field (thick black 
line), Stoner enhancement of Tc (thin red line), Tc af- 
ter correcting for correlated Mn orientation fluctuations 
using a spin-wave appr oximation (blue symbols). From 
iTungwirth e.t a/.Ll2nn5HV 



as the TBA/CPA (jjungwirth el all l2005b(l or ah initio 
appro aches based on either CPA or supercell ba nd struc- 
tures (jSandratskii et a;J . l2004ISato aZ.l.f2003t) . capture 
all these effects on an equal footing and can be used to es- 
timate trends in mean-field Tc beyond the virtual crystal 
approximation. 
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neighbor superexchange. The former point can be ex- 
pected to influence Tc at large hole densities, i.e., when 
the hole Fermi wavelength approaches inter-atomic dis- 
tances. Of course, the entire phenomenological scheme 
fails on many fronts when the Fermi wavelength ap- 
proaches atomic length scales since it is motivated by 
the assumption that all relevant length scales arc long; 
the k-p band structure, the use of the host material band 
parameters, and the neglect of momentum dependence in 
the Jpd parameter all become less reliable as the hole den- 
sity increases to very large values. The approximations 
are apparently not fatal, however, even for x ^ 10% and 
any degree of compensation. 

In the opposite limit of strongly compensated sys- 
tems, where the overall magnitude of the hole-mediated 
exchange is weaker, antiferromagnetic superexchange 
can dominate the near- n eighbor MnQa-Mnca coupling 
(|Kudrnovskv al\. l2004^■ leading to a reduced Curie 
temp erature ( Jungwirth et a/j . l2005bl : ISandratskii et all 
l2004j) . We emphasize that the k-p kinetic-exchange 
model cannot be applied consistently when nearest neigh- 
bor interactions dominate, since it implicitly assumes 
that all length scales are longer than a lattice con- 
stant. We also note that net antiferromagnetic coupling 
of near-neighbor MuQa-MuQa pairs is expected only in 
systems with large charge compensations. In weakly 
compensat ed (Ga,Mn)As the ferromagnetic contribution 
takes over (iDietl et g^J. l2001bt iKudrnovskv et"all.l2004t 
iMahadevan and ZnnseA 12004) 

In addition to the above effects related to random Mn 
distribution, Mn positional disorder can directly modify 
p — d interactions when the coherence of Bloch states be- 
comes significantly disturbed. Microscopic theories, such 



FIG. 16 Tc calculations within the microscopic TBA/CPA 
model: Tc versus hole density (left panel), Tc versus number 
of holes per Mnca (right panel). The overall theoretical T c 
trend is highlighted in grey. From (I Jungwirth et aZll2005tf) . 

The mean-field CPA Curie temperatures are obtained 
by evaluating the energy cost of flipping one Mnca mo- 
ment with all other moments held fixed in the ferromag- 
netic ground state. I t can be evaluated for any given 
chem ical composition (jLiechtenstein et aZ.l.ll987tlMasekl 
I1991D and defines an effective exchange field H^f / acting 
on the local moment. This energy change corresponds to 
Hmf in the kinetic-exchange model used in the previous 
section, i.e., 

ksTc = ^Hcff . (22) 

Results based on microscopic TBA band structure cal- 
culations are shown in Fig. 1161 as a function of hole 
density for sever a l Mnp a local moment concentrations 
(|Jungwirth et all l2005b|) . (The hole density is varied 
independently of Mnca doping in these calculations by 
adding non-magnetic donors or acceptors.) Comparison 
with Fig. El identifies the main physical origins of the 
deviations from the Tc ^ xp^^^ trend. Black dots in the 
left panel of Fig. E| correspond to a relatively low local 
MuGa moment concentration {x ~ 2%) and hole den- 
sities ranging up to p = 47Vjv/„ and show the expected 
suppression of Tc at large p. The effect of superexchange 
in the opposite limit is clearly seen when inspecting, e.g., 
the X = 10% data for p < 1 nm""^. The mean-field 
TBA/CPA Curie temperature is largely suppressed here 
or even negative, meaning that the ferromagnetic state 
becomes unstable due to the short-range antifcrromag- 
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nctic coupling. Note, that the neglect of Coulomb inter- 
actions in these TB A/CPA calculations likely leads to an 
overestimated strength of the antiferromagnetic superex- 
change. The inhomogeneity of the carrier distribution in 
the disordered mixed crystal may also contribute to the 
steep decrease of Tc with increasing compensation seen 
in Fig. Hi 

Although the Curie temperatures in the left panel of 
Fig. ^] appear to depart strongly for the Tc ^ xp^^^ 
dependence, the linearity with x is almost fully recov- 
ered when Tc is plotted as a function of the number of 
holes per MuQa local moment, p/Nmu (see right panel of 
Fig. [TB|l . Note that for compensations {I—p/Nmu) reach- 
ing 100% this property of the superexchange coupling is 
reminiscent of the be havior of (II,M n)VI diluted mag- 
netic semiconductors (|Furdvnal . Il988(l in which Mn acts 
as an isovalent magnetic impurity. The dependence on 
p in (Ga,Mn)As is expected to become very weak, how- 
ever, when approaching the uncompensated state. Sim- 
ilarly, the prospects for substantial increases in Tc by 
non-magnetic acceptor co-doping of weakly compensated 
material appear to be quite limited. 

In the left panel of Fig. El we show mean-field CPA 

Curie temperatures in uncompensated {p/Nmu = 1) 

(Ga,Mn)As DMSs as a function of Mn doping calculated 

using LDA and LDA+U ab initio methods. The LDA+U 

calculations, which give more realistic values of the p — d 

exchange coupling, confirm the linear dependence of Tc 

on x, showing no signs of saturation even at the largest 

doping X = 10% considered in these calculations. 
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FIG. 17 LDA/CPA (red) and LDA-hU/CPA (blue) calcu- 
lations of the Curie temperature as a function of Mn dop- 
ing in uncompen sated (one hole per Mn ) (Ga,Mn)As and 
(Ga,Mn)N. From iSandratskii et ai.ll200ll . 



2. Role of collective Mn-moments fluctuations and different 
II I- V hosts 

The potential influence of correlated Mn-moment fluc- 
tuations (corrections to mean field theory) on ferromag- 
netic ordering in (Ga,Mn)As can be recognized by consid- 
ering, within a simple parabolic band model, the RKKY 
oscillations effect which occurs as a consequence of the 
2kF anomaly in the wa yevector dependent susceptibil- 
ity of the hole system ijBrev and Gomez-Santod l2003t 



iDietl "oZI . Il997jl . In this theory, which treats the hole 
system perturbatively around the paramagnetic state, 
the sign of the hole-mediated Mnca-Mnca coupling varies 
as cos(2A:i?d), where d is the distance between Mnca 
moments, and its amplitude decays as d^ . Estimating 
the average Mnca-Mnca separation in a (Ga,Mn)As ran- 
dom alloy as J = 2(3/47riVj\/„)^/'^ for uncompensated 
(Ga,Mn)As systems and neglecting spin-orbit coupling 
and band warping, cos{2kFd) ~ —1, which means that 
the role of the RKKY oscillations cannot be generally 
discarded. In realistic valence bands the oscillations 
are suppressed due to non-parabolic and anisotropic dis- 
persions of the heavy- and li ght-hole bands and due 
to st rong spin-orbit coup ling (|Brev and Gomez-SantosL 
l200aiKonig- fi^^. 12001 ak 

More quantitatively, the range of reliability and cor- 
rections to the mean-field approximation in (Ga,Mn)As 
can be estimated by accounting for the suppression of 
the Curie temperature using the quantum theory of long- 
wavelength spin- waves or using Monte Carlo simulations 
which treat Mn-moments as classical variables. For weak 
p — d exchange coupling, SNMnJpd/ Ep *C 1 where Ep is 
the hole Fermi energy, the spin-polarization of the hole 
system is small and the RKKY and spin- wave approxima- 
tions treat collective Mn-moment fluctuations on a simi- 
lar level. The advantage of the spin-wave theory is that 
it can be used to explore the robustness of ferromagnetic 
states over a wider range oi p — d couplings, including 
the more strongly exchange-coupled as-grown materials 
with large Mn density and large hole compensation (See 
Section rVlILAI for a general discussion of magnetic inter- 
actions in the two coupling-strength limits). 

Calculations in metallic systems have been performed 



from the k-p kinetic-exchang 


e effective Hamil- 


jBrev and Gomez-Santosl 


1200,'^ 


iJunEwirth et al\. 



1 I 

band structur e calcula ti ons JB erggvist et al, 2004 



Bouzerar et all l2005at iHilbei't and Noltine, .2005 



Xu et all l2005(l . (Monte Carlo studies of T, 



the regime near the metal-insulator are reported in 
iMayr.eL^,|2Q0l.) Within a non-interacting spin- wave 
approximation, the magnetization vanishes at the tem- 
perature where the number of excited spin waves equals 
the total spi n of the ground state l)Jungwirth et all 
12002a I2005H) : 



knTr 



2S+1 
6 



klD{Tc) , 



(23) 



where ko = (BTr^iVMn)^/^ is the Debye cutoff and D(T) = 
Dq (S) (T) / S is proportional to the zero-temperature 
spin-wave stiffness parameter A {Dq — 2A/ SNmu) 
and the mean-field temperature-dependent averag e spin 
on Mn, {S){T) l|K5nig e t al. 2001a; Schlicniann et all 
l2001al) . (For a detail discussion of the micromagnetic 
parameter A see Section FV-CII . 

Comparing the spin-stiffness results obtained using 
the k-p kinetic-exchange model with a simple parabolic 
band and with the more realistic, spin-orbit coupled KL 
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Hamiltonian, the spin stiffne ss is observed to always b e 
much larger in the KL model ijKonig et adl2001all2000|) . 
For (Ga,Mn)As, the parabolic-band model underesti- 
mates Z? by a factor of ~10-30 for typical hole densities. 
This larger spin stiffness in the spin-orbit coupled va- 
lence bands is due to the heavy-hole - light-hole mixing. 
Crudely, the large mass heavy-hole band dominates the 
spin susceptibility and enables local (mean-field) mag- 
netic order at high temperatures, while the dispersive 
light-hole band dominates the spin stiffness and enables 
long-range magnetic order. The analysis highlights that 
the multi-band character of the semiconductor valence 
band plays an important role in the ferromagnetism of 
(Ga,Mn)As. 

Critical temperature estimates based on Eq. 123|) . the 
KL kinetic-exchange model, and including also the Stoner 
enhancement of Tc are summarized in Fig. llSl bv the blue 
symbols. These Tc estimates indicate that Tc will remain 
roughly proportional to x even at large dopings. The sup- 
pression of Tc due to spin- waves increases with increasing 
hole density relative to the local moment concentration, 
resulting in saturation of the critical temperature with 
increasing p at about 50% compensation. 

The suppression of Tc due to correlated Mn-moment 
fluctuations is a l so observed i n the LSDA ca l culations 
iBergnvist all EM 12003 iRouzerar et all MEM 
iHilbert and Noltind. l2005t IXu et all |2005|) . The trend 
is illustrated in Fig. ^| where collective fluctua- 
tions are accounted for using the spin-wave t heory 
or the Monte Carlo app roach (|Bergavist et all l2004t 
iHilbert and Noltind . l2005|) : similar trends of suppressed 
mean-field Tc due to collective Mn-moment fiuctuations 
have been predicted by a spin-wave theory using a more 
elab orate, self-consistent RP A technique on random lat- 
tice (|Bouzerar et a?J . l2005a|) . A larger suppression of the 
mean- field Tc seen in the ab initio calculations, com- 
pared to the KL kinetic-exchange model results, can be 
attributed partly to the simpler, three-fold degenerate 
LSDA valence band structure in theories that neglect 
spin-orbit coupling. Also, the stronger p — d exchange in 
the LSDA theories may result in a weaker spin stiffness 
of the magnetic system, as the holes are more strongly 
bound to the Mn acceptors and the hole mediated Mn- 
Mn coupling has a more short-range character. The en- 
hancement of fluctuation effects in stronger p ~ d cou- 
pled systems is clearly seen in Fig. 1181 when comparing 
the LSDA results for narrower gap (weaker p — d ex- 
change) (Ga,Mn)As and wider gap (stronger p — d ex- 
change) (Ga,Mn)N. 

The quantitative discrepancies between KL kinetic- 
exchange model and LSDA results for the mean-field Tc 
and for the suppression of ferromagnetism due to collec- 
tive Mn-moment fiuctuations partly cancel out, leading 
to similar overall predictions of the Curie temperatures in 
(Ga,Mn)As. Based on the Tc analysis alone it is therefore 
difficult to determine whether the magnetic interactions 
have a short-range or long-range character in (Ga,Mn)As 
DMSs. 
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FIG. 18 Comparison of the ab initio Curie temperatures 
of Gai-iMn^As (diamonds) and Gai-iMn^N (squares) ob- 
tained by the spin-wave approach (solid line, filled symbols), 
virtual-crystal-approximation-RPA (dashed line, filled sym- 
bols), mean-field-approximation (dotted line, filled symbols) 
and Monte Carlo (dash-do tted line, open symbols). From 
jHilhert a,nH Noltind 



Theoretically, the localization of the hole around the 
Mn impurity and the range of magnetic Mn-Mn in- 
teractions can be studied using microscopic TBA or 
ab initio calculations of the charge and moment dis- 
tribution in the lattice or by mapping the total en- 
ergy of the DMS cry st al to t he Heisenberg Ham il tonian 
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anion column in the periodic table from the nitride 
DMSs to antimonides, th e holes become more delocalized 
ijMahadevan and Zungeil l2004|) and, consequently, the 
Mn-Mn interactions are more long-range in these micro- 
scopic calculations. In (Ga,Mn)As the LSDA theory pre- 
dicts short-range magnetic coupling while the LDA+U 
or SIC-LSDA results suggest that the holes which medi- 
ate the Mn-Mn exchange interaction are more delocalized 
llvan Schilfgaarde and Mrvasov , l200lt ISchulthess et all 
l2005HWierzbowska et fflZ.l . l2004D . 

Combined theoretical and experimental studies of re- 
manent magnetization, micromagnetic parameters, and 
magneto-transport coefficients, discussed in detail in Sec- 
tions ^^^^1 indicate that in high quality (Ga.Mn)As 
fcrromagnets with metallic conductivities (conductivity 
increases with decreasing temperature) the holes are suf- 
ficiently delocalized to make the kinetic-exchange model 
approach applicable. It is natural to expect that the free- 
carrier-mediated ferromagnetism picture will also apply 
in narrower gap antimonidc DMSs, such as (In,Mn)Sb, 
with even large r conductivities due to the smaller hole 
effective mass ijWoitowicz et all . l2003|) . Smaller cffcc- 
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0.1 
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64 
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0.1 
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85 


88 
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85 


InSb 


0.1 


11 


12 


11 


11 
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37 
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TABLE I Mean-field {Tf^), Stoner enhanced (T^), spin- 
wave suppressed (T^°"), and estimated including both Stoner 
enhancement and spin- wave suppression {T°^^) ferromagnetic 
transition temperatures in III-V host semiconductors doped 
with 5% of Mn and with itinerant hole densit ies p = 0.1 and 
p = 0.5 nm-3. From JJungwirth et a^.U2002b^ . 
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FIG. 19 Temperature dependence of remanent mag- 
netization and inverse paramagnetic susceptibility for 
Gao.91Mno.09As sample with Tc = 173 K; inset: hysteresi s 
loop for the same sample at 172 K. From (IWang et aZ.U2005a^ . 



tive mass and larger unit cell volume, as compared to 
arsenide DMSs, explain the smaller Curie temperatures 
in (III,Mn)Sb. This trend is illustrated in Table |I| by 
comparing the respective mean-field KL kinetic-exchange 
model Tc's. In Mn-dopcd phosphides and nitrides the 
suppression of due to effects beyond the mean-field 
virtual-crystal approximation, seen in the LSDA calcu- 
lations and related to the shor t-range nature of mag - 
netic interactions, may prevail ijScarpuUa et al\ l2005j) . 
We note also that the above LSDA Curie temperature 
studies do not capture the possible transition of the Mn 
impurity state in the wide gap III-V's to the highly cor- 
related tri-valent (d*) center with four strongly local- 
ized d-clcctrons and an empty d-state deep in the gap, 
in which case the hole mediated fcrromag netism picture 
these calculations imp l y is iiot applicable llKreissl et alV 
I1996HLuo and Martini I200I ISchulthess et 0/11200511 



Experimentally, trends ha ve been mos t exten- 
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sively studied in (Ga,Mn)As 


IChiba et a/.l 


12003a 


Edmonds et al\. l2002al lJungwirth et alV 


12005b 


Ku et alV I200r4 lOhnd. Il998l 


Potashnik etall 


I2OO2 


Stone et alV 120031 lYu et oi.l I2003D. In Fig. 


d we 



show the temperature-dependent magnetization and 
inverse susceptibility of the current record Tc material 
(|Wang et all l2005aj) . The Brillouin- function character 
of the magnetization curve confirms that the mean-field 
theory is appropriate in these high quality DMS materi- 
als with metallic conductivities. Curie temperatures for 
a series of as-grown and annealed (Ga,Mn)As samples 
with experimentally characterized cha rge and moment 
compe nsations are plotted in Fig. |2D1 ijJungwirth et all 
l2005b() . The concentration of uncompensated Mn 
moments in the plot is x^ff = Xg — Xi, where it is 
assumed that the Mnj donors present in the system are 
attracted to MuQa acceptors and that these pairs couples 
antiferromagnetically, as discussed in Sections |II.C. H and 
IIV.AI (The consistency of this assumption is confirmed 



by independent magnetization studies reviewed in 
Section lV.Bp . The experimental Tc/xeff plotted against 
p/NIJjI in Fig. 1201 where 7V^/;[ = Axeff/af^, show a 
common Tc trend whi ch is co n sistent with theo retical 
expectations iBergqvist e^^. 120051: iBouzerar et all 
l2005at lJungwirth et an . l2005b() . In particular, theory 
and experiment agree on the very weak dependence 
of Tc/xcff on p/N'^^^ for low compensation and the 
relatively rapid fall of Tc/xeff with decreasing p/N^^^ 
for compensations of ^ 40% or larger. It should be noted 
that the maximum experimental Xeff is only 4.6% in the 
as grown sample and 6.8% after annealing for a total Mn 
concentration x ~ 9%. Hence the modest Tc's observed 
so far. Achieving Tc values close to room temperature in 
(Ga,Mn)As, which is expected to occur for Xeff ~ 10%, 
appears to be essential ly a material growth issu e, albeit 
a very challenging one ij.Iungwirth fl/ll2005lJ ^. 

Only a few experimental studies of LT- MBE grown 
(III,Mn)Sb DMSs have been r eported to date ijAbe et all 
l2000t IWoitowicz et oil l2003(l . Curie temperatures mea- 
sured in these materials are lower than Tc's in Mn-doped 
arsenides, consistent with the kinetic-exchange model 
predictions in Table The nature of fcrromagnetism 
and, therefore, the interpretation of experimental Curie 
temperatures observed in phosphide and nitride DMSs 
arc not established yet, as we have already pointed out 
in Section ILBI 



B. Magnetization 

In this section we focus on low-temperature ferro- 
magnetic moments in (Ga,Mn)As DMSs with metal- 
lic conductivities. Early experimental studies, report- 
ing large apparen t magn e tization deficits in (Ga , Mn)As 
([Korzhavvi et gj I l2002t lOh no and Matsukural l200lt 
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FIG. 20 Experimental T^/Xeff vs. hole density relative to 
effective concentration of Mnca moments. Deviations from 
linear dependence on Xeff are seen only for high compensa- 
tions (1 — paf^/Axeff = 1 — p/NlfJ^ > 40%) in agreement 
with theory. For weakly compensated samples Tc shows no 
signs of saturation with increasing x^ff. Theoretical (grey) 
Tc trend from F ig. 1161 is plotted for comparison. From 



I2002D . 

Some degree of non-coUinearity is inevitable as a com- 
bined consequence of positional disorder and spin-orbit 
coupling. Nevertheless it was argued theoretically that 
a large suppression of the ferromagnetic moment is not 
expected in metallic (Ga ,Mn)As samples with Mn con - 
centrations above ~ 1% (jTimm and MacDonaldl f2005(l . 
The minor role of non-coUinearity is due largely to 
the long-range character of magnetic interactions, which 
tends to average out the frustrating effect of anisotropic 
coupling between randomly d i stributed Mn imp urities 
IITimm and MacDonaldl l2005t IZhou et all \2004\ . In- 
deed, ab initio, microscopic TBA, and k-p kinetic- 
exchange model calculations of zero-temperature mag- 
netic moments in (Ga,Mn)As ferromagnets which neglect 
effects that would lead t o non-coUinearitv llDietl et all 
2001bt lJungwirth et al\. l2005al ISchulthess et all l2005t 
Wierzbowskae^fflZj . |2004|) are consistent with experi- 



llJungwirth et aLll2005^ 



13). 



ments reporte d in a series of high-quality (Ga,Mn)As 
ferroinagnets llEdmonds et all l2005al: lJungwirth et all 
l2005at IWang et all 12005^^ " It rules out any marked in- 
trinsic frustrations in the ground state of these DMSs. 
The substantial magnetization suppression seen in many 
early (Ga,Mn)As samples can be attributed primarily to 
the role played in those samples by interstitial Mn atoms 
and other unintentional defects. 



iPotashnik et all l2002() . motivated a theoretical search 
for possible intrinsic origins of frustrating magnetic 
interactions in this material. Using a wide spec- 
trum of com putational techniques, ranging from ab ini- 
tio methods Ip^o rzhavYj- ej_a4 120021: iKudrnovskv et all . 
120041: iMahadevan et aZJ. l2004ll and microscopic TBA 
ijTimm a n^^acDonaldf l2005tl to k -p kin e tic- exchange 
mode l s llBrev and Gomez-Santosl l2003t iFiete et al 



200f: : ISchliemannl. 1200.'^: ISchliemann and MacDonak , 
200l IZa,rand and Ja.nkni \2()Q'i\ the theoretical studies 



have identified several mechanisms that can lead to non- 
collinear ground states. The observation that long wave- 
length spin-waves with negative energies frequently oc- 
cur within the parabolic-band kinetic- exchange model 
illustrates that randomness in the distribution of Mn 
moments can resul t in an instability of the coUinea r 
ferromagnetic state (|Schliemann and MacDonaldl 1200^ . 
Frustration can be further enhanced when positional dis- 
order is combined with anisotropics in Mn-Mn inter- 
actions. The pd character of electronic states form- 
ing the magnetic moment leads to magnetic interac- 
tion anisotropics with respect to the crystallographic 
orientation of the vector connecting two Mn moments 



l|Bre^^n^^^niez-Saiito3.|2003*.'Kiidrnovskv et al. 



■12004 

200,^ 



IMahadevan et all I2OO4 I Timm and MacDonald. 
When spin-orbit coupling is taken into account, magnetic 
interactions also become anisotropic with respect to the 
relative orientation of th e Mn-Mn connecti n g vector and 
the magnetic moment (| Fictc et'oT. 2 0051 ISchliemanii , 



I2OO3I: iTimm and MacDonaldl i2005t iZarand and Jankc , 



1. Magnetization of an isolated Mn(rf'+hole) complex 

We start the discussion by identifying the key physical 
considerations that influence the ground-state magneti- 
zation of (Ga,Mn)As ferromagnets by focusing first on a 
single Mn(d^-l-hole) complex and approximating the to- 
tal magnetization in the coUincar state by a simple sum 
of individual (identical) Mn((i^-|-hole) complex contribu- 
tions. This crude model is used only to qualitatively clar- 
ify the connection between p — d hybridization and anti- 
ferromagnetic kinetic-exchange coupling, the sign of the 
hole contribution to total moment per Mn, and the ex- 
pected mean-field contribution to magnetization per Mn 
from the Mn local moments and from the antiferromag- 
nctically coupled holes. We also explain in this section 
that quantum fluctuations around the mean-field ground 
state are generically present because of antiferromagnetic 
character of the p — d kinetic exchange interaction. 

Magnetization at T = is defined thcrmodynamically 
by the dependence of the ground-state energy E on ex- 
ternal magnetic field B: 



m 



dE 
'dB 



(24) 



B=0 



To avoid confusion that may result from using the hole 
picture to describe magnetization of carriers in p-type 
(Ga,Mn)As materials, we recall first the relationship be- 
tween magnetizations evaluated using the physically di- 
rect electron-picture and magnetizations evaluated using 
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the indirect but computationally more convenient hole- 
picture. In mean-field theory the magnetization is re- 
lated to the change of single-particle energy with field, 
summed over all occupied orbitals. Orbitals that de- 
crease in energy with field make a positive contribution 
to the magnetization. For B j| +z, the c?-electron spins 
arc aligned along (— z)-dircction (down-spins) and the 
majority spin band electrons have spin-up due to an- 
tiferromagnetic p — d exchange coupling. Then, if the 
majority band moves up in energy with B and the mi- 
nority band moves down, as illustrated in the left part 
of Fig.|^ the band kinetic energy increases with B and, 
according to Eq. (|^ . the corresponding contribution to 
the magnetization is negative. In the hole-picture, we 
obtain the same respective sense of shifts of the major- 
ity hole and minority hole bands, as illustrated in the 
right part of Fig. |2] and therefore the correct (nega- 
tive in our case) sign of the magnetization. The car- 
toon shows that in order to circumvent the potentially 
confusing notion of the spin of a hole in magnetization 
calculations, it is safer to start from the full Hamilto- 
nian Ti.{B) in the physically direct picture of electron 
states, where the sign of the coupling of electron spin 
to the field and the exchange energy are unambiguously 
defined. The electron picture —^ hole picture transfor- 
mation {Ti{B) —H{B)) and the clearly defined notion 
of majority and minority bands in either picture guaran- 
tees the sign consistency of the calculated magnetization. 
Note that the language used here neglects spin-orbit in- 
teractions which lead to single-particle orbitals that do 
not have definite spin character. Although spin-orbit in- 
teractions are important t hey can be neglected in most 
qualitative considerations l)Jungwirth et aZ.l . l2005a|) . 




FIG. 21 Cartoon of Zeeman coupling with an external mag- 
netic field assuming g > 0, in the electron and in the hole 
picture. In the ferromagnetic state the valence band is spin- 
split at zero magnetic field (solid lines) . The majority band in 
both electron and hole pictures move up in energy when the 
field is applied (dashed lines), resulting in a negative band- 
contribution to the magnetization. 

The electron-electron exchange energy has a nega- 
tive sign and its magnitude increases monotonically 



when moving from the paramagnetic to the half-metallic 
(empty minority band) state. This together with Eq. 124|) 
implies that the magnetization contribution from the 
electron-electron exchange energy has the same sign as 
the contribution from the kinetic energy. Using the same 
arguments as above we see that in the electron-electron 
exchange energy case the sign of magnetization is also 
treated consistently by the electron picture — > hole pic- 
ture transformation. 

The mean-field ground state wavefunction of the 
Mn(d^-l-hole) complex is = -S)\jz = +j) and the 
magnetization per Mn equals tumf = {gsS — gji)l^B, 
where S and j are local d-electrons and hole moments 
and gs and gj are the respective Lande g-factors. The 
five d-electrons have zero total orbital angular momen- 
tum, i.e. gs = 2, and for spin j = 1/2 hole {gj = 2) 
we get TTiMF = ^fJ-B- Hole states near the valence 
band edge have p-character, however, so more realisti- 
cally we should consider gjj = 4/3*3/2 = 2 which gives 
toa/ f = SfiB- We show below that this basic picture of a 
suppressed niMF due to holes applies also to highly Mn- 
dopcd (Ga,Mn)As materials, although the magnitude of 
the mean-field hole contribution is weaker because of the 
occupation of both majority and minority hole bands 
and, partly, because of spin-orbit coupling effects. 

The two-spin, S and j model can be also used 
to demonstrate the presence of quantum fluctuations 
around the mean-field ground state, which is a conse- 
quence of the antifcrromagnetic sign of the S • j coupling 
(jJungwirth et al, 2005a). In the limit of -B ^ the 
two-spin Hamiltonian is given by, 

^5,= JS.j = ^(^L-^'-?)- (25) 

For antifcrromagnetic coupling (J > 0), Stot = S — j and 
the corresponding ground-state energy Eaf = ■^[(S' — 
j){S - J + 1) - S{S + 1) - j{j + 1)] = ~\J\{Sj + j) is 
lower than the mean-field energy, —\J\Sj. The mean- 
field ground state is not exact here and quantum fluctu- 
ation corrections to the magnetization will be non-zero 
in general. The difference between magnetizations of the 
exact and mean-field state is obtained from respective ex- 
pectation values of the Zeema n Hamiltonian, q cj UfjBSz + 
gjfisBjz, and from Eq. ijJungwirth et adl2005a|) : 

j 

m - ruMF = mQp ^ -^b . (gs - .gj) ■ (26) 

When j = 1/2 and gs = gj ~ 2 the quantum fluctuation 
correction to the magnetization vanishes even though the 
mean- field ground state is not exact. The correction re- 
mains relatively weak also in the case of j = 3/2 and 
gj = 4/3), for which rngF ~ — 0.25yU_B. 

2. Magnetization of (Ga,IVIn)As ferromagnets 

As in the Mn(d^-l-hole) complex, the magnetization 
of coupled Mn moment systems can be decomposed 
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into mean-field contributions from Mn local moments 
and valence band holes and a quantum fluctuations cor- 
rection. At a mean-field level, the TBA description 
of (Ga,Mn)As mixed crystals is particularly useful for 
explaining the complementary role of local and itin- 
erant moments in this p-type magnetic semiconductor 
and we therefore start by revie wing results of this ap- 
proach (pmig wirth et aZ.( . l2005aj) . In Fig. [211 the micro- 
scopic TBA/ CPA magnetic moments per Mn, ttitba, 
in (Ga,Mn)As ferromagnets are plotted as a function of 
p/Nmu- The value of ttitba is obtained here using the 
electron picture by integrating over occupied states up 
to the Fermi energy. Spin-orbit coupling is neglected in 
these TBA calculations and only the spin-polarization 
contribution to magnetization is considered in niTBA, 
which simplifies the qualitative discussion below. 

A common way of microscopically separating contribu- 
tions from local atomic and itinerant moments is by pro- 
jecting the occupied electron states onto Mn d-orbitals 
and sji-orbitals, respectively. In this decomposition, the 
resulting local Mn moments are smaller than 5fiB per Mn 
due to the admixture of d-character in empty states near 
the valence band edge. The effective kinetic-exchange 
model corresponds, however, to a different decomposi- 
tion of contributions, in effect associating one spectral 
region with local Mn moments and a different spectral re- 
gion with itinerant hole moments. The kinetic- exchange 
model, in which local moments have 5* = 5/2, is obtained 
from the microscopic models, e.g. from the TBA/CPA, 
by expressing vtltba as the difference between a contribu- 
tion m™^ji^ calculated by integrating over all electronic 
states up to mid-gap, i.e. including the entire valence 
band, and a contribution corresponding to the integral 
from Fermi energy to mid-gap. As long as the valence- 
conduction band gap is non-zero, the former contribution 
is independent of valence band filling and equals the mo- 
ment of an isolated Mn atom, The latter term, 
plotted in the lower inset of Fig. [5U represents magneti- 
zation of itinerant holes. 

The applicability of the kinetic- exchange model relies 
implicitly on the perturbative character of the micro- 
scopic p — d hybridization. The level of p ~ d hybridiza- 
tion over a typical doping range is illustrated in Fig. 
which shows the orbital composition of m™^^. The filled 
symbols are calculated including spectral weights from all 
spd orbitals while the half-open and open symbols are ob- 
tained after projecting onto the d and sp orbitals, respec- 
tively. If there was no hybridization, m^^^ projected on 
the d-orbitals would equal the total Ti™^^ and the sp- 
orbital projected m^xBA would vanish. In the TBA/CPA 
calculations, the d-orbital projected m^xBA is reduced by 
only 10% as compared to the total m^rPgj^ and, there- 
fore, the p — d hybridization can be regarded as a weak 
perturbation. The nearly constant value of the d-orbital 
projected tu^tba ^-Iso suggests that the kinetic-exchange 
coupling parameter Jpd in the effective kinetic-exchange 
Hamiltonian is nearly independent of doping over the 
typical range of Mn and hole densities. 




FIG. 22 Main panel: mean-field total magnetization per Mn 
as a function of hole density relative to the local Mn moment 
density. Lower inset: hole contribution to magnetization (see 
definition in the text) per Mn. Upper inset: hole contribution 
to magnetization per hole. Resul ts are obtai ned using the 
TBA/CPA model. From LJunewirth et adl2005aiV 
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FIG. 23 Integrated total and d- and sp-projected magnetiza- 
tions per Mn as a function of hole density relative to the local 
Mn moment density. See text for definition of m. 
Jjunawirth et a/.. .200'5ah . 



'tba- From 



The decrease of tutba in Fig. 1221 with increasing 
p/Nmu clearly demonstrates the antiferromagnetic p — 
d exchange. The initial common slope for data cor- 
responding to different Mn concentrations reflects the 
half-metallic nature of the hole system (only majority 
hole band occupied) when spin-orbit interactions are ne- 
glected. Here the hole contribution to magnetization per 
volume is proportional to p, i.e., magnetization per Mn is 
proportional to p/N^n- The change in slope of ttitba at 
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larger hole densities, which now becomes Mn-density de- 
pendent, reflects the population of the minority-spin hole 
band and, therefore, the additional dependence of hole 
magnetization on exchange splitting between majority- 
and minority-hole bands. Note that the maximum abso- 
lute value of the hole contribution to magnetization per 
hole (see upper inset of Fig. I22|) observed in the half- 
metallic state is IfiB in these TBA calculations which 
assume j = 1/2 and gj = 2 holes. 

Similar conclusions concerning the character of contri- 
butions to the magnetization of (Ga,Mn)As have been 
inferr ed from LDA-I-U and SIC-LSDA supercell calcula - 
tions l)Schulthess et aZ.I . l2005HWierzbowska et allhOQd^ . 
(These microscopic calculations also neglect spin-orbit 
coupling.) The half- metallic LDA+U band structure 
in the case of zero charge compensation (jp/Nmu = 
1) results in a total ma gnetization per Mn of 4/is 
(|Wierzbowska et al\ . l2004l) , in agreement with the cor- 
responding rriTRA value s. In the SIC-LSDA calculations 
ijSchulthess et all l2005(l . the system is not completely 
half-metallic and, consistently, the total moment per Mn 
is larger than 4fiB- The LDA+U and SIC-LSDA local 
moments on Mn are and 4.5/iB, respectively, in 

good agreement again with the d-projected m™^/! '^^^ 
ues. In both ab initio calculations the oppositely aligned 
moment on the As sublattice extends over the entire su- 
percell, confirming the delocalized character of the holes 
and the antiferromagnetic sign of the p — d exchange. 

The KL kinet ic- exchange mode l calcu lations 
l|Dietl et all l2001ht lJungwirth et all l2005a|) have 
been used to refine, quantitatively, predictions for the 
total magnetization based on the above microscopic 
theories. In particular the number of minority holes at a 
given total hole density is underestimated in these TBA 
and ab initio approaches. This is caused in part by the 
quantitative value of the exchange spin-splitting of the 
valence band which, e.g., in the TEA/CPA calculations 
is a factor of 1.5-2 larger than value inferred from 
experiment. Also neglecting the spin-orbit interaction 
results in three majority bands that are degenerate at 
the F-point, instead of only two bands (heavy-hole and 
light-hole) in the more realistic spin-orbit-coupled band 
structure. (This deficiency is common to all calculations 
that neglect spin-orbit coupling.) In addition to having 
more states available in the majority band which leads 
to underestimating the minority hole density, these 
microscopic calculations also omit the reduction of the 
mean spin-density in the majority band caused by the 
spin-orbit coupling. The total magnetization values 
would be underestimated due to these effects. On the 
other hand, assuming only the spin contribution to 
the hole magnetization leads to an overestimated total 
magnetization, as already illustrated in Section lV.B.ll 

In the kinetic-exchange effective model the T = lo- 
cal moment contribution to the magnetization per Mn 
is 5iiB- As emphasized above, this is not in contra- 
diction with the smaller d-electron contribution to the 
magnetic moment in microscopic calculations. The ki- 



netic band energy contribution to the mean-field mag- 
netization per Mn is obtained by numerically integrat- 
ing over all occupied hole eigenstates of the Hamiltonian 
(|21|l and by finding the coefficient linear in B of this ki- 
netic energy contribution to the total energy ijPietl et all 
l2001bt lJungwirth et a/.ll2005a|) . Results of such calcula- 
tions are summarized in Fig. 1241 which shows the spin 
and orbital contributions to the magnetization of holes, 
and in Fig. [53 showing hole moment per Mn, m'^^^/p, 
for several local Mn moment and hole densities. Note 
that the decoupling of the hole magnetization into spin 
and orbital terms is partly ambiguous in the spin-orbit 
coupled valence bands and that only the total moment, 
TTT-'tVl', has a clear physi cal meaning l)Dietl et all l2001bt 
lJungwirth et aZ.l . l200 5al). As expected holes give a nega- 
tive contribution to magnetization, i. e. they suppress the 
total magnetic moment. The magnitude of the mean- 
field magnetization per hole, |m^}'J^| A^7\/„/p, is smaller 
than 2/iB obtained in Section lV.B.ll for the isolated spin- 
orbit coupled hole bound to the Mn impurity. It is due 
to occupation of both majority and minority heavy- and 
light-hole bands at these typical (Ga,Mn)As hole densi- 
ties (see inset of Fig. I25|l . Data shown in Fig. [21 and 
in the main panel of Fig. 1251 indicate a ~0.2 to 0.4/iB 
suppression of the mean-field moment per Mn due to the 
hole kinetic energy contribution to magnetization. 
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FIG. 24 Magnetization of tlie hole liquid (squares) in 
Gai_iMni,As computed as a function of the hole concentra- 
tion for the spin splitting parameter Bg = —30 meV ( corre- 
sponding to a; = 0.05 at T = 0). The crosses show spin and 
orbital contribution to the hole magnetization. The results 
were obtained using the six-band Kohn-Luttinger parameter- 
ization of_the_vaJ«n£e_te^ and the kinetic-exchange model. 
From (iDietl et all l2001bD . 

The hole exchange energy contribution to the total 
mean-field magnetization was found to be negative and 
nearly independent of x and p in the typical doping range, 
and its magnitude is about a factor of 5 smaller than the 
magnitude of the term originating fro m the hole kinetic 
band energy ijJungwirth et alj l2005a|) . Quantum fluctu- 
ation corrections lead to a 1% suppression of the mean- 
field moment per Mn (jjungwirth et a^j . l2005a|l . (Details 
of these calculations, using the imaginary time path- 
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FIG. 25 Mean-field kinetic energy contribution to the hole 
magnetization per Mn as a function of hole density. The 
results were obtained using the six-band Kohn-Luttinger pa- 
rameterizatio n of the valence band an d the kinetic-exchange 
model. From (),Tungwirth et aLLl20^5a^ . 



integral formulation of quantum many body theory com- 
bined with the Holstein-PrimakofF bosonic r epresentation 
of the Mn lo cal moments can be foun d in jKonig et al\ . 
l2001atl and l|.Tungwirth et oil l2005at) .) Combining all 
these considerations the T = magnetization per Mn in 
the effective kinetic-exchange model has a positive mean- 
field contribution equal to 5fj,B from the Mn local mo- 
ments and a negative contribution from band holes and 
quantum fluctuations which suppress the moment per Mn 
by -5-10%. 

The total magnetization per nominal Mn density and 
per effective density of uncompensated Mnca local mo- 
ments has been measured by a superconducting quantum 
interference device (SQUID) i n a series of as-grown an d 
annealed (Ga,Mn)As samples l).Iungwirth et aZ.! . l2005a^ . 
The characterization of these materials has already been 
discussed in the previous section (see Fig. 1201 and the 
related text). Within experimental uncertainty, the 
SQUID magnetization was found to be independent of 
the magnetization o rientation, in ag:reement with theo- 
retical expectations (|Jungwirth et all l2005a|) . The mo- 
ment decreases with increasing nom inal Mn concen- 
tration, and increases on an nealing l|,Iungwirth et al\ . 
l2005al: IPotashnik et oil lloHS)- This is consistent with 
the anticipated f ormation of inter s titial M n for dop- 
ing above ~2% (|Jungwirth et all l2005b(l . given the 
antiferromagnetic coupl i ng b e tween Mnj an d Mng a 
l)Blinowski and Kacmanl l2003t lEdmonds efdl . l2005a() . 
and with breaking of this coupling by low-temper ature 
annealing ^Edmonds et all l2nn4aUYu et d\. 12001. In 
agreement with the above theoretical calculations, the 
total magnetization per effective density of uncompen- 
sated MuGa, "^sgfy/Lij falls within the range 4-5^b for 
all samples studied. Furthermore, although there is ap- 
preciable scatter, it can be seen that samples with lower 



hole densities tend to show higher m'ggjjjjj, consistent 
with a negative contribution to magnetization from an- 
tiferromagnetically coupled band holes. 

The local and d-statc projected contribution from Mn 
to the magnetic moment in (Ga,Mn)As has been probed 
experimentally by meas uring the x-ray rnagnetic circu- 
lar dichroism (XMCD) llEdmonds et oil l2004b[ l2005at 
lJungwirth et allhoOba^ . In agreement with the SQUID 
measurements and theoretical expectations, the XMCD 
data are independent, within experim ental uncertainty , 
of the direction of magnetization (jJungwirth et alj 
The data are listed in Table ITU for two annealed 
samples with low and high Mn doping. In both cases, 
magnetic moments of around 4.5/zb were obtained, show- 
ing a negligible dependence on the hole density. Sim- 
ilar results were found for samples with intermediate 
Mn doping l).Tungwirth et all l2005a(l . The experimen- 
tal XMCD results are in good agreement with the corre- 
sponding TBA values indicated by half-open symbols in 
Fig. 123 an d with the LDA+U and SIC-LSDA Mn local 
moments (|Schulthess et oil l200f)l: IWierzbowska et all . 
I2OO4I) . (Note that all these microscopic calculations ac- 
count only for the spin angular momentum contribution 
to the local Mn 3c? moment since spin-orbit coupling ef- 
fects were neglected.) 
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FIG. 26 SQUID magnetization per efi'ective density of un- 
compensated MuGa local moments in as-grown (open sym- 
bols) and annealed (filled symbols) (Ga,Mn) As materials plot- 
ted as a function of hole density per efi'ecti ve density of un- 
compe nsated MuGa local moments. From (|,Tungwirth et ail 
l2005al) . 



C. Micromagnetic parameters 

A small set of parameters is often sufficient to phe- 
nomenologically describe the long-wavelength properties 
of ferromagnets. The description, which usually captures 
all properties that are relevant for applications of mag- 
netic materials, starts from the micromagnetic energy 
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spin rn°'^'> rn'^P^" ^ 

■'' '^XMCD ^XAICD ^XMCD'^^XMCD 



(%) (±0.3^5 ) (±0.03^5 ) 



(±0.3mb) 



2.2 



4.3 



0.15 



4.5 



8.4 



4.3 



0.16 



4.5 



TABLE II Mn 3d moments obtained from XMCD and de- 
composed into the spin and orbital contributions in annealed 
samples with nominal Mn dopings 2.2 and 8.4%. From 
(|,Tungwirth et aLLl2005ah . 



functional e[n] of the s patially depende nt magnetization 
orientation n = M/M ()AliaroniL l200l|) . The micromae- 
netic energy functional treats the long-ranged magnetic 
dipole interactions explicitly and uses a gradient expan- 
sion for other terms in the energy. The magnitude M of 
the magnetization is one of the micromagnetic parame- 
ters characterizing the material. Zeroth order terms in 
the energy functional arc the magnetic anisotropy energy, 
Gam Mi ^^^^ thc Zccman coupling to an external mag- 
netic field if present, — /xqH • nM. The leading gradient 
term, e^x [n] , referred to as exchange energy in micromag- 
netic theory, represents the reduction in magnetic con- 
densation energy when the magnetization orientation is 
not spatially constant. Micromagnetic parameters used 
to characterize the magnetic anisotropy energy depend 
on the symmetry of the system. For example, in a fer- 
romagnet which possesses uniaxial anisotropy with the 
easy-axis aligned along the 2;-direction, Cani = Kuul, 
where < is the uniaxial anisotropy constant. As we 
explain below, thc magnetic anisotropy of (Ga,Mn)As 
fcrromagnets is a combination of the cubic term, the 
in-planc uniaxial anisotropy, and the uniaxial term in- 
duced by thc growth-direction lattice-matching strain 
which often dominates in (Ga,Mn)As epilayers. The large 
anisotropics and relatively small magnetic moments of 
these dilute magnetic systems rank (Ga,Mn)As DMSs 
among hard fcrromagnets (magnetic hardness parameter 
K ~ |i<r„//ioM^|^/^ > 1) with outstanding micromag- 
netic properties, including frequently observed single- 
domain-like characteri stics of field-induced magnetiza- 
tion reversals (.Abolfath adl2001btlpietl et ad. 120011: ; 



Goennenwein et all l2005t lOhnd. Il998l: iPotashnik et al. 
2003HWang et ffl^.l . l200fid^ . 

The anisotropy of the exchange term in thc micromag- 
netic functional is often neglected in which case it can be 
written as, Cex = A(Vn)^, where A is the spin stiffness 
constant. Collective magnetization dynamics is described 
by the Landau-Lifshitz-Gilbert equation. 



dM 
~dt 



where /ioHg = — 9e/(9M and a is the Gilbert damp- 
ing micromagnetic parameter. In this section we review 
KL kinetic-exchange model calculations of the micromag- 
netic parameters of m e tallic (Ga,Mn)As fcrromagnets 
(lAbolfath et all l2001at iBrev and Gomez-Santol 1200.1 " 



. - . - -^^ G6mez-Sa,nto4 l2003| 

Dietl et all l2001albl : iKonig et all l2001at ISinova et all 

and predictions based on the microscopic values 



of these parameters for anisotropy fields, the character- 
istic size of domains, and the crit ical current in current - 
induced magnetization switching JAbolfath et all l2001at 
iDietl et g^J . l2001alal: IShlova et aZ.Ll2004b|) . 



1. Magnetocrystalline anisotropy 

Magnetocrystalline anisotropy, which is the depen- 
dence of the energy of a ferromagnet on the magneti- 
zation orientation with respect to crystallographic axes, 
is a spin-orbit coupling effect often associated with lo- 
calized electrons in magnetic d or /-shells. Local Mn 
moments in (Ga,Mn)As, however, are treated in the KL 
kinetic exchange model as pure spins S = 5/2 with 
angular momentum L = and, therefore do not con- 
tribute to thc anisotropy. The physical origin of the 
anisotropy energy i n this model i s spin-orbit coupling in 
thc valence band (|Abolfath et al . ,2001a.: Dlctl , 2001, ). 
Even within the mean-field approximation to the KL 
kinetic-exchange model, magnetic anisotropy has a rich 
phenomenology which has explained a number of ex- 
perimental observations in the (III,Mn)V DMSs, includ- 
ing easy axis reorientations as a function of hole den- 
sity, temperature, or strains i n the lattice l^ iu 



2nn4bll2nn3lMa,sma,nidis_^^ gilliinilOhno et aiillflflB^ 
Sawicki et oil 12004. .200513) . 



The remarkable tunability of thc magnetic properties 
of (Ga,Mn)As DMSs through lattice matching strains 
is an important byproduct of LT-MBE growth of fer- 
romagnetic films with lattices locked to those of their 
substrates in the plane perpendicular to the growth axis. 
X-ray diffraction studies have established that thc re- 
sulting strains are not relaxed by dis locations or other 
defects even for ~ 1/xm thick epilayers l|Shen et a^J . ll99ii 
IZhaoIiriDilioni)- We mentioned in Section IrTbI that 
the lattice constant of relaxed (Ga,Mn)As is larger than 
the lattice constant of GaAs, especially so if interstitial 
Mni or AsGa antisites are present in the DMS crystal. 
(Ga,Mn)As grown on the GaAs substrate is therefore 
under compressive strain. Tensile strained (Ga,Mn)As 
DMSs have beeii produced by using (Ga,In)As substrate s 
•^.120031 lOhno et giJ.ll 996aUShono et alVM)^. 

[001] growth direction strain breaks the cubic symme- 
try of the (Ga,Mn)As resulting in a combined cubic and 
uniaxial anisotropy form of the energy functional, 



TV' ( 2 2, 2 2, 2 2\ 

= Kci {n^Uy + n^n^ + nyU^} 
+ Kc2 {nxnyn^f + Kuul . 



(28) 



Here Kd and Kc2 are the two lowest order cubic 
anisot ropy constants (|Abolfath et giJ . l2001atrDietl et all . 
l2001bh . Strain must be included in the k-p description 
in order to evaluate the uniaxial anisotropy constant K-^. 
For small strains this is done by expressing the posi- 
tional vector r' in the strained lattice in terms of r in 
the unstrained lattice as = ra + X]/3 ^a/s?"/?, and ex- 
panding thc KL Hamiltonian in lowest order of the strain 
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const ants ea0 ijChow and KochLll999l:lJones and MarchL 
Il973j) . In (Ga,Mn)As epilayers grown along the [001] di- 
rection, the strain constant Cxx = Byy can be tuned from 
approximately —1% to +1%. For the larger strain values, 
the uniaxial term dominates the total ani sotropy energy 
llAholfath fi^ r7i J. l2001aUDiet1 p.t aimomH) . InFig.EHwe 
show mean-field KL kinetic-exchange model calculations 
of the uniaxial anisotropy field fiQHu — \2Ku/M\ rela- 
tive to the T = magnetization M = g^sNMnS ■ ^j-qHu 
corresponds to the minimum external magnetic field nec- 
essary to align magnetization M along the hard axis. The 
figure illustrates the dependence of the easy axis orienta- 
tion on the hole density and on the sign of Cxx- In par- 
ticular, the easy-axis is in-plane for compressive strain 
{cxx < 0) and out-of-plane for tensi le strai n (e^x > 0) , 
consis tent with experiment (jLiu et a/. . .2003:.Qhno et al\ . 
Il996a|) . The ability to manipulate the easy axis orienta- 
tion from in-plane to out-of-plane has many implications 
for fundamental research on DMS materials and is very 
attractive also from the point of view of potential appli- 
cations in the magnetic recording technologies. 

(Ga,Mn)As epilayers have also a rel atively strong in- 
plane uniaxia l anisotropy c omponent JLiu et all l2003t 
ISawicki et oil 12004 l2005bl: iTang et a^J^2003^ . The in- 
plane easy axis is consistently associated with particular 
crystallographic directions and can be rotated from the 
[TlO] direction to the [110] direction by low-temperature 
annealing. Although the origin of the in-plane uniax- 
ial anisotropy has not been established, its dependence 
on hole concentration (varied by annealing) and temper- 
ature was modeled successfully within the KL kinetic- 
exchange model by assuming th at it is associated wit h a 
small shear strain Cxy ~ 0.05% (|Sawicki et a/.L l2005b|) . 

2. Spin stiffness 

Thermal and quantum magnetization fiuctuation ef- 
fects in metallic (Ga,Mn)As DMSs have been described 
by Holstein-Primakoff boson representation of the Mn 
local spin operators. Assuming that fiuctuations around 
the mean-field orientation are small, the relationship be- 
tween spin-raising and lowering operators and boson cre- 
ation and annihilation operators is, 5+ = b\/2NMnS and 
S- = b U/2NMr,.S iS^. = {S+ + S- ) I2. k, = {S+ - 
S~)/2i) (lAuerbachL ll9M iKonig etdl 1200 lai l2000ft . Af- 
ter integrating out the itinerant hole degrees of freedom 
in the coherent-state path-integral formalism of many- 
body problem, the partition function, 

Z = Jv[zz]cxp{-Seff[zz]), (29) 

depends only on the bosonic degrees of freedom (repre- 
sented by the complex numbers z and z in the effective 
action Seff)- The independent spin- wave theory is ob- 
tained by expanding Seff[zz] up to quadratic order in z 
and z, i.e., spin excitations are treated as noninteracting 
bosons. The spin stiffness parameter A is then calcu- 
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FIG. 27 Computed minimum magnetic field Hun (divided 
by M) necessary to align magnetization M along the hard 
axis for compressive (a) and tensile (b) biaxial strain in 
Gai_a:Mn2,As film for various values of the spin-splitting pa- 
rameter Bg- The easy axis is along [001] direction and in 
the (001) plane at low and high hole concentrations for com- 
pressive strain, respectively (a). The opposite behavior is ob- 
served for tensile strain (b). The symbol [100] — > [001] means 
that the easy axis is along [100], so that Hun is applied along 
[001] {Ba = —30 meV correspo nds to the saturatio n value of 
M for Gao.g5Mno.05As). From iPietl et aLLl2001bD 



lated by fitting the microscopic spin-wave dispersion at 
long wavelength to the form. 



^^ = ^ + ^^' + ^(^') (30) 

to match the conventions used for exchange and 
anisotropy constants in micromagnctic theory. Typical 
values of A in (Ga,Mn)As derived from the above many- 
body formalism and the KL kinetic-exchange description 
of the hole bands are shown in Fig. The values are 
consistent with the experimental spin stiffness parameter 
in a G an.949Mnn.n.qi As measured by ferromagnetic reso- 
nance ijGoennenwein et a/.l . 12003^1 . 
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FIG. 28 Theoretical spin stiffness parameter in 
Gao.95Mno.05As as a function of the hole density cal- 
culated using the Holstein-Primakoff representation of 
fluctuating Mn local spins and the KL kinetic-exchange 
description of the hole bands. 

3. Gilbert damping of magnetization precession 

The damping of small cone-angle magnetization pre- 
cession in a ferromagnet is parameterized by the Gilbert 
coefBcient. For small fluctuations of the Mn magneti- 
zation orientation in (Ga,Mn)As around the easy axis, 
Eq. (|27() can be used to derive an expression for the linear 
response of a magnetic system to weak transverse fields 
in terms of the phcnomcnological constants of micromag- 
netic theory. For zero external static magnetic field and 
zero wavevector (uniform rotation) the corresponding in- 
verse susceptibility reads: 

^ fi f Ku - iat^ -iuj \ 

{g^-LBYNMnS V «w Ku - iau} J 

where Ku = Ku/ (hNMnS) and uj is the frequency of the 
external rf field perturbation. 

Microscopically, Gilbert damping in (Ga,Mn)As DMSs 
was attributed to the p ~ d exchange-cou pling between 
local Mn moments and itinerant holes ijSinova et alV 
l2004b() . The elementary process for this damping mech- 
anism is one in which a local-moment magnon is anni- 
hilated by exchange interaction with a band hole that 
suffers a spin-flip. This process cannot by itself change 
the total magnetic moment since the exchange Hamil- 
tonian commutes with the total spin S -I- s. Net relax- 
ation of the magnetization requires another independent 
process in which the itinerant hole spin relaxes through 
spin-orbit interactions. A fully microscopic theory of the 
kinetic-exchange contribution to the Gilbert coefficient 
was derived by comparing Eq. 131() with microscopic lin- 
ear response theory and by identifying the Gilbert coeffi- 
cient with the dissipative par t of the quantum-mechanica l 
susceptibility ijSinova et l2004bt iTserkovnvak et aH 



120041) . 

X5(r,i|r',t') = {gtiBf'-{[S,{r,t),S,{v\t')])e{t-t'). 

(32) 

Here i = x,y and Si{r,t) = Mi[r,t)/{g^B) are the Mn 
transverse spin-density operators. As in the microscopic 
theory of the spin-stiffness, the correlation function H32|) 
for the uniform (k = 0) precession mode was evaluated 
using the the long-wavelength non -interacting spin- wav e 
form of the partition function H29|l ijSinova et a/J . l2004bfl . 
Gilbert coefficients in (Ga,Mn)As DMSs obtained within 
this formalism agree quantitatively with experimental 
values of a in homogeneous (annealed) systems measured 
from the width of the ferromagnetic resonance curves, as 
shown in Figs. 1291 and I3UI 




FIG. 29 Theore tical Gilber t dam ping coefficient in 
(Ga,Mn)As. From jSinova ef aZ.ll2004bD . 



4. Domains and spin-transfer magnetization switching 

Theoretical values of magnetic anisotropy and spin 
stiffness have been combined to estimate the typical do- 
main size in tensile str ained (Ga,Mn)As e pilayers with 
out-of-plane easy axis ijPietl et alV l2001at) . The calcu- 
lated low-temperature width of a single domain stripe 
of 1.1/im compares favorably with the experimental 
value of 1.5/xm seen in the micro-Hall probe experiments 
ijShono et al\ . l2000|) . Near the discrepancy between 
theoretical and experimental domain sizes becomes large, 
however, and has been attributed to critical fluctuations 
effects not included in the mean-field theory. 

The calculated Gilbert damping coefficient and mag- 
netic anisotropy constants were used to predict crit- 
ical current s for s pin-transfer magne tization switch- 
ing l)Bergeil Il996t Isionczewskil Il996l) in (Ga ,Mn)As 
based tunneling structures ijSinova et ai\ . \2QM\i) . Spin- 
polarized pcrpcndicular-to-planc currents in magnetic 
multilayers with non-parallel spin configurations can 
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FIG. 30 Experimental peak-to-peak ferromagnetic resonance 
linewidth and corresponding Gilbert coefficient in as-grown 
(filled symbols) and annealed (open symbols) Gao.92Mno.08As 
samples measured as a function of temperature for [001] and 
[110] dc magnetic-field orientations (main pl ot) and as a func- 
tion of the field ang le at 4K (inset). From ih'm et all 12003 " 
ISinova et all 120041 



transfer spin betwe en magnetic l ayers and exert current- 
dependent torques (|Slonczewski Il99^ . The critical cur- 
rent for magnetization switching is obtained by adding 
the torque term to the Landau-Lif s hitz-G ilbert Eq. H27(l 
l)Sinova all . l2004bl: ISlonczewski Il996(l . Critical cur- 
rents ~ 10^ A cm~^ obt ained in these calcu lations and 
confirmed in experiment (IChiba et a^J ■ l2004b^ are two or- 
ders of magnitude smaller than those observed typically 
in metals. The small moment densities explain a large 
part of the orders of magnitude reduction in the criti- 
cal current. This finding suggests that DMS materials 
have the potential to be particularly useful for exploit- 
ing the current-induced magnetization reversal effect in 
magnetic tunnel junctions. 



VI. MAGNETO-TRANSPORT 

Studies of the temperature dependent resistivity, 
anisotropic magnetoresistance, and anomalous and ordi- 
nary Hall effects have been used to characterize DMS 
materials and to test different theoretical models de- 
scribing these ferromagnets. In this section we review 
dc magneto-transport properties of (Ga,Mn)As focusing 
mainly on the metallic regime. 



T ^ 0, although for practical reasons this adjective is of- 
ten used to describe a material whose resistivity decreases 
with temperature over most or all of the range of tem- 
peratures studied in a particular series of experiments.) 
In optimally annealed samples with low density of un- 
intentional defects, metallic behavior is observed for Mn 
doping larger than approximately 1.5% ijCampion et all 
l2003at .Potashnik et oL, 2002). In Section HLDI we al- 
ready discussed this observation as a consequence of the 
Mott metal-insulator transition due to doping with sub- 
stitutional MuGa acceptors. We also introduced in Sec- 
tions III . D1 and lTlLDl some of the theoretical work that has 
qualitatively addressed magnetic and transport proper- 
tics of DMS systems near the metal-insulator transition. 
To our knowledge, a systematic experimental analysis has 
not yet been performed which would allow a reliable as- 
sessment of the theory predictions in this complex and 
intriguing regime. On the other hand, the dc transport 
in metallic (Ga,Mn)As DMSs has a rich phenomenology 
which has been explored in a number of experimental 
works and microscopic understanding of many of these 
effects is now well established. 

LSDA/CPA band structure calculations combined 
with Kubo linear response theory were used to study cor- 
relations between the low temperature conductivity and 
the density of various defects in the l attice, the hole den - 
sity, and Tc in metalhc (Ga,Mn)As (|Turek et all l^QQ^). 
The theory tends to overestimate the conductivity in low- 
compensation materials but the overall range of values 
~ 100 — 1000 il^^cm^^ for typical material parameters is 
consistent with measured data. As illustrated in Fig. ISTl 
and I32L the theory models capture at a qualitative level 
the correlation between large conducti vities and high 
Curie temperatures seen in experiment fC^^gionei^Zj, 
l2003at Ed mond s et ai, 2002a; Potashni k et all l2002(l . 
Comparable absolute values of the T = conductance 
and similar trends in the dependence of the conductance 
on the density of impurities were obtained using the KL 
kinetic-exchange model and the semiclassical Boltzmann 
descr i ption of the dc transport llHwang and Das Sarma , 
2005 ; lJungwirth et all l2002at iLopez-Sancho and Brev . 
2003). We note, however, that there are differences in 



the detailed microscopic mechanisms limiting the con- 
ductivity in the two theoretical approaches. Scattering 
in the ab initio theory is dominated by the p — d ex- 
change potential on randomly distributed Mn atoms, and 
by the local changes of the crystal potential on the impu- 
rity sites. Long-range Coulomb potentials produced by 
MuGa acceptors and other charged defects are omitted in 
the CPA approach. 



A. Low temperature conductivity 

(Ga,Mn)As materials can exhibit insulating or metallic 
behavior depending on the doping and post-growth an- 
nealing procedures. (Strictly speaking, a material is de- 
fined as being metallic if its resistivity is finite in the limit 



In the kinetic-exchange effective Hamiltonian model 
the elastic scattering effects were included using th e 
first order Born approximation ijjungwirth et al\ . \2002s^ . 
The corresponding transport weighted scattering rate 
from the MuGa impurity has a contribution from both 
the p—d exchange potential and the long-range Coulomb 
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FIG. 31 Mean-field LSDA/CPA Curie temperatures versus 
T = do conductivities calculated for (Gao.gssMno.osAsj; 
)As alloys with varying As antisite content y (full dots) and 
the experimental values obtained in as-grown and annealed 
(Ga,Mn)A s thin films with the c orresponding level of com- 
pensation iCampion et aLll2003ar) (open squares). Note that 
these are illustrative calculations since in the experiment the 
change in the unintentional impurity concentration upon an- 
nealin g is mostly due to the out-diffusion of interstitial Mn. 
From l|Turek et all l2004l . 




200 400 600 

conductivity (CI cm ) 

FIG. 32 Curie temperature versus room temperature con- 
ductivity for (Ga,Mn)As films with x=0.08 (squares), x=0.06 
(circles), x=0.05 (triangles). Straight lines are to guide the 
eye. The data points correspond to samples annealed at dif- 
ferent temperatures and times. Note that similar trends are 
obtained for Curie temperatures plotted versus low tempera- 
ture conductivities, whose values are systematically larger by 
~ 20— 40% compared to the ro om temperature conductivities. 
From jEdmonds et all l2002ar) . 
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Here e^ost is the host semiconductor dielectric con- 
stant, \z^^) is the multi-component eigenspinor of the KL 
Hamiltonian H21() and j is the corresponding eigenen- 
ergy, and the Thomas- Fermi screening wavevector qtf = 
v/e2DOS(£;F)/(e/iosteo), where DOS(£;f) is the densit y 
of states at the Fermi energy (jjungwirth et g^j . l2002a|) . 
(Analogous expressions apply to scattering rates due to 
other defects.) 

The relative strengths of scattering off the p — d 
exchange and Coulomb potentials can be estimated 
by assuming a simple parabolic-band model of the 
valence band characterized by the heavy-hole effec- 
tive mass m* = 0.5me. The p — d kinetic- 
exchange contribution in this approximation is Tpd = 
^A/n./pd<S'^TO*V2m*i?F/(47r7i^) and the scattering rate 
due to the screened Coulomb potenti al, Fp, is given by 
the standard Brooks-Herring formula ijBrookj. Il955|) . In 
(Ga,Mn)As with p = 0.4 nm""^ and x = 5%, these es- 
timates give hTpii ^ 20 meV and TiTc ~ 150 meV. For 
lower compensation (p « 1 nm^'^), the screening of the 
Coulomb potential is more efficient resulting in values of 
Tc below 100 meV, but still several times larger than 
Tpd- (Note that these elastic scattering rates are smaller, 
although by less than a factor of 10, than other charac- 
teristic energy scales such as the Fermi energy and the 
spin-orbit coupling strength in the GaAs valence band 
which partly establishes the consistency of this theoreti- 
cal approach.) 

The dominance of the Coulomb potential in Born ap- 
proximation scattering rates for typical chemical compo- 
sitions is confirme d by calculations based on the six-band 
KL Hamiltonian l)Jungwirth et all l2002a(l . The good 
agreement between T = conductivity values obtained 
using the ab initio and the kinetic-exchange model the- 
ories should therefore be taken with caution as it may 
originate, to some extent, from the stronger local p — d 
exchange in the LSDA/CPA theory which partly com- 
pensates the neglect of long-range Coulomb potentials in 
this ab initio approach. 



B. Anisotropic magnetoresistance 

Boltzmann transport theory combined with the KL 
kinetic-exchange model of the (Ga,Mn)As band struc- 
ture is a practical approach for studying magneto- 
transport effects that originate from the spin-orbit cou- 
pling. In Section IV. CI we reviewed the predictions 
of the model for the role of spin-orbit interaction in 
magnetic properties of (Ga,Mn)As. Here we focus on 
the anisotropic magnetoresistance (AMR) effect which 
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The AMR effect can be regarded as the first spin- 
tronic functionahty implemented in microelectronic de- 
vices. AMR magnetic sensors replaced simple horse-shoe 
magnets in hard-drive read heads in the early 1990's. 
With the introduction of giant magnetoresistance based 
devices in 1997, magnetoresistive sensing launched a new 
era era in the magnetic memory industry. In ferromag- 
netic metals, AMR has been known for well over a cen- 
tury. However, the role of the various mechanisms held 
responsible for the effect has not been fully sorted out 
despit e renewed interest motivated by pract ical applica- 
tions l)Jaoul et all Il977t iMalozemofl Il985j) . The difh- 
culty in metal ferromagnets partly stems from the rela- 
tively weak spin-orbit coupling compared to other rel- 
evant energy scales and complex band structure. In 
(Ga,Mn)As ferromagnets with strongly spin-orbit cou- 
pled holes occupying the top of the valence band, mod- 
eling of the AMR effect can be accomplished on a semi- 
quantitative level without using free parameters, as we 
review below. 

(Ga,Mn)As epilaycrs with broken cubic symmetry due 
to growth-direction lattice-matching strains are usually 
characterized by two AMR coefficients. 



strain 



AMRop = 



AMRi 



p^^(M II x)- p^^jm II z) 
pxx[y^ II z) 

p^^(M 11 x) -p^^(M II y) 
P..(M II y) 



(35) 



where z is the growth direction and the electrical current 
/ II X. The three different experimental configurations 
used to determine AMRop and AMRip are illustrated 
in Fig. The coefficient AMRip is given by combined 
effects of spin-orbit coupling and the current induced bro- 
ken symmetry between the two in-plane cubic axes. The 
difference between AMRip and AMR^p is a consequence 
of the lowered symmetry in strained samples. Theoretical 
predictions for AMRip and AMRop were compared with 
measurements in a series of (Ga,Mn)As epilayers grown 
under compressive strain ijjungwirth et aZ.l . l2003bj) . The 
calculated sign of the AMR coefficients and the magni- 
tude which varies form ^ 1% in highly Mn doped mate- 
rials to ^ 10% in samples with low Mn concentration are 
consistent with experiment. 

Reminiscent of the magnetocrystalline anisotropy be- 
havior, the theory predicts rotations as a function of 
strain of the magnetization direction corresponding to 
the high (or low) resistance state. Calculations illustrat- 
ing this effect arc shown in Fig. 1341 and the experimen- 
tal demonstration in (Ga,Mn)As epilayers with compres- 
sive and tensile strains is presented in Fig. 1351 In the 
top panel of Fig. corresponding to (Ga,Mn)As ma- 
terial with compressive strain, M _L plane (M || z) is 
the high resistance configuration, the intermediate re- 
sistance state is realized for in-plane M _L I (M || y), 
and the low resistance state is measured when M || I 
(M II x). In the sample with tensile strain, the in-plane 
M _L I and M _L plane curves switch places, as seen in the 




FIG. 33 Schematic illustration of AMR configurations: elec- 
trical current I flows along i-direction and magnetization is 
aligned by an external magnetic field in either of the two in 
plane directions, M || I (M || x) and in-plane M _L I (M || y), 
or in the out-of-plane direction perpendicular to the current, 
M _L plane (M || z). The non-zero coeflicient AMRip, de- 
flned as the relative diflFerence in resistivities in the in-plane 
M ± I and M || I configurations, originates from spin-orbit 
coupling. In strained (Ga,Mn)As materials, AMRip is in gen- 
eral not equal to AMRop, defined as the relative difference in 
resistivities in the M J. p lane and M || I configurations. From 
llMatsukura et a/.ll200ll . 



bottom panel of Fig. (JSl Consistent with these experi- 
mental observations the theoretical AMRip and AMRop 
coefficients are negative and \AMRip\ < \AAIRop\ for 



compressive strain and \AMRipi 
strain. 



> \AMRop\ for tensile 



C. Anomalous and ordinary Hall effects 

The anomalous Hall effect (AHE) is another trans- 
port phenomenon originating from spin-orbit coupling 
which has been used to study and characterize fer- 
romagnetic films for more than one hundred years 
ijChien and Westeatd. Il980(l . The difficulties that have 
accompanied attempts at accurate microscopic descrip- 
tion of the effect in metal ferromagnets are reminiscent 
of those in the AMR. The success of AHE modeling in 
(Ga,Mn)As materials llDietl et aIl.l2003tlEdmonds et d\ . 



l2003t lJungwirth et all l2002cl l2003b() . reviewed in this 
section, has had implications also beyond the field of 
DMSs. It helped to motivate a reexamination of the AHE 
in transition metals and and in a scries of more complex 
ferromagnetic compounds which has led to a significant 
progress in resolving the microscopi c mechanism respon- 
sible for the effect in these i 
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The Hall resistance Rxy = Pxy/d of a magnetic 
thin film of thickness d is empirically observed to con- 
tain two distinct con t ributi ons, Rxy = RqB + RaM 
dCh ien and Westgatd. Il980(l . The first contribution 
arises from the ordinary Hall effect (OHE) which is pro- 
portional to the applied magnetic field B, and the sec- 
ond term is the AHE which may remain finite at S = 
and depends instead on the magnetization. Rq and Ra 
arc the OHE and AHE coefficients respectively. To put 
the AHE studies in DMSs in a broader perspective, the 
following paragraphs offer a brief excursion through the 
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FIG. 34 Theoretical AMR efltects in (Ga,Mn)As obtained 
from the KL kinetic-exchange model and Boltzmann trans- 
port theory. The curves correspond to total Mn doping x\ 
and compensation due to As-antisites or total Mn doping x-i 
and compensation due to Mn-interstitials. Main panel: out- 
of-plane AMR coefficient as a function of strain for several 
Mn dopings. Lower inset: out-of-plane and in-plane AMR 
coefficients as a function of strain. Upper inset: out-of-plane 
AMR coefficient as a f unction of the hole density. From 
iJungwirth et aLll2002ah . 



history of AHE theory and an overview of microscopic 
mechanisms discussed in this context in the literature 
(a mo re detailed survey is given, e.g., in ijSinova et al\ . 

liooii)). 

The first detailed theoretical analysis of the AHE was 
given by Karplus and Luttingcr ijKarplus and Luttingeil 
Il954j) , where they considered the problem from a pertur- 
bative point of view (with respect to an applied electric 
field) and obtained a contribution to the anomalous Hall 
conductivity in systems with spin-orbit coupled Bloch 
states given by the expression: 



2p^ I t 



.2 



n.k 



du 



(36) 



where /^^ g is the Fermi occupation number of the Bloch 
state |u„j:). This contribution is purely a property 
of the perfect crystal and has become known in re- 
cent years as the intrinsic AHE. It leads to an AHE 
coefficient Ra proportional to p^^ and therefore can 
dominate in metallic fcrromagnets that have a rel- 
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FIG. 35 Upper and middle panel: experimental field-induced 
changes in the resistance of Gao.95Mno.05As grown on GaAs 
substrate under compressive strain for current along the [110] 
crystal direction (upper panel) and [100] crystal direction 
(middle panel) for three different orientations of the mag- 
netization as indicated in Fig. 1331 Lower panel: experimen- 
tal magnetoresistance curves in Gao.957Mno.o43As grown on 
(In,Ga)As substrate under te nsile strain (current alo ng the 
[110] crystal direction). From jMatsukura et aZ.L [20041) . 



atively large resistivity. The intrinsic AHE is re- 
lated to Bloch state Berry phases in momentum space 
and depends non-perturbatively on spin-orbit interac- 
tion strength when degenerac ies in momentum space are 
lifted by spin- o rbit c o upling llBurkov and Balentsl 1200^ ; 



iDugaev et all l2005t iHaldand. |2004 lJungwirth et al 
'2002d[ lOnoda and Nagaosal l2002t ISundaram and Niu . 
il999j) . This point is particularly relevant for ferromag- 
netic semiconductors because all carriers that contribute 
to transport are located near particular points in the Bril- 
louin zone, often high symmetry points at which degen- 
eracies occur. 

Shortly after the seminal work of Karplus and Lut- 
tingcr, Smit proposed a different interpretation of 
the AHE based on a picture of asymmetric spin- 
dependent skew scattering off i mpurit y potential in- 
volving spin-orbit coupling llL eroux-HAjEjo^^ GhazaliL 
Il972t iNozieres and Lewineil Il973l ISmiti . Il95fil) . Analyt- 
ically, the skew scattering appears in the second order 
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Born approximation applied to the collision term of the 
Boltzmann transport equation. This mechanism gives a 
contribution to Ra oc pxx, proportional to the den- 
sity of scatterers and dependent on the type and range 
of the impurity potential. 

The AHE conductivity has a number of contributions 
in addition to the skew scattering and intrinsic contribu- 
tions, that can originate either from spin-orbit coupling 
in the disorder scattering or spin-orbit coupling in the 
Bloch bands. Among these, side jump scatte ring has 
been identified has an important contribution (jBergeil 
Side jump due to spin-orbit coupling in the Bloch 
bands appears as a ladder diagram vertex correction to 
the intrinsic anomalous Hall effect and its importance 
depends on the nature of that coupling ijDueaev et all 
l2nn,^ISinitsvn p.t r7i.l . l2f)0l . 

The ratio of intrinsic and skew-scattering contributions 
to the AHE conductivity can be approximated, assum- 
ing a single spin-orbit coupled band and scatter i ng off 
ionized impurities, by the expres sion l|Chazalvie]L Il975t 
iLeroux-Hugon and Ghaza'iilll972() : 



lOnoda and Nagaosal 1200^ 



—int 



'AH 



cN 



(37) 



where N/p is the ratio of the density of ionized impu- 
rities and the carrier density, is the average distance 
between carriers in units of Bohr radius, / is the mean free 
path, and c ~ 10, varying slightly with scattering length. 
For the sh ort range s catter ing potential considered by 
Luttinger (|Lijttingeit ^58|) and Nozieres and Lewiner 
l)Nozieres and Lewinei , Il97^ . V{r) ~ VoS{f~ fl): 



—int 
'^AH 



'AH 



Tr\Vo\BOS{EF)kFl' 



(38) 



The estimate is a useful first guess at which mechanism 
dominates in different materials but one must keep in 
mind the simplicity of the models used to derive these 
expressions. 

In (Ga,Mn)As, Eq. (|37|l gives a ratio of intrin- 
sic to skew scattering contribution of the order of 
50 an d the intrinsic A HE is therefore likely to domi- 
nate l)Dietl et all l2003j) . Consistently, the experimen- 
tal Ra in metallic (Ga , Mn)As DMSs is proportional 
to ([Edmonds et all l2002bfl . Microscopic calcula- 
tions of the intrinsic low-temperature AHE conductiv- 
ity in (Ga,Mn)As were performed by taking into ac- 
count the Berry phase anomalous velocity term in the 
semiclassical Bol t zmann equation leading to Eq. H3()|) 
l|Jungwirth et all l2002c^) . and by applying t he fully 
quantu m mechanical Kubo formalism jjungwir th et all 
l20031j) . In both approaches, the KL kinetic-exchange 
model was used to obtain the hole band structure. 

In the Kubo formula, the dc Hall conductiv- 
ity for non-interacting quasip articles at zero external 
magnetic field is given by jjungwirth et alj l2003bt 
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The real part of Eq. (|39|) in the limit of zero scattering 
rate (T — > 0) can be written as. 



to2 J (27r)3 /^^•'"'.'^^ -'"^k) 
lui[{n' k\px\nk) {nk\py\n' k)] 

^ (E r-E ,r)2 ■ 

V nk n'k' 



(40) 



Realizing that the momentum matrix elements 
{n'k\pa\nk) = {ralK){n'k\dH{k)/dc,\nk), Eq. ^ 
can be shown to be equivalent to Eq. (|36|l . The 
advantage of the Kubo formalism is that it makes it 
possible to include finite lifetime broadening of the 
quasiparticles in the simulations. (See Section IVI.BI 
for the discussion of quasiparticle scattering rates in 
(Ga,Mn)As. Whether or not the lifetime broadening is 
included, the theoretical anomalous Hall conductivities 
are of order 10 fi"^ cm"^ for typical (Ga,Mn)As DMS 
parameters. On a quantitative level, non-zero F tends to 
enhance a ah at low Mn doping and suppresses a ah at 
high Mn concentrations where quasiparticle broadening 
due to disorder becomes comparable to the strength of 
the kinetic-exchange field. 

A systematic comparison between theoretical 
and experimental AHE data is shown in Fig. 1^ 
(|Jungwirth et all l2003b() . The results are plotted versus 
nominal Mn concentration x while other parameters of 
the measured samples arc listed in the figure legend. 
Experimental a ah values are indicated by filled squares 
and empty triangles correspond to theoretical data 
obtained for F = 0. Results shown in half-open triangles 
were obtained by solving the Kubo formula for a ah 
with non-zero F due to scattering off Mnca and As- 
antisites, or Mnca and Mnj impurities. The calculations 
explain much of the measured low-temperature AHE 
in metallic (Ga,Mn)As DMSs, especially so for the 
Mni compensation scenario. The largest quantitative 
discrepancy between theory and experiment is for the 
X = 8% material which can be partly explained by a non 
mean-field-like magnetic behavior of this specific, more 
disordered sample. 

In DMSs the AHE has played a key role in establishing 
ferromagnetism and in providing evidence for the hole- 
mediated coupling betw een Mn local moments ((Ohiiol 
Il998t lOhno et gZj 119921) . The dominance of the AHE 
in weak field measurements (see Fig. ^ allows the Hall 
resistance to serve also as a convenient proxy for the 
magnetization. On the other hand, the same property 
can obscure Hall measurements of the hole density. If 
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FIG. 36 Comparison between experimen tal and theoreti- 
cal an omalous Hall conductivities. From l),Tungwirtli et all 
I2003U) . 



the magnetization is not fully saturated at low fields, 
for example, then pah = o'ahPxx '^iH increase with in- 
creasing external field B through the dependence of (Tah 
on M{B), and hole densities derived from the slope of 
Pxy{B) will be too low. Accurate determination of hole 
densities in DMSs is essential, however, and the Hall 
effect is arguably the most common and accurate non- 
destructive tool for measuring the level of doping in semi- 
conductors. Hall experiments performed in high mag- 
netic fields to guarantee magnetization saturation seem 
a practical way for separating AHE contributions, es- 
pecially in s amples showing weak longitudinal m agne- 
toresistance (|Edmonds et ad l2002bt IOhnolll999(l . Hole 
density measurements performed using this technique as- 
sume that the Hall factor, th = {pxy{B) — pah) / {B / ep) 
with PAH = Pxy{B = 0), is close to 1 despite the 
multi-band spin-orbit coupled nature of hole dispersion 
in (Ga,Mn)As ferromagnets. In the following paragraphs 
we briefly review a theore tical analysis of this assumption 
l|Jungwirth et g^J . l2005h|) . 

Microscopic calculations in non-magnetic p-type GaAs 
with hole densities p ^ 10^^ — lO'^^ cm~'^ have shown that 
th can vary between 0.87 and 1.75, depending on doping, 
scattering mechanisms, and on details of the model use d 
for the GaAs valence band teim and Maierfeldi ITool . 
An estimate of the influence on th from the spin-splitting 
of the valence band and from the anomalous Hall term is 
based on the KL kinetic-exchange model description of 
the hole band structure. The Hall conductivity has been 
obtained by evaluating the Kubo formula at flnite mag- 
netic fields that includes both intra-band and inter-band 
transitions. The approach captures th e anomalous and 
ordinary Hall terms on equal footing (jjungwirth et aO . 

iooii). 

Many of the qualitative aspects of the numerical cal- 
culations, shown in Fig. 1371 can be explained using a 
simple model of a conductor with two parabolic uncou- 



1.5 - 



1- 



0.5 



1.5 - 



1- 
0.5- 



intra-band transitions 



— - A3Q=0, 1/x=50 meV 

1/1=1 SOmeV 

A;p=341meV, 1/T=50mev| 

1/1=1 50meV 



x=4%, 1/i=50meV 
1/1=1 SOmeV 
x=8%, 1/i=50me 
1/1=1 50m 



=341 me 



intra&inter-band transitions 



B(T) 



10 



FIG. 37 Theoretical HaU factors for p = 0.8 nm~=*; Ti/t = 
50 meV (dashed lines), Ti/t = 150 meV (solid lines). Top 
panels: only intra-band transitions are taken into account. 
Bottom panels: intra- and inter-band transitions are taken 
into account. Left panels: GaAs (x = 0); zero spin-orbit cou- 
pling (red lines), Aso = 341 meV (blue lines). Right panels: 
(Ga,Mn)As with Mnca concentration 4% (green lines), 8% 
(brown lines). pxy{B = 0) = in all panels except for the 
bottom left panel whe re pxv(B = 0) 7^ due t o the anoma- 
lous Hall effect. From Jjungwirth et aZ.ll2005b|) . 



pled bands. Note that the typical scattering rate in 
(Ga,Mn)As epilaycrs is h/r ^ 100 mcV and that the 
cyclotron energy at _B = 5 T is huj ^ 1 meV, i.e., the 
system is in the strong scattering limit, ujt <C 1. In this 
limit, the two band model gives resistivities: 
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where the indices 1 and 2 correspond to the 1st and 2nd 
band respectively, the total density p ^ pi + p2, and 
the zero-fleld conductivity ctq = e^rp/m* . Eq. |(1T|) sug- 
gests that in the strong scattering limit the multi-band 
nature of the hole states in (Ga,Mn)As should not re- 
sult in a strong longitudinal magnetoresistancc. This ob- 
servation is consistent with the measured weak depen- 
dence of Pxx on B for magnetic fields at which mag- 
netization in the (Ga,M n)As ferromagnet is saturated 
l)Edmonds et g/J . l2002b() . 

The simple two-band model also suggests that the Hall 
factor, rn, is larger than one in multi-band systems with 
different dispersions of individual bands. Indeed, for un- 
coupled valence bands, i.e. when accounting for intra- 
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band transitions only, the numerical Hall factors in the 
top panels of Fig. 1371 are larger than 1 and independent 
of T as also suggested by Eq. (|41|1 . The suppression of rn 
when spin-orbit coupling is turned on, shown in the same 
graphs, results partly from depopulation of the angular 
momentum j = 1/2 split-off bands. In addition to this 
"two-band model" like effect, the inter-Landau-level ma- 
trix elements are reduced due to spin-orbit coupling since 
the spinor part of the eigenfunctions now varies with the 
Landau level index. In ferromagnetic Gai_2;Mn2;As the 
bands are spin-split and higher bands depopulated as x 
increases. In terms of rjy, this effect competes with the 
increase of the inter-Landau-level matrix elements since 
the spinors are now more closely aligned within a band 
due to the exchange field produced by the polarized Mn 
moments. Increasing x can therefore lead to either an 
increase or a decrease in rn depending on other param- 
eters, such as the hole density (compare top right panels 
of Fig. EH). 

Inter-band transitions result in a more single-band like 
character of the system, i.e. rn is reduced, and the slope 
of the pxy{B) curve now depends more strongly on r. 
Although the AHE and OHE contributions to pxy can- 
not be simply decoupled, the comparison of numerical 
data in the four panels confirms the usual assumption 
that the AHE produces a field-independent off-set pro- 
portional to magnetization and p^^ . The comparison also 
suggests that after subtracting Pxy{B = 0), rn can be 
used to determine the hole density in (Ga,Mn)As with 
accuracy that is better than in non-magnetic GaAs with 
comparable hole densities. For typical hole and Mn den- 
sities in experimental (Ga,Mn)As epilayers the error of 
the Ha ll measurement of p is e stimated to be less than 
±20% (jJungwirth et a^J . l2005al) . 



D. Conductivity near and above Tc 

Typical Fermi temperatures, Tp = Ep/kB, in ferro- 
magnetic (Ga,Mn)As are much larger than the Curie 
temperature relegating direct Fermi distribution effects 
of finite temperature to a minor role in transport. The 
carrier-mediated nature of ferromagnetism implies, how- 
ever, strong indirect effects through the temperature de- 
pendence of the magnetization. A prime example is the 
AHE which from the early studies of (III,Mn)V DMSs has 
served as a p ractical tool to acc urately measure Curie 
temperatures (|Ohno et all 1X993 ). A rough estimate of 
Tc can be inferred also from the temperature dependent 
longitudinal resistivity which exhibits a shoulder in the 
more metallic (optimally annealed) samples and a peak in 
the less metallic (as-grown) mate r ials near the ferromag- 
netic transition jEd niond s et al. . l2002at iHavashi et all 
1997t iMatsukurae^ ali . Il998t IPotashnik et all l200l[ 
^nEsdr'er'aini997j) . An example of this behavior is 

% nom- 



shown in Fig. l38lfor a (Ga,Mn)As material with I 
inal Mn-doping ijPotashnik et a/J . 12001)) . 

The shoulder in p^x (T) has been qualitatively modeled 
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FIG. 38 Experimental resistivity of Gai_a;Mna; As for x = 
8% vs. tem perature for various annealing times. From 
iPotashnik et a/.^.200ll) . 



using the mean-field, KL kinetic-exchange Hamiltonian 
iven by Eq. I|21|l. Solut i ons to the Boltzmann equation 



■Hwang and Das Sarmal l2005t iLopez-Sancho and Brevl 
120031) are shown in Fig.|211 The temperature dependence 
of the longitudinal conductivity follows in this theory 
from variations in the parameters derived from the spin- 
polarized hole band structure {e.g. Fermi wavevector) 
and from variations in screening of impurity Coulomb 
potentials. 

The peak in resistance near Tc has been discussed in 
terms of scattering effects beyond the lowest order Born 
approximation and by using a network resistor model 
l)Timm et aZ.l.l2005|) . It has also been suggested that this 
transport anomaly in more highly resistive DMSs is a 
consequence of the change in lo calization length cau sed 
by the ferromagnetic transition ijZarand et aZ.! . 12005)) . 

Above the Curie temperature, measurements of pxx 
have been used to estimate the value of Jpd. Assuming 
scattering off the p — d exchange potential at randomly 
distributed paramagnetic MuQa impurities and parabolic 
hole bands, the corresponding contribution to the resis- 
tivity is approximated by. 
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(42) 

where X-L and X|| are t he transverse and longitudinal 
magnetic susceptib il ities ijDietlL Il994t IMatsukura et all 
l2002t lOmiva et all l2000|) . This theory overestimates 
critical scattering, particularly near the Curie tempera- 
ture where the susceptibility diverges. Far from the tran- 
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FIG. 39 KL kinetic-exchange model resistivity, normalized 
with respect to the value above Tc, vs. T for a Mn doping of 
5% and hole densities p = 0.2n m~^ and p = 0.2nm^'' From 
llLopez-Sancho and Brevl l2003t) . 
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FIG. 40 Temperature dependence of the resistivity for a 200- 
nm thick film of Gai-^Mn^^As with x = 0.053 in the high- 
temperature paramagnetic region. Solid squares and open 
circles show experi mental data and th e fit using Eq. 1421 . 
respectively. From iOmiva et aLl 120001) . 



sition on the paramagnetic side, however, fitting Eq. 
to experimental magnetoresistance data gives an esti- 
mate of the Jpd ijOmiva et alV l2000(l which is consistent 
with values inferred from spectroscopical measm'cmcnts 
l|Okabavashi et a/.llT998(l . 



VII. MAGNETO-OPTICS 

The prospects for new technologies based for example 
on materials in which the ferromagnetic transition can 
be controlled by light or on (iff ,Mn)V Faraday isolators 
monolithically integrated with existing semiconductor 
lasers, have motivated research in magneto-op tical prop- 
erties ofDMSs ( | Koshihara et a/.. 1997: Matsukura et al\. 
l2002HMunekata et aZ.l . ll997HSugano and Koiimall2000j) . 
Apart from these applied physics interests, ac probes 
have been used to study DMS materials by a wide 
range of experimental techniques. In Sections III.B.ll 
and IV.B~2l we mentioned x-ray spectroscopies (core- level 
photoemission and XMCD) used to characterize Mn 
id states and detect the sign and magnitude of the 
p — d exchange coupling. Dispersions of hole bands in 
DMSs have been studied by angl e-resolved photoemis- 
sion with ultraviolet excitations (|Asklund et all . l2002t 
IOkabavashi~^a/. . 2001. 200^ and infrared-to-ultraviolet 
spectroscopic ellipsometrv ijBurch et all l2004|) . Raman 
scattering induced by excitations in the visible range was 
used as an alte r native means of estimating hole densities 
llLimmer et all l2002t ISapeea et all . 120011: ISeong et all 
I2002r . Spectroscopic studies of isolated Mn(d^-l-hole) 
impurities in the infrared region provided key infor- 
mation on the valence of Mn in (Ga,Mn)As, as dis- 
cussed in Section lll.fjl and cyclotron resonance mea- 
surements were used to study hig hly Mn-doped DMS 
mater ials in this frequency range (iK hodaparast et 
2004: .Mitsumori et a/.l . l2003HSanders et aL . .2GGa) . The 
microwave EPR and FMR experiments, mentioned in 
Section III.BI and IV. CI have been invaluable for un- 
derstanding the nature of Mn in III-V hosts at low 
and high dopings, and for characterizing magnetocrys- 
talline anisotropics and magnetization dynamics in fer- 
romagnetic materials. In this section we re view studies 
of sq n ie of the magnet o -optic a l responses (lAndo et al. , 
I998t iBeschoten et all Il999t iHraboyskv ad 1200^ : 
Kimel et all l2005t iLang et aZ.l 120051 iMatsukura et al . 
2002t ISzczvtko et aZj . Il999a|) . Darticularly magnetic cir- 
cular dichroism ( MCD), in the visib l e range and in- 
frared absorption l|Bur ch et al\ l2005t iHirakawa et all 
20021 l200lUNagai et oi.Ll200lUSinglev et aiJ . l2003l 120021: 
Szczvtko et a/lll999a|) . 



A. Visible magnetic circular dichroism 

Optical absorption due to electron excitations across 
the band gap is a standard characterization technique in 
semiconductors. In (Ga,Mn)As, the absorption occurs 
in the visible range and the position of its edge on the 
frequency axis depends on the circular polarization of 
the incident light. Analysis of this magneto-optical effect 
provides information on the p—d exchange i nduced band 
splitting and on doping in the DMS material l)Diet]| . l2002t 
IMatsukura et a/J . l200^ . 

The schematic diagram in Fig. shows that for a 
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FIG. 41 Schematic diagrams of electron excitations across 
the band-gap induced by circularly polarized light absorption 
in undoped (left) and p-type (right) DMSs with spin-split 
bands. Only the heavy-hole band is shown for illustration. 
For a given circular polarization of the absorbed light, optical 
selection rules allow transitions between the heavy-hole va- 
lence band and the conduction band with one spin orientation 
only and therefore make it possible to spectrally resolve the 
band splitting. The corresponding MCD signal can change 
sign in the doped system due to the Moss-Burnstein effect, as 
illustrated in this diagram. 



given sign of the exchange coupUng, the order of absorp- 
tion edges corresponding to the two circular photon po- 
larizations can reverse in p-type materials, compared to 
systems with a completely filled valence band. Calcu- 
lations for (Ga,Mn)As that include this Moss-Burnstein 
effect were carrie d out using the me an-field KL kinetic- 
exchange model (|Dietl et aZ.1 . 1200111) . The resulting ab- 
sorptions or^ of t he (7^ circularly po l arized Ught, and 
MCD, defined as jSugano and Koiimall2000(l 



MCD = 



Re[cr^:r(tj)] ' 



(43) 



are shown in Fig. 021 As suggested in the cartoon, the 
sign of the MCD signal in (Ga.Mn)As is opposite to the 
one obtained in bulk (II,Mn)VI DMSs where the sense of 
band splittings is the same as in (Ga,Mn)As but Mn sub- 
stituting for the grou p-II element is an isovalent neutral 
impurity l)DietJ . ll994j) . (The Moss-Burnstein sign change 
in MCD was also observed i n co-d oped p-type (II,Mn)VI 
quantum well (Haurv a/.i Il997|) .'i 

Experimentally, the incorporation of several per cent of 
Mn in GaAs strongly enhances MCD , as sh own in Fig. 1431 
l|Ando et aZ.1 . 11998 iMatsukura et a;J . l2002() . and the sign 
of the signal near band-gap frequencies is consistent with 
the above theory w hich assumes antiferromagnetic p — d 
excha nge coupling llAndo et all IT 998t iBeschoten et all 
ITflfl^lDietl ail 12001 bUSzczvtko et aiJ.ll999aF 



B. Infrared absorption 



3.0 
„ 2.5 

3 

1. 2.0- 
o 

1.5. 

1 1.0. 

S 0.5 
IT 

o.oH 



-5 meV 
fl^ = -10meV 
flg = -20 meV 
30 meV 




1.5 1.6 1.7 1.8 

PHOTON ENERGY [eV] 



1.9 



2.0 



FIG. 42 Top panel: Theoretical absorption edge for two circu- 
lar polarizations in p-(Ga,Mn)As computed for spin-splitting 
parameter Bq = — 10 meV (corresponding to T = mean- 
field from local Mn moments of density x = 1.7%) and hole 
concentration 3.5 x 10^° cm~^. Inset shows how the Fermi 
sea of the holes reverses the relative positions of the edges 
corresponding to cr"*" and a~ polarizations in agreement with 
experimental findings. Bottom panel: Spectral dependence of 
magnetic circular dichroism in the optical range in (Ga,Mn) As 
computed for hole concentration 3.5 x lO'^^ cm^^ and vari- 
ous spin-splitting parameters Bq- The magnitudes of MCD 
at given Bp are n ormalized by its value at 1.78 eV. From 
llDietl et a;.ll200lH) . 



impurity ba i id exc itations jAlvarez and Dagottol l2003t 
ICraco et all l2003t |Hwang et o,ll l2002^ in the more in- 
sulating inatej^ls_and_djae to in t ra-valence-band e xei- 
tations ^Sinova et all l2003l 120021: lYang et all l2003(l in 

the more metallic systems, as illustrated schematically in 
Fig. ^1 The infrared absorption associated with substi- 
tutional Mn impurities is spectrally resolved from higher 
energy excitations to donor levels of the most common 
unintentional defects, such as the Mnj interstitials and 
AsGa antisites in (Ga,Mn)As, and therefore represents 
another valuable probe into intrinsic properties of these 
systems. Since infrared wavelengths are much larger than 
typical (sub-micron) DMS epilayer widths, the absorp- 
tion is related to the real part of the conductivity by 
ijSugano and Koiimair2000|) 



In (III,Mn)V DMSs, light absorption can occur at 
sub-band-gap frequencies due to valence-band to Mn 
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FIG. 43 Experimental MCD spectra of (a) undoped 
semi-insulating GaAs substrate and (b), (c) of epitaxial 
Gai_j;Mna;As films at T = 5 K and B = 1 T. The spectrum 
of GaAs is magnified ten times b ecause the signal is weaker 
than that of Gai_i,Mna;As. From llAndo et a^.. .1998,) . 



where Y and Yq arc the admittances of the substrate and 
free space, respectively. 

Model Hamiltonians (jSJl (see Section Fill. Dp combined 
with dynamical mean-field-theory or Monte Carlo sim- 
ulations were used to stud y the role of the irnpurity 
band in infra r ed ab sorption l)Alvarez and Dagottd . l20nit 
iHwang et all 12002}) . The impurity band forms in this 
theory when the strength of the model effective exchange 
interaction J is comparable to the width of the main 
band, characterized by the hopping parameter t. A non 
Drude peak is observed in the frequency-dependent con- 
ductivity, associated with transitions from the main band 
to the impurity band. The behavior is illustrated in 
Fig. together with the predicted temperature depen- 
dence of the absorption spectra obtained by the Monte 
Carlo technique. 



As discussed in Section IIILDI these 
model calculations are expected to apply to systems with 
strong p—d exchange coupling, like (Ga,Mn)P and possi- 
bly also to (III,Mn)V DMSs which are strongly compen- 
sated due to the presence of unintentional donor defects. 
(Impurity-band mediated ferromagnetism does not occur 
in uncompensated samples.) 

Theoretical infrared absorption spectra calculated us- 
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FIG. 44 Schematic diagram illustrating absorption due to 
electron excitations to the impurity band (left panel) and due 
to intra- valence-band excitations (right panel). 
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FIG. 45 Impurity-band model calculations of (j{ui) vs. ui for 
a 10x10 periodic system with 26 spins (a; ~ 0.25), J/t = 2.5, 
p = 0.3, and for different temperature s, as indicated. Inset : 
Drud e weight, D, vs. temperature T (lAlvarez and Daeottd . 
I2003D. 



ing the k • p model for (Ga,Mn)As DMSs with delocal- 
ized holes in the semiconductor valence band, plotted in 
Fig. 1461 show a similar non-Drude characteristics with a 
peak near the excitation energy of 220 meV. The underly- 
ing physics is qualitatively different, however, as the peak 
in these KL kinetic-exchange model calculations orig- 
inates from heavy-hole to ligh t-hole intra-valence-band 
transitions ijSinova et all |2002() . These results were ob- 
tained by evaluating the Kubo formula for ac conductiv- 
ity assuming non-interacting holes and modeling disorder 
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within the first order Born approximation (see Eq. (|33|l in 
Section rVl.A|l . In Fig.^lwe show theoretical predictions 
of exact diagonalization studies based on the KL kinetic- 
exchange Hamihonian but treating disorder effects ex- 
actly in a finite size system. The results correct for the 
overestimated dc conductivity in the former model, which 
is a quantitative deficiency of the Born approximation as 
already mentioned in Section [VI.AI At finite frequencies, 
the theoretical absorption in these metallic (Ga,Mn)As 
DMSs is almost insensitive to the way disorder is treated 
in the simulations, as see from Figs. and IT7I 
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FIG. 46 KL kinetic-exchange model calculations of infrared 
conductivity Re[(T(ci;)] and absorption coefiicient a(tj) for car- 
rier densities from p = 0.2 to O.Snm"'' i n the direction ind i- 
cated by the arrow, for Gao.95Mno.05As l)Sinova et aLLl2002ri . 
Disorder is treated within the first order Born approximation. 
The dot-dashed line is the experime ntal absorption curve for 
a sample with 4% Mn doping from ijHirakawa et al\. 120021) . 

Experimental infrared absorption studies in ferro- 
magnetic (Ga,Mn)As epilayers exhibit several com- 
mon features summarized i n Fig . 1481 jBurch et all . 
2005t iHirakawa et all I200I l200ft iNagai al\ l200l[ 
Singlev et alV l2003l l2002t ISzczvtko et all Il999a^. Fer- 
romagnetic materials {x = 5.2% curves in Fig. I48|l show 
a non-Drude behavior in which the conductivity increases 
with increasing frequency in the interval between and 
220 mcV, a broad absorption peak near 220 — 260 meV 
that becomes stronger as the sample is cooled, and a fea- 
tureless absorption up to approximately 1 eV. As seen 
in Fig. 0H| the peak is absent in the reference, LT-MBE- 
grown GaAs sample confirming that the infrared absorp- 
tion in ferromagnetic (Ga,Mn)As is related to changes in 
the band structure near the Fermi energy induced by Mn 
impurities. 

The presence of a finite-frequency peak in both im- 
purity band and KL kinetic-exchange models for the in- 
frared conductivity leads to an ambiguity in the inter- 
pretation of existing data which has for the most part 
been taken in as grown, presumably heavily compensated 
material. The metallic behavior of the x = 5.2% ma- 
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FIG. 47 Infrared absorption and conductivity of a metallic 
(Ga,Mn)As computed using the KL kinetic-exchange model 
and exact diagonalization technique in a finite size disordered 
syste m. Here p = 0.33 nm~^, UMn ~ 1 nm~^ {x « 4.5%). 
From (lYang ei adl2003D . 

terial below Tc, seen in the lower panel of Fig. ^HJ fa- 
vors the inter-valence-band absorption scenario. On the 
other hand, the large compensation likely present in as- 
grown low-Tc (Ga,Mn)As suggests that many holes may 
be strongly localized and that both absorption mecha- 
nisms may contribute to the measured absorption peak. 
Experiments in a series of samples interpolating between 
as grown and optimally annealed limits, analogous to the 
resistance-mo nitored annealing studies ijEdmonds et alV 
I2004ail2002al) . should enable a clear interpretation of in- 
frared absorption spectra in (Ga,Mn)As DMSs. These 
studies will hint towards necessary refinements of the 
simplified theories used so far, 6.17. inclusion of the energy 
dependence of Jpd and a more quantitative theory of the 
impurity band model. The support for either scenario by 
these experiments has to be considered in conjunction 
with other available data in a self-consistent picture, e.g. 
in an impurity band picture Tc is predicted to approach 
zero as the system reaches zero compensation whereas 
the KL kinetic-exchange model has an opposite trend. 

VIM. DISCUSSION 

A. Magnetic interactions in systems with coupled local and 
itinerant moments 

Systems with local moments coupled to itinerant elec- 
trons are common in condensed matter physics and ex- 
hibit a wide variety of behaviors. Ferromagnets are far 
from the most common low-temperature states. For that 
reason it is useful to ask how (III,Mn)V materials, and 
(Ga,Mn)As with its robust ferromagnetic order in par- 
ticular, fit in this larger context. This general qualitative 
analysis can help to identify some of the key factors that 
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FIG. 48 Top panel: Real part of the conductivity derived 
from measured transmission spectra for a paramagnetic x — 
1.7% sample (dashed line), a ferromagnetic x = 5.2% sample 
(thick solid line), and the LT GaAs film (thin black line). 
Bottom panel; Temperature dependence of conductivity for 
X = 5.2% sample on a log scale for T > Tc and T <Tc ■ From 
llSinglev et ffl/.tl2003D . 



might limit the strength of ferromagnetic interactions in 
the highly doped and strongly p ^ d coupled DMS ferro- 
magnets for which mean-field theory predicts the highest 
Curie temperatures. 

An important class of materials that has been very ex- 
tensively studied is heavy fermions. in which f-clectron 
local moments are exchange coupled to band electrons 
fStewart , 198^ . Kondo lattice models, which are 
believed to qualitatively describe heavy fermion sys- 
tems, assume that local moments exist at each lat- 
tice site. Models of DMS systems which make a vir- 
tual crystal approximation (see Section lil.B.ll ind lV.A.l|) 
place moments on all lattice sites and are therefore 



Kondo lattice models, often with specific details that 
attempt to capture some of the peculiarities of specific 
DMS materials. Theories of Kondo lattice model of- 
ten start from the comparison of the R KKY (see Sec- 
tionlll. All and Kondo temperature scales l)Degiorglll999t 
iDoniachl [l977t iTsunetsueu et adll997(l . The character- 
istic RKKY temperature refers to the strength of interac- 
tions between local moments mediated by a weakly dis- 
turbed carrier system and is proportional to the mean- 
field Tc given by Eq. H19|l . The Kondo scale refers to 
the temperature below which strong correlations are es- 
tablished between an isolated local moment and the car- 
rier system with which it intera cts. Standard scale es- 
timate formulas ijDoniachl Il977|) applied to the case of 
DMS ferromagnets imply that the Kondo scale is larger 
than the RKKY scale when the mean-field exchange cou- 
pling, SNMnJpd, is larger than the Fermi energy Ep of 
the hole system, in other words in the strong coupling 
regime. The Kondo scale falls rapidly to small values at 
weaker coupling. (In heavy-fermion materials the Kondo 
temperature scale is larger than the RKKY temperature 
scale.) 

Optimally annealed metallic (Ga,Mn)As materials 
are on the weak-coupling side of this boundary, but 
(Ga,Mn)N and possibly (Ga,Mn)P may be starting to 
reach toward the strong coupling limit (if the simple 
S=5/2 local moment model still applies in these mate- 
rials). As the strong coupling limit is approached, quan- 
tum fluctuations (see discussion below Eq. (|25|) ) will play 
a greater role, reducing the saturation moment per Mn 
and eventually driving down the ferromagnetic transition 
temperature. When the Kondo scale is much larger than 
the RKKY scale, the local moments are screened out by 
strongly correlated band electron spin-fluctuations and 
effectively disappear before they have the opportunity to 
couple. 

On the weak coupling side, RKKY interactions in 
Kondo lattice models tend to lead to ferromagnetism only 
when the number of itinerant electrons per moment is 
small, i.e., only when at least the near-neighbor RKKY 
interaction is ferromagnetic ijTsunetsueu et aZl Il997j) . 
One of the surprising features of (Ga,Mn)As is the prop- 
erty that ferromagnetism still occurs when the number of 
itinerant electrons per moment is 1. As we have men- 
tioned in Section IV. A. 21 this property follows from the 
specific multi-band electronic structure and spin-orbit 
coupling at the top of the valence band. Nevertheless, 
frustrating antiferromagnetic RKKY interactions, and 
exchange interactions that promote non-collinear mag- 
netic states (see Section IV.BII , will eventually become 
important for sufficiently large carrier densities. As this 
regime is approached from the ferromagnetic side, the 
transition temperature will be suppressed. These ten- 
dencies are summarized schematically in Fig. 1491 

Similar considerations apply in assessing the robust- 
ness of ferromagnetism in the ordered state as char- 
act erized by the spin-stiffness micromagnetic parame- 
ter (|Konig et a;J . l2001bll2003lSchliemann et a/.ll200ia(l . 
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FIG. 49 Top panel: A schematic qualitative diagram illus- 
trating the requirements for ferromagnetic ordering in systems 
with coupled local and itinerant moments. The y-axis repre- 
sents the strength of the exchange coupling relative to the itin- 
erant system Fermi energy and the x-axis the ratio between 
carrier and local moment densities. In the weak coupling 
regime the interaction between local moments is described by 
RKKY theory in which polarization of band electrons (holes) 
due to the interaction with the local moment at one site is 
propagated to neighboring sites. The mechanism corresponds 
to out-of-phase Friedel oscillations of carriers scattering off 
an attractive impurity potential for one carrier spin orienta- 
tion (maximum density at the impurity site), and a repulsive 
potential for the opposite carrier spin polarization (minimum 
density at the impurity site). With increasing carrier den- 
sity, the ferromagnetic RKKY state gradually approaches a 
regime of frustrated RKKY interactions. The RKKY picture 
does not apply to strong couplings. Here the tendency to 
order ferromagnetically is weakened by correlated (flip-flop) 
quantum fluctuations of the interacting local and itinerant 
moments making the moment disappear eventually. This is 
the so-called Kondo singlet regime. Bottom panel: The ro- 
bustness of the ferromagnetic state is viewed from the ferro- 
magnetic low-temperature side. With increasing hole density, 
the ferromagnetically stiff mean-field state will start to suf- 
fer from the occurrence of lower energy spin-wave excitations 
until finally the energy of spin-waves will become negative sig- 
naling the instability of the ferromagnetic state. In the strong 
coupling regime, hole spins will tend to align locally and in- 
stantaneously in opposition to the fluctuating Mn moment 
orientations. Since the hole Fermi energy is relatively small 
in this regime, the kinetic energy cost of these fluctuations 
is small, resulting in soft spin-wave excitations for the mag- 
netic system. The long-range ferromagnetic order disappears 
eventually at temperatures much smaller than the mean-fleld 
Tc, while short-range order may still exists above the Curie 
temperature. 
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We have used this approach in Section IV. A. II when 
analyzing the limitations of the mean-field theory in 
(Ga,Mn)As. Starting from the ferromagnetic state, 
long-wavelength spin-orientation modulation will tend 
to lower the energy of some Mn-Mn interactions for 
sufficiently high carrier densities. The spin-stiffness 
will weaken as the frustrated magnetism regime is ap- 
proached, until finally the energy of ferromagnetic spin- 
waves will become negative signaling the in stability of 
this state (|Schliemann and MacDonaldLl200^ . Similarly 
for sufficiently strong coupling, the band system will be 
(nearly when spin-orbit coupling is included) fully spin- 
polarized and the cost of spin-orientation spatial varia- 
tion will be borne mainly by the kinetic energy of the hole 
system and no longer increase with exchange coupling. 
As shown schematically in the lower panel of Fig. 1491 
hole spins in this regime are locally antiferromagnetically 
locked to the fluctuating Mn moment orientations. For 
relatively small hole Fermi energies the kinetic energy 
cost of these fluctuations is small, resultin g in soft spin- 
wave e xcitations of the magnetic system (|Konig et all 
l2001tJl . In this regime long-range ferromagnetic order 
disappears at temperatures smaller than the mean-field 
Tc, i.e., short-range order still exists above the Curie tem- 
perature. 

The white bottom-left area in the panels of Fig. 1^ 
qualitatively depicts the parameter range in which the 
ferromagnetic RKKY mean-field state applies. Here Tc 
increases with the carrier density, local moment den- 
sity and (quadratically) with the strength of the p ~ d 
exchange coupling. For a fixed ratio of p/Ni\jn, in- 
creasing Nmti corresponds to moving only slowly (as 

1/3 

^Mn) upwards in the diagrams. This may explain why 
(Ga,Mn)As materials with larger Mn doping and with 
similar hole compensations as in the ferromagnetic low 
local moment density systems do not show any marked 
weakening of the ferromagnetic state. With p/Nmu still 
fixed, attempts to increase Tc in (III,Mn)V DMSs by in- 
creasing the Jpd constant, in e.g. ternary host alloys 
of Ga(As,P), might at some point reach the boundary 
of the soft spin-wave (Kondo screened) state. Similarly, 
the Kondo lattice model allows only a limited space for 
enhancing the robustness of the ferromagnetic state by 
tweaking the carrier and local moment densities indepen- 
dently. In this case moving horizontally from the bound- 
ary of the frustrated RKKY (ferromagnetically unsta- 
ble) state is accompanied by approaching vertically the 
soft spin-wave (Kondo screened) regime, and vice versa. 
Viewed from the opposite perspective, however, it is as- 
tonishing that a window in this parameter space has been 
found by the material research of (III,Mn)V compounds 
for robust DMS ferromagnets with Curie temperatures 
close to 200 K. The diagrams do not imply any general 
physical mechanism that limits Tc in these materials be- 
low room temperature. 

Our remarks on the cartoons in Fig. 0^1 refer to the 
properties of Kondo lattice models which have been ad- 
justed to reflect peculiarities of the zinc-blende semicon- 
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ductor valence band. We have so far neglected the im- 
portance of disorder and of Coulomb scattering in DMSs, 
and these can modify some parts of the simple quali- 
tative picture depicted in Fig 1491 This is particularly 
true in the very low-density isolated Mn limit, i.e., very 
strong exchange coupling limit. Because of Coulomb at- 
traction between valence band holes and the charged Mn 
ion which carries the local moment, a total angular mo- 
mentum F = I isolated bound state is created, as we 
have explained in Section 111. B. II instead of the strongly 
correlated Kondo singlet. The importance of Coulomb 
interactions and disorder are lessened by screening and 
Pauli exclusion principle effects when both the Mn den- 
sity and the carrier density are high. 



B. Ferromagnetism in the diluted magnetic semiconductor 
family 

In this section we will narrow down the discussion of 
magnetism in local moment systems to semiconducting 
compounds, focusing on the phenomenology of ferromag- 
netic DMSs other than the (lll,Mn)V materials. Almost 
any semiconducting or insulating compound that con- 
tains elements with partially filled d or f shells (local 
moments) will order magnetically at a sufficiently low 
temperature. Semiconductors and insulators with high 
density of magnetic moments usually order antiferromag- 
netically however, although ferromagnetism does occur 
in some cases. A famous example of a ferromagnetic sys- 
tem that can be regarded as a doped semiconductor is 
provided by the manganite family {e.g. Lai-ajSrj^MnOa) 
whose ferromagnetism is favored by the double-exchange 
mechanism and occurs over a wide range of transition 
temperatures from below 100 K to nearly 400 K. The 
onset of magnetic order in these systems is accompanied 
by a very large i ncrease in conduc tivity. For a review, 
see for example (|Coev et all Il99^ . Other well known 
dense moment (of order one moment per atom) ferro- 
magnetic semiconductors with strong exchange interac- 
tion between it inerant and local spins include Eu- and Cr- 
chalc o eenides jBaltzer et all 1 19661 Kasuva and Yanas3 . 
Il968t iMaueer and Godartl 119861 Ivan Staoelel Il982|) . 
such as rock-salt EuO and spinel CdCr2Se4 with Curie 
temperatures 70 K and 130 K, respectively. 

DMS systems in which magnetic atoms are introduced 
as impurities have moments on only a small fraction of all 
atomic sites. The mechanisms that control magnetic or- 
der are therefore necessarily associated with the proper- 
ties of these impurities. The coupling between moments 
will generally depend on the locations of the dilute mo- 
ments in the host lattice, on the doping properties of 
the magnetic impurities, and on other dopants and de- 
fects present in the material. It seems plausible therefore 
that when DMS systems are ferromagnetic, their mag- 
netic and magnetotransport properties will be more sen- 
sitive to engineerable material properties. 

This review has concentrated on (Ga,Mn)As and re- 



lated materials in which, as we have explained, substi- 
tutional Mn acts both as an acceptor and as a source of 
local moments. Ferromagnetism is carrier mediated and 
it has been demonstrated that it persists to surprisingly 
high temperatures. More may be achieved in the future 
by tweaking these materials. On the other hand there 
is a vast array of alternate DMS materials that could 
be contemplated. Research to date has only scratched 
the surface of the volume of possibilities - we are truly 
still at the beginning of the road in studying diluted mo- 
ment magnetism in semiconductors. Each system brings 
its new challenges. The interpretation of simple magnetic 
and transport characterization measurements is often not 
immediately obvious, in particular because of the possi- 
bility that the moments will segregate into crystallites 
of one of a variety of available dense moment minority 
phases which are often thermodynamically more stable. 
In addition, magnetic properties will certainly depend in 
general not only on the dilute moment density, which 
normally is well controlled and variable, but also on the 
partitioning of local moments among many available sites 
in the host crystal which is not always known and is 
usually much harder to control. The search for promis- 
ing DMS materials would be simplified if ab initio DFT 
methods had reliable predictive power. Unfortunately 
this luxury appears to be absent in many cases because 
of extreme sensitivity of magnetic properties to details 
of the electronic structure and because of strong correla- 
tion effects that are often present in these systems. We 
mention briefly in the following paragraphs some of the 
other classes of diluted magnetic semiconductors which 
have been studied. 

The class of DMS ferromagnets that is closest to 
(III,Mn)V materials is (II,Mn)VI compounds co-doped 
with group V element acceptors. Examples include p- 
Zn,Mn)Te:N llFerrand et aZ.l . l200l|) and p-(Be,Mn)Te:N 
Hansen et all l200l|) .~ These materials differ from 
(Ga,Mn)As mainly because the local moments and holes 
are provided by different types of impurities and can 
be controlled independently. Although the physics be- 
hind ferromagnetism seems to be very similar in the two 
classes of materials, the highest ferromagnetic transition 
temperatures that have been achieved are much smaller 
in the case of co-doped (II,Mn)VI materials, ~ 2K rather 
than ^ 200K. The difference is explained partly by dif- 
ficulty in achieving the same extremely high hole doping 
(^ lO^^cm"'^) in (II,Mn)VI that has been achieved in 
(Ga,Mn)As and partly by a favorable interplay between 
electrostatic and magnetic effects in the (III,Mn)V ma- 
terials. In (III,Mn)V materials, unlike (II,Mn)VI ma- 
terials, the Mn moment is charged and attracts holes. 
The tendency of holes to have a higher density near 
Mn sites tends to increase the effective strength of the 
p — d exchange interaction. This effect is magnified 
when two Mn moments are on neighboring cation sites. 
In (Ga,Mn)As the interaction between Mn moments on 
neighboring cation positions is ferromagnetic, compared 
to the strongly antiferromagnetic interactions seen in 
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(II,Mn)VI materials. (For a discussion of the inter- 
play between ele c trosta tic and magnetic interactions see 
l|Sliwa and Dietll . l200,'ij) .) In p-doped (II,Mn)VI materi- 
als, competition between antiferromagnetic near neigh- 
bor interactions and the longer range carrier-mediated 
ferromagnetic interactions suppresses the magnetic or- 
dering temperature. This competition apparently docs 
not occur in (III,Mn)V ferromagncts with large hole den- 
sities. 

(Zn,Mn)0 is an i nteresting II- VI counterpart of the 
nitride III-V DMS l)Liu et all l2005(l . With advances 
in oxide growth techniques (Zn,Mn)0, can be consid- 
ered to be much like other (II,Mn)VI DMS materials 
and its inves tigation was origin ally motivated by theo- 
retical work l)Dietl et all l2000|) that extrapolated from 
experience with (III,Mn)V DMS ferromagnets and pre- 
dicted large Tc's. Studies of this material have pro- 
vided clear evidence of strong p — d exchange but, so 
far, have led to inconsistent conclusio ns about the oc- 
currences of l ong range magnetic order (iFiikumura et 
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Tctrahcdral DMS materials doped with transition 
metal atoms other than Mn have shown promising re- 
sults. For instance (Zn,Cr)Te is apparently homoge- 
neous, has th e required c o upling between local mom ents 
and carriers l)Mac et all 119961 ISaito et all |2003|) and 
Curie temperature as high as 300 K has been reported 
for this material. It may be, though, that the fcrromag- 
nctism is due to superexchange interactions rather than 
being carrier mediated since it occurs at very small ratios 
of the carrier density to the moment density (|Saito et all 
l2002j) . Another interesting material with Cr moments 
is (Ga,Cr)N whic h exhibits ferromagnetism at ^ 900K 
ijLiu et aZ.l . l2004a|) . The question still at issue in this ma- 
terial is the possible role of dense moment precipitates. 

Traditional groups of ferromagnetic DMSs also include 
riV,Mn)VI solid solutions wit h the rock-salt structure 
fe :Kgcnkamp et al , 1993; Storv et all [r992. .198^ . Al- 
though the band structures of IV- VI and III-V semicon- 
ductors are quite different, these DMS ferromagnets {e.g. 
(Pb,Sn,Mn)Te) appear to have a carrier-mediated mech- 
anism quite similar to that of (Ga,Mn)As. Holes with 
densities up to 10^^ cm~'^ are supphed in these materials 
by cation vacancies, rather than by Mn substitution for 
the divalent cations. The repor ted Curie temperatures in 
(Pb,Sn,Mn)Tc arc below 40 K (|Lazarczvk et all\l99T^ . 

DMS ferromagnetism with Si or Ge as the host 
semiconductor is obviously attractive because of their 
greater compatibility with existing silicon based tech- 
nology. In Si; Mn impurities favor interstitial po- 
sition which significantly complicates the synthesis of 
a uniform DMS system. Mn in Ge, on the other 
hand, is a substitutional impurity and ferromagnetism 
has be en report ed in MBE grown Ge.;,M ni_,; thin film 
DMSs l|Li et all l2005t iPark et all 120021) . Careful stud- 
ies of Gc^Mni_^ (|Li et a;J . l2005|r have demonstrated that 
slow low-temperature growth is required to avoid the 



formation of thermodynamically stable dense-moment 
ferromagnetic precipitates; it is likely that the high- 
temperature ferromagnetism sometimes found in these 
materials is du e to precipitates. The latest studies 
()Li et all l2005j) appear to indicate that true long-range 
order in Gca^Mni-a, emerges only at low temperatures 
^ 12K and that weak coupling between remote mo- 
ments is mediated by holes which arc tightly bound 
to Mn acceptors. Further work is necessary to de- 
termine whether this picture of magnetism, reminis- 
cent of the polaronic physics discussed in the context 
of (Ga,Mn)P or low-carrier-density (III,Mn)V system s 
()Kaminski and Das Sarmail2003tlScarDulla et a/J . l200,^ . 
applies to Ge^Mui-a; DMSs. 

The possible presence of dense- moment thermodynam- 
ically stable precipitates has also confused studies of ox- 
ide semiconductor DMS systems. More consistent ev- 
idence of above room temperature ferromagnetism has 
been reported in Co-doped Ti02 although the origin 
of ferroma gnetism in this material is still under debate 
see, e.g.. llMatsu motq et aZj. l200 ll) and re cent reviews 
Fukumura et a;J . l2005HPrellier et ai\ . \200^ ). Mn-doped 
indium-tin-oxide (ITO) is another promising candidate 
for a transparent ferromagnetic semiconductor which 
could be easily integrated into magneto-optical devices. 
Particularly encouraging is the observation of a large 
anomalous Hall effect showing that charge transport and 
magnetism a r e inti mately connected in this oxide DMS 
l)PhiliD et aZ.l.l2004tl . 

Other interesting related materials are the Mn-doped 
II-IV-V2 chalcopyrites surveye d theoretically in a firs t 
principles calculation study ijErwin and Zuti j . l2004|) . 
Three of these compounds, CdGeP2, ZnGeP2 and 
ZnSnAs2, have shown ferromagnetism experimentally. 
The origin of this ferromagnetic behavior has not been 
explored extensively yet. 

Finally we mention recent observations of ferromag- 
netic order up to « 20 K in a layered semicon ductor 
Sb2Te3 doped with V or Cr llDvck et all l200,1 [2002). 
These highly anisotropic materials combine DMS behav- 
ior with strong thermoelectric effects. The character of 
the ferromagnetic coupling in these compounds is unclear 
at present. 



IX. SUMMARY 

This article is a review of theoretical progress that has 
been achieved in understanding ferromagnetism and re- 
lated electronic properties in (III,Mn)V DMSs. The ma- 
terials we have focused on have randomly located Mn(c?^) 
local moments which interact via approximately isotropic 
exchange interactions with itinerant carriers in the semi- 
conductor valence band. (Ga,Mn)As is by far the most 
thoroughly studied and the best understood system in 
this class. Some (III,Mn)V materials may exhibit fluc- 
tuations in the Mn valence between Mn(c?^) and Mn(c?^) 
configurations or have dominant Mn(c?^) character, possi- 
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bly (Ga,Mn)N for example. Magnetic and other proper- 
ties of materials in the latter class will differ qualitatively 
from those of (Ga,Mn)As and this review makes no at- 
tempt to discuss the theory that would describe them. 
When we refer to (III,Mn)V ferromagnetism below, it 
should be understood that any materials that prove to 
be in the latter class are excluded. 

Interest in DMS ferromagnetism is motivated by the 
vision that it should be possible to engineer systems that 
combine many of the technologically useful features of 
ferromagnetic and semiconducting materials. This goal 
has been achieved to an impressive degree in (III,Mn)V 
DMSs, and further progress can be anticipated in the 
future. The goal of high temperature semiconductor fer- 
romagnetism flies in the face of fundamental physical lim- 
its, and the fact that so much progress has nevertheless 
been achieved, is due to a serendipitous combination of 
attributes of (III,Mn)V materials. We have reserved the 
term ferromagnetic semiconductor for materials in which 
the coupling between local moments is mediated by car- 
riers in the host semiconductor valence or conduction 
bands. Then magnetic properties can be adjusted over 
a broad range simply by modifying the carrier system 
by doping, photo-doping, gating, heterojunction band- 
structure engineering, or any technique that can be used 
to alter other semiconductor electronic properties. Most 
of these tuning knobs have already been established in 
(III,Mn)V ferromagnets. 

The progress that has been made in achieving 
(Ga,Mn)As ferromagnetism and in understanding its 
phenomenology has a few lessons. The analysis of any 
DMS should start with understanding the properties of 
isolated defects associated with the magnetic element. 
In the case of (Ga,Mn)As materials the desirable mag- 
netic defect is substitutional MuQa, because Mn then 
both introduces a local moment and acts as an accep- 
tor. The holes doped in the system by Mnca impurities 
provide the glue that couples the moments together. Un- 
derstanding the role of other defects that arc present in 
real material is also crucial. Substitutional Mnj is par- 
ticularly important in (Ga,Mn)As because it reduces the 
number of free moments and reduces the density of the 
hole gas that mediates ferromagnetism. Learning how to 
remove defects that are detrimental to strong magnetic 
order is key to creating useful materials. These two steps 
have been largely achieved in (Ga,Mn)As. There is ev- 
ery reason to believe that if the same progress can be 
made in other DMS materials, some will be even more 
magnetically robust. 

We have reviewed in this article a number of theoreti- 
cal approaches that shed light on what controls key prop- 
erties of ferromagnetic (III,Mn)V semiconductors. First 
principles electronic structure calculations give a good 
overview of fundamental material trends across the se- 
ries and explain many of the structural characteristics of 
these alloys. Semi-phenomenological microscopic tight- 
binding models provide a convenient way to use exper- 
imental information to improve the quantitative accu- 



racy of the description. Another phenomenological de- 
scription that successfully models magnetic, magneto- 
transport, and magneto-optical properties is a single pa- 
rameter theory that adopts a k • p description for the 
host semiconductor valence bands, and assumes that the 
exchange interaction between local moments and band 
electrons is short-ranged and isotropic. The single ex- 
change parameter that appears in this theory can be de- 
termined by fitting to known properties of an isolated 
Mn local moment, leading to parameter- free predictions 
for ferromagnetism. Qualitative models which focus on 
what kind of physics can occur generically for randomly 
located local moments that are exchange coupled to ei- 
ther localized or itinerant band electrons, also provide 
useful insights for interpreting experiments. 

The most important properties of (III,Mn)V materials 
are their Curie temperature and ferromagnetic moment 
which reflect both the strength of the coupling between 
Mn local spins and its range. The highest ferromagnetic 
transition temperatures in (Ga,Mn)As epilayers have so 
far been achieved with substitutional Mnca fractions in 
the neighborhood of 5% by post-growth annealing which 
eliminates most of the interstitial Mn ions. Achieving Tc 
values close to room temperature in (Ga,Mn)As, which 
is expected to occur for 10% MuQa doping, appears to 
be essentially a material growth issue, albeit a very chal- 
lenging one. In optimally annealed samples, experiment 
and theoretical considerations indicate that the Mn-Mn 
exchange interactions are sufficiently long range to pro- 
duce a magnetic state that is nearly collinear and insensi- 
tive to micro-realization of the Mnca spatial distribution. 
Magnetization and Curie temperatures in these systems 
are well described by mean-field theory. 

The magnetic and transport properties of high-quality 
(Ga,Mn)As materials are those of a low-moment-density, 
low-carrier-density metallic ferromagnet, with a few spe- 
cial twists. Because of the strong spin-orbit interactions 
in the valence band, metallic (Ga,Mn)As shows a large 
anomalous contribution to the Hall effect and the source 
of magnetic and transport anisotropics is the more itin- 
erant electrons, unlike the transition metal case in which 
anisotropics originate primarily in d-clcctron spin-orbit 
interactions. The small moment densities lead to a large 
magnetic hardness reflected in a single-domain-like be- 
havior of many (Ga,Mn)As thin fllms. They also ex- 
plain a large part of the orders of magnitude reduction 
in the current densities required for transport manipu- 
lation of the magnetic state through spin-momentum- 
transfer effects. The low-carrier density of the itiner- 
ant holes responsible for magnetic coupling means that 
they are concentrated around a particular portion of the 
Brillouin-zone in the valence band which has a large os- 
cillator strength for optical transitions to the conduction 
band. This property opens up opportunities for optical 
manipulation of the magnetic state that do not exist in 
transition metal ferromagnets and have not yet been fully 
explored. The research which we have reviewed here that 
is aimed at an understanding of the optical properties of 
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(Ga,Mn)As fcrromagnets is still incomplete, particularly 
for the most interesting ideal annealed materials, and will 
be important in setting the ground work for the explo- 
ration of new effects. 

These conclusions do not necessarily apply to all 
(III,Mn)V ferromagnets. For example, material trends 
suggest that wider band-gap hosts would have stronger 
exchange scattering that would lower the conductivity, 
shorten the range of Mn-Mn exchange interactions, and 
increase the importance of quantum fluctuations in Mn 
and band hole spin orientations. This could eventually 
lead to Curie temperatures significantly below the mean- 
field estimates. In the opposite limit, when the exchange 
interaction is weak enough to be treated perturbatively, 
sign variations in the RKKY Mn-Mn interaction are ex- 
pected to lead to frustration and weaken ferromagnetism 
at large carrier densities. One of the important miracles 
of (Ga,Mn)As ferromagnetism is that this effect is much 
weaker than would naively be expected because of the 
complex valence band structure. In (Ga,Mn)As, incipi- 
ent frustration that limits magnetic stiffness may be re- 
sponsible for the weak dependence of ferromagnetic tran- 
sition temperature on carrier density. This property of 
(Ga,Mn)As suggests that little progress on the Tc front is 
likely to be gained by non-magnetic acceptor co-doping. 

Disorder is an inevitable part of the physics of all 
DMS ferromagnets because of the random substitution 
of elements possessing moments for host semiconduc- 
tor elements. Even in metallic, ideal annealed sam- 
ples that have only substitutional MuQa impurities, 
randomness in the Mn microstructure leads to both 
Coulomb and spin-dependent exchange potential scat- 
tering. For (Ga,Mn)As, Coulomb scattering dominates 
over exchange potential scattering, limiting the conduc- 
tivity to ~ 100— 1000 r2~^cm~^. Frustration and disorder 
are certainly very important near the onset of ferromag- 
netism at low Mn density, where the network of exchange 
interactions that lead to long range order is still tenuous. 
Studies of well characterized materials with low Mn frac- 
tion near the metal-insulator transition are now possible 
because of progress in understanding the role of defects 
and are likely to exhibit complex interplay between dis- 
order, and Coulomb and exchange interactions. 
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